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NOVEL PROTEINS WITH INSUUN-UKE ACTIVITY 
USEFUL IN THE TREATMENT OF DIABETES 

This application Is a continuing application of U.S.S.N. 60/134,930. filed May 19, 1999. 

RELD OF THE INVENTION 

5 The invention relates to novel insulin activity OA) proteins and nucleic adds. The invention further 
relates to the use of the lA proteins in the treatment of insulin-related disorders such as type 1 
diabetes and type 2 diabetes. 

BACKGROUND OF THE INVENTION 

Insulin is a honmone that plays a mayor role in the regulation of growth and metabolism jq vertebrates. 

10 A deficiency of Insulin is the most important factor in diabetic disease states. Absence of insulin leads 
to severe metabolic disorders resulting from the failure to nonmally m^abolize cartx)hydrate, fat and 
proteins at a normal rate. These disorders include, for example diabetes mellitus (DM), a complex 
chronic metabolic disorder. Diabetes mellitus is characterized In two broad groups based on clinical 
manifestations, namely the non-insulin-dependent diabetes (NiDDM) or maturity onset fomn. also 

15 known as Type 2 and the insulin-dependent diabetes (IDDM) or juvenile onset form, also known as 
Type 1. In the general population, diabetes mellitus occurs with a prevalence of approximately 1%, 
with one fourth of these being the Typel. In its nnost fully expressed dinlcal form, diabetes melGtus 
manifests itself as a series of hormone-induced metabolic abnormalities which eventually lead to 
serious, long4erm and deblfitating complications such as glucosuria, ketonuria, growth arrest, and 

20 negative nitrogen balance. These conditions can ultimately lead to death. FamiGal 

hyperproinsuGnenia is a genetic disorder characterized by a mariced increase in serum prolnsulin-like 
molecules. The cause of this disease is an amino add sut)stitution which results in incomplete 
cleavage of proinsulln by the proteases which fonn insuQn. 

Type 1 diabetes arises for examfrfe. when patients lack beta-cells producing Insufin in their pancreatic 
25 glands or when the produced insulin is inactive due to mutation(s). Type 2 diabetes occurs in patients 
with an impaired beta cell function. Type 1 diabetic patients are currently treated with insufin, while 
the majority of type 2 diabetic patients are treated either with sulfonylureas that stimulate beta cell 
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function or with agents that enhance the tissue sensitivity of the patients towards insuOn (e.g., 
metfomnin) or with insulin. 

Today, insulin administration to diabetic patients Is the primary therapeutic means for controlling the 
disease. In the treatment of diat>etes meliitus, many varieties of insulin preparations have been 
suggested and used. Some of the preparations are fast acting and other preparations have more or 
less prolonged actions. Such a prolonged action may be obtained by administering the insulin as a 
suspension of insulin crystal wtuch can be otitained by crystalHzation of insulin in the presence of zinc 
(such as LENTE; Novo Terapeutisk Laboratorium) or by ciystallizafion of insulin in the presence of 
zinc and protantine (such as NPH-insuHn). 

The human insulin monomer, a 6000 datton protein, is composed of two chains, the 21 amino acid A- 
chain and the 30 amino acid B-chain. Insulin is synthesized in pancreatic beta cells located within the 
isl^ of Langerhans as a precursor form that Is post-translatlonatiy processed to the mature two 
polypeptide chain active hormone. In the biologically active human Insulin, the A and B chains are 
linked with one another via two cysteine bridges, and a further cysteine bridge occurs within the A- 
chain. The following cysteine residues are linked with one another in human insulin: A6-A-1 1. A7-B7, 
and A20-B-19 (the letters A and B stand for the amino acid chain, and the numbers for the posjtion of 
the cysteine residues counted from the amino to carfooxyl end of each chain; see Figure 1). 

Based on the functional analysis of insulins from various species and several insulin analogs or 
mutants, some pertinent properties of insulin have emeiged with respect to the anrnno acid sequence 
(see also Figure 1): biologically active insulin has three disulfide bonds; Bl-Phe is present in all known 
mammalian insulins; Al-Gly; a terminal tripeptide sequence (A19-21) Tyr-Cys-Asn (removal of A21- 
Asn by carboxypeptidase digestfon results in a loss of activity of >90%); an invariant sequence at B24- 
26. Phe-Phe-Tyn B12-Val is highly conserved; A2-A3, lle-Val is highly conserved; B5-His and B22-Ang 
are invariant in insulins of high potency; Invariant surface residues (highly conserved) Al-Gly, A4-Glu, 
A5-Glh. A7-Cys, A19-Tyr, A21-Asn. and B7-Cys. 

To this end, variants of insuGn sequences, appKcations. production procedures and assays are known; 
see for example U.S. Patent Nos. 4,421 ,685 (reports process for producing insulin); 4,992,417 
(reports superactive insuGn analogues); 5.008,241 (reports insuOn analogs characterized by amino 
add residue sutrstitutions of N21 in the A-chain, with a resulting improvement in stability of insulin 
solutions at addle pH levels); 5,506.202 (reports preparation and use of insulin derivatives comprising 
various anruno add substitutions); 5.514,646 (reports analogs of human InsuDn nu>dified at position 29 
of the B-chain that have nnodified physico-cheriiicat and pharmacokinetic properties and are useful in 
the treatment of hyperglycentia); 5,559,094(reports analogs of human insuOh containing an aspartic 
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add at position 1 of the B-chain and 5.618,913 (reports rapid acting human insulin analogues having 
anino add residue substitution and less tendency to self-assodate into dimers, tetramers, hexaniers 
or polymers); 5,621,073 (reports a process for purification of insulin or an insulin acetylated at position 
A9); 5.663,291 (reports a process fr obtaining insulin having con-ectly linked cysteine bridges from a 
corresponding proinsulin amino add chain): 5,700,662 (reports process for preparing insulin analogs 
comprising a modification of position 29 of the B-chaln);6,034.054 (reports a monomeric insulin analog 
fonnulation stabifized against aggregation): all of which are expressly, incorporated by reference and 
further by Marid et al. JHoppe Seylers Z, Phsiol. Chem. 360(1 1):1 61 9^2 (1979)1; Hu et al. 
{Biochenistry 32(10)2631-5 (1993)1; Schwartz et al. proc. Natl. Acad. Sd. U.S.A. 84(1 8):6408.1 1 
(1987)1; Kitagawa et al. piochemistry 23(7):1405-13 (1984)1; Kobayashi et al. [Biochem. Biophys. Res. 
Commun. 107(1):329-36 (1982)1; Shoelson et al. IBiochemistiy 31(6):1757-67 (1992)1; and references 
dted therein, all of which are expressly incorporated as reference. 

Human insulin in solution is known to exist In many molecular fonns, such as the monomer, the dimer, 
the tetramer and the hexamer [Blundell et al.. in Advances in Protein Chemistry. Academic Press, New 
Yoric and London, Vol. 26. pp279-330 (1972)]. with the oligomer forms bdng favored at high insulin 
concentrations and the monomer being the active fpmi of insulin. Insulin in the bloodstream Is highly 
dilute, being 10'^' to 10^ M and is primarily in monomer fbmn. The much more concentrated Insulin 
stored in the beta cell of the pancreas and in the usual administeraWe solution is largely in the non- 
acUve hexamer forth, as the well-known 2 zinc hexamer (see betow). The delayed absorption 
phenomena [Binder. Diabetes Care 7(2):188-99 (1984)1 is in some large part attributable to the time 
required for the insulin to disassociate from hexamer. tetramer and dimerfomi into the acUve 
monomer form. 

In the presence of zinc (Zn). natural human insulin associates to a 2 Zn-hexamer that functions as an 
altosteric protein. Phenolic ligands or certain salts are capable of indudng a conformational transitton. 
resulting in the N-temtinal 8 amino ackJs of the B-chain converting from an extended confbnnation to 
an Q-helix. This confomiational state Induced by phenolic ligjands has been refen^ed to as the R state 
and the apoinsulin form as the T state. The R slate is more compact, less flexible, and the Zn 
exchange is retanled compared to the T state [Derewenda et aL. Nature 338(6216):594.596 (1989)1. 
A stable intennediate state. T3R3 has been identified that has one trimer in the T state and the other in 
the R state [Chothia et al. Nature 302(5908):500-505 (1983)1, The T3R3 state was fonrelly known as 
the 4-Zn insulin structure wtiich is induced by salts (e.g. chloride) or by limited amounts of phenolics 
IKmger. et al., Btol. Chem. Hoppe-Seyler 371(8):669-673 (1990)]. 

The different altosteric states of insulin hexamer have been best characterized in the crystal state by 
X-ray crystaltography pentley et al.. Nature 261(5556):166.168 (1976); Smith & Dodson. Btopolymers 
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32{4):441-445 (1992)]. in solution by proton NMR. circular dichroism [Renscheidt et al.. Eur. J. 
Biochem. 142(1):7-14 (1984)]. and visible absorption spectroscopy of Co2^ substituted Insulins [Brader 
et al., Btochenustry 30(27):6636-6645 (1991)]. The bioioglcal significance of insulin allosterism has not 
been fully elucidated. The biologically active form of insulin is thought to be a nnonomer due to the 
dilute concentrations of insulin in the blood circulation [Frank et al.. Diabetes 21(2):Suppl. 2:486-491 
(1972)]. A receptor-mediated conformational change in the insulin confomiation is thought to be 
required for binding [see Hua et al. Nature 354(6350):238-41 (1991); Bao et al. Proc. Natl. Acad. Sci. 
\J.SA. 94(7):2975-80 (1997)]. For the medicinal use of insulin, the T->R confomiations have 
important consequences. Most fonmulations of insulin are solutions or suspensions that contain 
phenolics that function as preseiyatives against bacterial contamination. The phenolic concentrations 
in insulin fonmulations are 2-10 times that necessary to induce the R conformation (Kruger et al., 
supra). Hie presence of phenolics in insulin fommjlations has also important consequences on the 
shelf-life stability [Brange et al.. Phann. Res. 9(6):715-726 (1992); Brange et al., Phami. Res. 
9(6):727-734 (1992); Brange & Langkjaer, Acta Pharm. Nord..4{3): 149-1 58 (1992)] and possibly the 
time action profile. Minimizing degradation of insulin fomriulations is extremely important in reducing 
undesirable side effects of insulin therapy. 

The crystal structures of vaiiotis recombinant insulin molecules are solved. The structures can be 
obtained frohi the Research CollaboFatory for Stmctural Bioinfomiatics as entries into the Protein Data 
Bank (PDB). Insulin entries include wild type insulin (e.g.. see PDB entries 1ZEH, 1ZNJ, 1ZN1. 1XDA. 
41NS, and 9INS). insulin analogs or mutants (e.g.. PDB entries 1IOH, 1I0G, 1B9E, 1BZY. 12EI, 1A7F, 
1HUI, 1LPH. and IIZA). R6 (wild type or analog) insulin hexamers (e.g., PDB entries 5AIY, 4AIY, 
3AIY, 2AIY, 1AIY, 1AI0. 1QIZ, 1QJ0, and 1QiY); insulin compl«ced with phenol (e.g., PDB entry 
1ZEG), Insufin complexed with 4-hydroxyben2amide (e.g.. PDB entry 1BEN). insulin (wild type or 
analog) complexed with Zn ions (e.g., PDB entries 17YM, 1TYL. 7INS, 1TRZ, and 1IZB). and other 
variants (e.g., PDB entries IMP J. 2TCI. 3MTH, 6INS, and 2INS) and insulins in various pH solutions 
(e.g.. PDB entries 1DPH. 1CPH, 1BPH, and 1APH), all of which are expressly incorporated by 
reference. 

When carrying out protein engineering to nuxiify protein properties, usually one had to select from the 
following options: (i) site-specific mutagenesis and Ql) random nruitagenesis of the nucleic add 
encoding the protein, or (iii*) post-translational cherracal modifications. No matter which method of 
protein engineering Is used, a key aspect is determining which amino acids to modify, because few 
choices will improve the properties of the protein. The available crystal structure of insufin allows a 
completely different approach by using computational protein design and the generation of more 
stable proteins or protein variants with an altered activity. Several groups have applied and 
experimentally tested systematic, quantitative methods to protein design with the goal of devek>ping 
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general design algorithms (Hellinga et at.. J. Mol. Biol. 222: 763-785 (1991); Huriey et at.. J. Mol. Biol. 
224:1 143-1 154 (1992): DesjarlaisI et al.. Protein Science 4:2006-2018 (1995): Harfoury et at.. Proc. 
Natl. Acad. Sd. U.SA 92:8408-8412 (1995); Klemba et at., Nat. Stmc. Biol. 2:368-373 (1995); 
Nautlyal et a!,. Biochemistry 34:1 1645-1 1651 (1995); Betzo et at.. Biochemistry 35:6955-6962 (1996); 
Dahiyat et al.. Protein Science 5:895-903 (1996); Dahiyat et at., Science 278:82-87 (1997); Dahiyat et 
at., J. Mol. Biol. 273:789-96; Dahiyat et at.. Protein Sci. 6:1333-1337 (1997); Jones, Protein Science 
3:567-574 (1994); Konoi, et a!.. Proteins: Stmcture, Function and Genetics 19:244-255 (1994)). Tliese 
algorithms consider the spatial positioning and steric complementarity of side chains by explicitly 
modeling the atoms of sequences under consideration. In particular . WO98/47089, and U.S.S.N. 
09/127,926 describe a system for protein design; both are expressly inoofporated by reference. 

A need still exists for proteins exhibiting both significant stability and insulin activity. Accordingly, it is 
an object of the invention to provide insulin activity (lA) proteins, nudeic adds and antibodies for the 
treatment of insulin-related disorders. 

SUMMARY OF THE INVENTION 

In accordance with the objects outlined above, the present invention provides non-naturalty occurring 
insulin activity (lA) proteins (e.g. the proteins are not found in nature) comprising amino add 
sequences that are less than about 98% Identical to human Insulin. The lA proteins have at least one 
altered biological property of an insulin protein; for example, the I A proteins will be more stable than 
insulin and bind to a cell comprising an insulin receptor. Thus, the invention provides lA proteins with 
amino add sequences that have at least about 1-20 annsno acid substitutions as compared to the 
human insulin sequence shown in Figure IB. 

In a further aspect, the present Invention provides non-naturalty occurring lA protein conformers that 
have three dimensional backbone stmctures that substantially correspond to the three dimensional 
l>ackbone structure of insufin. The amino add sequence of the lA protein conformer and the an^no 
add sequence of insulin are less than about 98% identical. 

In an additional aspect, the changes are selected from the amino add residues at positions selected 
frorri positions A2, AS, AS, A7, All. A15. A16. A19. A20, B2. B7, 811. B15. B18, B19, 822. and 
B24. In a preferred aspect of this embodiment, the changes are sutjstitutions selected from the group 
of A7-S, A7-E, B2-E, B2-T, B4-Y, B7-Y. B4-F, B7-E. and B7-D. 
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In a preferred aspect, the changes are selected from the amino add residues at positions selected 
from positions B5 and B14. In a preferred aspect of this embodiment, the changes are substitutions 
selected from the group ofSS-F, B5-W, B14-F, B14-W, B14-Y, and B14-I. 

In an additional aspect, the changes are selected from the amino add residues at positions selected 
from positions A1. A10. A16. A17, A19, 81. B2. 84, 88, 811, 812, 814. 825, 826, 827, and 828. In a 
prefen-ed aspect of this embodiment, the changes are substitutions selected from the group of A1-N, 
A10-Q. A16-Y, A17.Y, A19-F, B1-D. B2-K. B4-F, B8^ BIH, B12-R, B14-W, B25-N, 826-F. B27-D, 
and B28-N. 

In a further aspect, the invention provides recombinant nucleic adds encoding the non-naturally 
occurring lA proteins, expression vectors comprising the recombinant nucleic adds, and host cells 
comprising the recombinant nucleic adds and expression vectors. 

in an additional aspect, the invention provides methods of prbdudng the lA proteins of the invention 
comprising culturing host cells comprising the recombinant nucleic adds under conditiohs suitable for 
expression of the nudeic acids. The proteins may optionally be recovered. In a further aspect, the 
invention provides pharmaceutical connpositions comprising an lA protein of the invention and a 
pharmaceutical canler. 

In an additional aspect, the invention provides methods for treating an insulin responsive condition 
comprising admnistering an lA protein of the invention to a patient The insuGn responsive condition 
indudes a disorder of cartK)hydrate metabolism, type 1 diabetes, and type 2 diabetes. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1A depicts the amino add sequence of the insulin precursor (GenBank accession #P01308, 
#AAA59173). Amino add residues 1-24 represent the signal peptide; amino add residues 25-54 
represent the mature B-chain; an^no add residues 90-1 1 0 represent the nnature A-cha!n. 

Figure IB depicts the amino add sequence of human insulin [GenBank accession ^9122; Nicol and • 
Smith, Nature 187:483-485 (I960)], which is used herein for PDA design and for reference of amino 
add positions. The A-chain comprises amino add residues 1-21 (sometimes refen'ed to as A1 
through A 21) and the B-chaIn comprises residues 22-51 (sometinnes referred to as B1 through 830). 
respectively. 
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* Figure 1 C depicts the amino add sequences of the A-chains of human insulin (ITRZrA and ITRZiC) 
and B-chains of human Insulin (1TRZ:B and 1TRZ:D) as used in detenmination of the structure of 
insufin (T3R3) complex with two zinc ions (PDB entry ITRZ; Ciszak and Smith, Biochen^stry 
33(6):1512-7 (1994)] and secondary structure elements. Secondary structure element legend: H, 
5 alpha helix (44ie{ix); 6. residue in Isolated beta bridge; E, extended strand, participates In beta ladder. 
G. 310 helix (3-henx); I. pi heOx (S-hefix); T. hydrogen bonded tum; S. bend. 

Rgure 2A depicts the structure of a wild type insulin monomer with side chains for cBsulfide bonds (A6- 
All, A7-B7. and A20-B19) and the B-Ala14 side chain shown. 

Figure 2B depicts the structure of a wild type insulin hexamer. 

1 0 Figure 2C depicts a doseup of the B1 4. B5 design region in the insulin hexamer. 

Figure 3 shows preferred lA protein sequences from PDA designs involving disulfide replacement, 
Amino add changes when compared to wild type human insulin are indicated in bold and are 
underiined. 

Rgure 3A shows a preferred I A protein sequence from PDA design 'cysl', 

15 Rgure 38 shows a preferred lA protein sequence from PDA design 'cys77a'. 

Figure 3C shows a preferred lA protein sequence from PDA design •cys77b'. 

Rgure 3D shows a preferred lA protein sequence from PDA design 'cys77d'. 

Rgure 3E shows a prefeaed lA protein sequence from PDA design *cys77dV. 

Figure 3F shows a preferred lA protein sequence from PDA design 'hefix 24\ 

2 0 Figure 3G shows a prefen-ed lA protein sequence from PDA design *cys-4\ For this design, a In the 
sequence indicates a deletion. 

Figures 4A-4G show preferred lA protein sequences from PDA designs involving mutations which 
promote insulin hexamer fomiation. Aituno add changes when compared to wild type human insufin 
are Indicated In bold and are underflned. 
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Figure 5 shows preferred lA protein sequences from PDA designs involving global redesigns for 
improved stability. Amino add changes when compared to wild type human insulin are indicated in 
bold and are underilned. 

Rgure 5A shows a preferred lA protein sequence from PDA design *trz_06'. 

Rgure SB shows a preferred lA protein sequence from PDA design 1rz_7b*. 

Rgure 5C shows a preferred lA protein sequence from PDA design *trz_b8'. 

Rgure 6 depicts the synthesis of a fulMength gene and ail possit>le mutations by PGR. Overlapping 
ofigonudeotides corresponding to the fulMength gene (black bar. Step 1) and comprising one or more 
desired mutations are synthesized, heated and annealed. Addition of DMA polymerase to the 
annealed ofigonudeotides results in the 5* to 3* synthesis of DMA (Step 2) to produce longer DMA 
fragments (Step 3). Repeated cycles of heating, annealing, and DNA synthesis (Step 4) result in the 
production of longer DNA, including some fulWength molecules. These can be selected by a second 
round of PGR using prinrwrs indicated by an^ows) con-esponding to the end of the fulMength gene 
(Step 5). 

Rgure 7 depicts a preferred scheme for synthesizing an lA protein library of the invention. The wild 
type gene, or any starting gene, such as the gene for the global minima gene, can be used. 
Oligonudeotides comprising sequences that encode different amino adds at the different variant 
positions Ondicated in the Rgure by box 1 , box 2, and box 3) can be used during PGR. Those primers 
can be used in combination with standard primers. This generally requires fewer oligonudeotides and 
can result In fewer errors. 

Rgures 8A and 8B depict an overtapping extension method. At the top of Rgure 8A is the template 
DNA showing the locations of the regions to be mutated (black boxes) and the binding sites of the 
relevant primers (aaows). The primers R1 and R2 represent a pool of primers, each containing a 
different mutation; as described herein, this may be done using different ratios of primers if desired. 
The variant position Is flanked by regions of homology suffident to get hybridization. Thus, as shown in 
this example, ofigos R1 and F2 comprise a region of honrx)logy and so do oligos R2 and F3. In this 
example, three separate PGR reactions are done for step 1 . The first reaction contains the template 
plus oGgos Fl and R1 . The second reaction contains template plus otigos F2 and R2, and the third 
contains the temfdate and ofigos F3 and R3. The reaction products are shown. In Step 2, the products 
from Step 1 tube 1 and Step 1 tube 2 are taken. After purification away from the primers, these are 
added to a fresh PGR reaction together with F1 and R4. During the denaturation phase of the PGR, 
the overtapping regions anneal and the second strand is synthesized. The product is then amplified by 
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the outside primers, F1 and R4. In Step 3, the purified product from Step 2 is used in a third PGR 
reaction, together with the product of Step 1 , tube 3 and the primers F1 and R3. The final product 
corresponds to the full length gene and contains the required mutations. Alternatively, Step 2 and 
Step 3 can be perfonned in one PGR reaction. 

Figures 9A and 9B depict a ligation of PGR reaction products to synthesize the libraries of the 
invention. In this technique, the primers also contain ah endonudease restriction site (RE), either 
generating blunt ends, 5' overtianging ends or 3* overhanging ends. We set up three separate PGR 
reactions for Step 1. The first reaction contains the template plus ollgos F1 and R1 . The second 
reaction contains template plus oligos F2 and R2, and the third contains the template and ollgos F3 
and R3. The reaction products are shown. In Step 2. the products of Step 1 are purified and then 
digested with the appropriate restriction endonudease. The digestion products from Step 2, tube 1 and 
Step 2, tube 2 are Rjgated togther with DIMA Ogase (Step 3). The products are then amplified in Step 4 
using oligos F1 and R4. The whole process is then repeated by digesting the amplified products, 
ligating them to the digested products of Step 2. tube 3, and then amplifying the final product using 
oligos F1 and R3. It would also be possible to ligate all three PGR products from Step 1 together in 
one reaction, providing the two restriction sites (REIand RE2) were different 

Rgure 10 depicts Wunt end figation of PGR products. In this technique, ofigos such as F2 and R1 or 
R2 and F3 do not overiap. but they abut. Again three separate PGR reactions are performed. The 
products from tube 1 and tube 2 (see Figure 9A, Step 1) are ligated. and then ampfified with outside 
primers F1 and R4. This product is then ligated with the product from Step 1 . tube 3. The final 
products are then amplified witii primers F1 and R3. 

DETAILEO DESGRIPTION OF THE INVENTION 

The present invention is directed to novel proteins and nucleic adds possessing insulin activity 
(sometimes referred to herein as lA proteins" and "lA nudeic adds"). The proteins are generated 
using a system previously described in WO98/47089 and U.S.S.Nos. 09/058.459, 09/127.926. 
60^104.612. 60/158,700, 09/419,351. 60/181.630. 60/186.904, and U.S patent application, entitied 
Protein Design Automation For Protein Ubraries (RIed: /\pril 14. 2000; Inventor: Bassil Dahiyat), all of 
which are expressly incorporated by reference in their entirety, ttiat is a computational modeling 
system ttiat allows the generation of extremely stable proteins without necessarily disturbing the 
biological functions of the protein itself. In this way. novel lA proteins and nudeic adds are generated, 
that can have a plurality of mutations In connparison to the wild-type enzyme yet retain significant 
activity. 
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Generally, there are a variety of computational methods that can be used to generate the I A proteins 
of the invention. In a preferred embodiment, sequence based methods are used. Alternatively, 
structure based methods, such as PDA, described in detail below, are used. 

Similarly, molecular dynamics calculations can be used to computationally screen sequences by 
Individually calculating mutant sequence scores and compiOng a rank ordered list 

In a preferred embodiment, residue pair potentials can be used to score sequences (Miyazawa et al.. 
Macromolecules 18(3):534-552 (1985). expressly Incorporated by reference) during computational 
screening. 

In a preferred embodiment, sequence profile scores (Bowie et aL. Science 253{5016):164-70 (1991), 
incorporated by reference) and/or potentials of mean force (Hendlich et al.. J. Mol. Biol. 216(1):167- 
180 (1990). also Incorporated by reference) can also be calculated to score sequences. These 
methods assess the match between a sequence and a 3D protein structure and hence can act to 
saeen for fidelity to the protein structure. By using different scoring functions to rank sequences, 
different regions of sequence space can be sampled in the computational screen. 

Furthermore, scoring functions can be used to screen for sequences that would create metal or co- 
factor binding sites in the protein (HelUnga, Fold Des. 3(1):R1-8 (1998), hereby expressly incorporated 
by reference). Similarty. scoring functions can be used to screen for sequences that would create 
disulfide bonds in the protein. These potentials attempt to specifically modify a protein structure to 
Introduce a new structural motif. 

m a preferred embodiment, sequence and/or structural alignment programs can be used to generate 
the lA proteins of the invention. As is known in the art. there are a number of sequence-based 
angnment programs; including for example, Smith-Waterman searches. Needleman-Wunsch. Double 
Affine Srrith-Waterman. frame search. Gribskov/GCG profile search. Gribskov/GCG profile scan, 
profile frame search. Bucher generalized profiles. Hidden Maricov models. Hframe. Double Frame, 
Blast, Psi-Blast, Clustal, and Gene\A/ise. 

As is known in tiie art. there are a number of sequence aKgnment methodologies that can be used. 
For example, sequence homology based afignment methods can be used to create sequence 
alignmenls of proteins related to the target structure (Wlschul et al.. J. Mol. Biol. 215(3):403-410 
(1990). Altschul et al.. Nucleic Adds Res. 25:3389-3402 (1997). both incorporated by reference). 
These sequence alignments are then examined to determine ttie observed sequence variations. 
These sequence variations are tabulated to define a set of lA proteins. 
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Sequence based alignmenls can be used in a variety of ways. For example, a number of related 
proteins can be aligned, as is known In the art, and the Variable" and 'consented" residues defined; 
that is, the residues that vaiy or remain identical between the family members can be defined. These 
results can be used to generate a probability table, as outlined below. Similarly, these sequence 
variations can be tabulated and a secondary Obrary defined from them as defined below. Alternatively, 
the allowed sequence variations can be used to define the anino adds considered at each position 
during the computational screening. Another variation is to bias the score for antino adds that occur In 
the sequence alignment, thereby increasing the Gkelihood that they are found during computational 
screening but still allowing consideration of other amino adds. This bias would result in a focused 
library of lA proteins but would not eliminate from consideration amino adds not found in the 
afignment. In addition, a number of other types of bias may be introduced. For example, diversity may 
be forced; that is, a "consewed* residue Is chosen and altered to force diversity on the protein and 
thus sample a greater portion of the sequence space. Altematively, the positions of high variability 
between fantily members (i.e. low conservation) can be randomized, either using all or a subset of 
arnno adds. Similarty, outlier residues, either positional outiiers or side chain outliers, may be 
eliminated. 

Sin^larly, structural aOgnment of structurally related proteins can be done to generate sequence 
aOgnments (Orengo el al.. Stnicture 5{8):1093-108 (1997); Holm et al. Nucleic Adds Res. 26{1):316-9 
(1998). both of which are incorporated by reference). These sequence alignments can then be 
exan^ned to deterrnne the observed sequence variations. Ubraries can be generated by predicting 
secondary structure from sequence, and then selecting sequences that are compatible with the 
predicted secondary stnicture. There are a nun*er of secondary stnidure prediction methods such 
as heOx-coil transition theory (Munoz and Serrano. Biopolymers 41:495, 1997). neural networks, local 
structure aOgnment and ottiers (e.g.. see in SelWg et al.. Bidnfomnatics 15:1039^6. 1999). 

Sin^'larty. as outiined above, other computational methods are known, including, but not Omited to, 
sequence profifing [Bowie and Bsenberg. Sdence 253(5016):164.70. (1991)1, rotamer Obrary 
selections [Dahiyat and Mayo, Protein Sd. 5(5):895-903 (1996); Dahiyat and Mayo, Sdence 
278(5335):82-7 (1997); Desjartais and Handel, Protein Sdence 4:2006-2018 (1995); Harbury et al, 
Proc. Natl. Acad. Sd, U,SA 92(18):8408-8412 (1995); Kono et al.. Proteins: Structure, Function and 
Genetics 19:244-255 (1994); Helfinga and Richards, Proc. Natl. Acad. Sd. U.SA 91:5803-5807 
(1994)1; and residue pair potentials [Jones, Protein Sdence 3: 567-574. (1994)1; PROSA IHeindfich el 
ai;. J, Mol. BioK 216:167-180 (1990)]; THREADER [Jones et al.. Nature 358:86-89 (1992)]. and other 
inverse folding methods such as those described by Simons et al. [Proteins, 34:535^, (1999)], 
Levitt and Gerstdn [Proc. Nati. Acad. Sd. U.SA, 95:5913^5920, (1998)]. Godzik and Skolnidc [Proc. 
Nati. Acad. Sd. U.SA. 89:12098-102. (1992)]. Godzik et al. (J. Mol. Biol. 227:227-38. (1992)] and two 
profile methods (Gribskov et al. Proc. Nati. Acad. Sd. U.SA 84:4355-4358 (1987) and Fischer and 
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Bsenberg. Protein Sd. 5:947-955 (1996). Rice and Eisenberg J. Mol. Biol. 267:1026-1038(1997)]. all 
of which are expressly incorporated by reference. In addition, other coniputational methods such as 
those described by Koehl and Levitt (J. Moh BloL 293:1161-1181 (1999); J. Mol. Biol. 293:1183-1193 
(1999); expressly incorporated by reference) can be used to create a protein sequence library which 
can optionally then be used to generate a smaller secondary library for use in experimental screening 
for improved properties and function. In addition, there are computational methods based on 
forcefield calculations such as SCMF that can be used as well for SCMF, see Delarue et al. Pac. 
Symp. Biocomput. 109-21 (1997); Koehl et al.. J. Mol. Biol. 239*249-75 (1994); Koehl et a!.. Nat 
Stmct, Biol. 2:163-70 (1995); Koehl et al.. Cun*. Opin. Struct. Biol. 6:222-6 (1996); Koehl et al., J. Mol. 
Biol. 293:1183-93 (1999); Koehl et al.. J. Mol. Biol. 293:1161-81 (1999); Lee J.. Mol. Biol, 236:918-39 
(1994); and Vasquez Biopolymeis 36:53-70 (1995); all of which are expressly incorporated by 
reference. Other forcefield calculations that can be used to optimize the confonnation of a sequence 
within a computational method, or to genenite de novo optimized sequences as outlined herein 
include, but are not Umrted to, OPLS-AA [Jorgensen et al., J. Am. Chem. Soc. 1 1 8:1 1 225-1 1 236 
(1996); Jorgensen,. W.L.; BOSS, Version 4.1; Yale University: New Haven. CT (1999)]; OPLS 
[Jorgensen et al., J. Am. Chem. Soc.1 10:1657ff (1988); Jorgensen et al.. J Am. Chem. Sod 12:4768ff 
(1990)1; UNRES (United Residue Forcefield; Liwo et al., Protein Science 2:1697-1714 (1993); Liwo el 
al.. Protein Science 2:1715-1731 (1993); Liwo et al.. J. Comp. Chem. 18:849-873 (1997); Uwo et al., J. 
Comp. Chem. 18:874-€f84 (1997); Liwo et al.. J. Comp. Chem. 19:259-276 (1998); Forcefield for 
Protein Structure Prediction (Uwo et al.. Proc. Natl. Acad. Sd. U.S.A . 96:5482-5485 (1999)1; 
ECEPP/3 [Uwo et al., J Protein Chem. 13(4):375-80 (1994)]; AMBER 1.1 force field (Weiner et al.. J. 
Am. Chem. Soc. 106:765-784); AMBER 3.0 force field lU.C. Singh et al.. Proc. Natl. Acad. Sd. U.S.A,. 
82:755-759 (1985)1; CHARMM and CHARMM22 (Brooks et al., J. Comp. Chem. 4:187-217); cvff3,0 
(Dauber-Osguthorpe et al., Proteins: Structure. Function and Genetics, 4:31-47 (1988)1; (Maple 
et al.. J. Comp. Chem. 15:162-182); also, the DISCOVER (cvff and cfl91) and AMBER forcefields are 
used in the INSIGHT molecular modeling package (Biosym/MSI. San Diego Catifomia) and HARMM is 
used in ttie QUANTA molecular modeDng package (Biosym/MSI. San Diego California), ail of which 
are expressly incorporated by reference. In fact, as is outlined below, these forcefieW methods may 
be used to generate the secondary fibrary directly; that is, no primaiy fibraiy is generated; rather, these 
methods can be used to generate a probability table from which the secondary library is directly 
generated. 

In a preferred embodiment, the computational method used to generate the primary library is Protein 
Design Automation (PDA), as Is described in U.S.S.N.S 60/061 .097. 60/043,464, 60/054,678. 
09/127,926, 60/104.612, 60/158,700. 09/419,351, 60/181630, 60/186,904, U.S patent appOcation, 
entitied Pmtein Design Automation ForPmte'tn Utmries (Filed: April 14, 2000; Inventor. Bassil 
Dahiyat) and PCT US98/07254, all of wrtiich are expressly incorporated herein by reference. Briefly. 
PDA can be described as follows. A known protein structure is used as the starting point. The 
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residues to be optinruzed are then identified, which may be the entire sequence or subset(s) thereof. 
The side chains of any positions to be varied are then removed. The resulting structure consisting of 
the protein backbone and the remaining side chains is called the template. Each variable residue 
position is then preferably classified as a core residue, a surface residue, or a boundary residue; each 
classification defines a subset of possible amino add residues for the position (for example, core 
residues generally will be selected from the set of hydrophobic residues, surface residues generally 
will be selected from the hydrophilic residues, and boundary residues may be either). Each arnno 
add can be represented by a discrete set of all allowed confbrmers of each side chain, called 
rotamers. Thus, to arrive at an optimal sequence for a backbone, all possible sequences of rotamers 
must be screened, where each backbone position can be occupied either by each amno add in all its 
possible rotameric states, or a subset of anruno adds, and thus a subset of rotamers. 

Two sets of interactions are then calculated for each rotamer at every position: the interaction of the 
rotamer side chain with all or part of the backbone (the 'singles' energy, also called the 
rotamerAemplate or rotamer/backbone energy), and the interaction of the rotamer side chain with all 
other possible rotamers at every other position or a subset of the other positions (the "doubles" energy, 
also called the rotamer/rotamer energy). The energy of each of these interactions is calculated 
through the use of a variety of scoring functions, which include the energy of van der Waal's forces, 
the energy of hydrogen bonding, the energy of secondary structure propensity, the energy of surface 
area solvation and the electrostatics. Thus, the total energy of each rotamer interaction, botfi with the 
backbone and other rotamers. Is calculated, and stored in a matrix form. 

The discrete nature of rotamer sets allows a simple calculation of the number of rotamer sequences to 
be tested. A backbone of length n with m possible rotamers per position will have irf possible rotamer 
sequences, a number which grows exponentially with sequence length and renders the calculations 
eittier unwieldy or Impossible In real time. Accordingly, to solve tiiis combinatorial search problem, a 
•Dead End Efiriiination" (DEE) calculation is performed. The DEE calculation is based on ttie fact that 
if the worst total Interaction of a first rotamer is still better than tiie best total interaction of a second 
rotamer. then the second rotamer cannot be part of ttie global optimum solution. Since the energies of 
all rotamers have already been calculated, the DEE approach only requires sunrs over the sequence 
lengtii to test and efinranate rotamers, which speeds up the calculations considerably. DEE can be 
rerun comparing pairs of rotamers, or combinations of rotamers, which will eventually result in the 
detennination of a single sequence wNch represents the global optimum energy. 

Once the global solution has been found, a Monte Carlo search may be done to generate a rank- 
ordered list of sequences in the neighborhood of the DEE solution. Sterting at ttie DEE solution, 
random positions are changed to otiier rotamers, and the new sequence energy is calculated. If the 
new sequence meets the criteria for acceptanpe, it is used as a starting point for another jump. After a 
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predetermined number of jumps, a rank-ordered fist of sequences is generated. Monte Carlo 
searcNng is a sampling technique to explore sequence space around the global minimum or to find 
new local minima distant in sequence space. As is more additionally outlined below, there are other 
sampling techniques that can be used, including Boltzman sampling, genetic algorithm techniques and 
simulated annealing. In addition, for all the sampling techniques, the kinds of jumps allowed can be 
altered (e.g. random junnps to random residues, biased jumps (to or away from wild-type, for 
exanple). jumps to tiased residues Qo or away from dmilar residues, for example), etc.). Similarty. 
for all the sampling techniques, the acceptance cnteria of whether a sampfing jump is accepted can be 
altered. 

As outlined in U.S.S.N. 09/127.926. the protein backbone (comprising (for a naturally occumng 
protein) the nitrogen, the cartsonyl carbon, the a-carbon, and the carbonyl oxygen, along with the 
direction of the vector from the a-carbon to the p-carbon) may be altered prior to the computational 
analysis, by varying a set of parameters called supersecondary structure parameters. 

Once a protein structure backbone is generated (with alterations, as outlined above) and input into the 
computer, expBdt hydrogens are added if not induded within the structure (for exanfiple, if the structure 
was generated by X-ray crystallography, hydrogens must be added). After hydrogen addition, energy 
ntinlmization of the structure is run, to relax the hydrogens as well as tiie other atoms, bond angles 
and bond lengths. In a preferred embodiment, this Is done by doing a number of steps of conjugate 
gradient minimization (Mayo et al., J. Phys. Chem. 94:8897 (1990)] of atomic coordinate positions to 
minimize the Dreiding force field with no electrostatics. Generally from about 10 to about 250 steps is 
preferred, witii about 50 being most prefen^ed. 

The protein backbone structure contains at least one variable residue position. As is known in the art. 
the residues, or amino adds, of proteins are generally sequentially numbered starting with the N- 
teminus of the protein. Thus a protein having a methionine at it's N-teminus is said to have a 
methionine at residue or an^no add position 1 . witii the next residues as 2. 3. 4. etc. At each position, 
the wild type (i.e. naturally occurring) protein nnay have one of at least 20 amino adds, in any number 
of rotamers. By "variable residue portion* herein Is meant an anrino add position of the protein to be 
designed that Is not fixed in the design method as a spedfic resWue or rotamer. generally the wild-type 
residue or rotamer. 

In a preferred embocfiment. all of the residue positions of the proton are variable. That is. every anrano 
add side chain may be altered In the methods of the present Invention. This is particularty desirable 
for smaller proteins, although the present methods allow the design of larger proteins as well. While 
there is no ttieoretical fimrt to the length of the protein which may be designed tWs way. there is a 
practical computational fimit. 
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In an alternate prefeired emlwdiment, only some of ttie residue positions of tlie protein are variable, 
and the remainder are •fixed', that is, they are identified in the three dimensional structure as being in 
a set confonnatlon. In some embodiments, a fixed position is left in its original conformation (which 
may or may not con-elate to a specific rotamer of the rotamer Bbraiy being used). Altemafively. 
residues may be fixed as a norHwId type residue; for example, when known site-directed mutageneas 
techniques have shown that a particular residue is desirable (for example, to efintinate a proteolytic 
site or alter the substrate spedfidty of an enzyme), the residue may be fixed as a particular amino 
add. Altematwely. the methods of the present invention may be used to evaluate mutations de novo, 
as Is discussed below. In an altemate preferred embodiment, a fixed position may be floated"; the 
amno add at that position is fixed, but different rotamers of that amino add are tested. In this 
embodiment, the variable residues nray be at least one. or anywhere from 0.1 % to 99.9% of the total 
number of residues. Thus, for example, it may be possible to change only a few (or one) residues, or 
most of the residues, with all posslWfities In between. 

In a preferred embodiment, residues whidi can be fixed indude. but are not fimited to. structurally or 
biologically functional residues: alternatively, biologically fundional residues may spedfically not be 
fixed. For example, residues which are known to be Important for biological activity, such as the 
resMues v*ldi the binding site for a binding partner (Bgand/receptor. antigen/antibody, etc.). 
phosphorylation or glycosylation sftes which are crudal to biological fundion. or strudurally important 
residues, sudi as disuWde bridges, metal binding sites, critical hydrogen bonding residues, residues 
critical for backbone conformation sudi as Foline orglydne. residues critical for paddng interactions, 
etc. may all be fixed in their amino add Identity and a single rotamer conformation, or floated". wWch 
only fixes the identity but not the rotamer confomwtion. 

Similariy. residues whldi may be diosen as variable residues may be those that confer undesirable 
biological attributes, such as susceptibility to proteolytic degradation, dimerization or aggregation sites, 
glycosylation sites whidi may lead to immune responses, unwanted binding activity, unwanted 
allostery . undesirable enzyme activity but with a presereaBon of bindng. etc. 

In a preferred embodiment, each variable position Is dassified as either a core, surface or boundary 
residue position, although In some cases, as explained below, the variable position may be set to 
glydne to minimize badcbone strain. In addition, as outfined herein, residues need not be dassified. 
they can be diosen as variable and any set of amino adds may be used. Any combination of core, 
surface and boundary positions can be ufiBzed: core, surface and boundary resWues; core and surface 
residues: core and boundary residues, and surface and boundary residues, as well as core residues 
alone, surfece residues alone, or boundary residues alone. 
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The classification of residue positions as core, surface or boundary may be done in several ways, as 
will be appreciated by those in the art. In a preferred embodiment the classification is done via a 
visual scan of the ori^nat protein bac1<bone structure, including the side chains, and assigning a 
classification based on a subjective evaluation of one skilled in the art of protein modeOng. 
/Vltematively . a preferred embodiment uUBzes an assessment of the orientation of the Ca-CP vectors 
relative to a solvent accessible surface computed using only the template Ca atoms, as outDned in 
U.S.S.N5 60/061.097. 60/043.464. 60/054.678. 09/127.926 60/104.612. 60/158.700. 09/419.351. 
60/181630. 60/186.004. U.S patent application, entilled Protein Design Automation For Protein 
Ubraffes (Filed: /^pril 14. 2000; Inventor. Bassil Dahiyat) and PCT US98/07254. /Utemativety. a 
surfece area calculation can be done. 

Suitable core, boundaiy and surface positions for lA proteins are outlined below. 

Once each variable position is classified as either core, surface or boundary, a set of amino add side 
chains, and thus a set of rotameis. is assigned to each position. That is. the set of possible amino add 
side chains that the program will allow to be considered at any particular position is chosen. 
Subsequently, once the possible anino add side drains are chosen, the set of rotamers that win be 
evaluated at a particular position can be detennined. Thus, a core residue wiU generally be selected 
from the group of hydrophobic residues consisting of alanine, valine, isoleudne. leudne. 
phenylalanine, tyrosine, tryptophan, and methionine On some embodiments, when the a scafing factor 
of the van der Waals scoring function, described below, is low. methionine is removed from the set), 
and the rotamer set for each core position potentially indudes rotamers for these eight amino aad side 
chains (all the rotamers if a backbone independent Ibrary Is used, and subsets if a rotamer dependent 
backbone is used). SInilarty. surface positions are generally selected from the group of hydrophrbc 

residues consisting of alanine, serine, threonine, aspartic add, asparagine. glutamine. glutamic aad. 

aiglnine. lysine and histidine. The rotamer set for each surface position thus includes rotamera for 

these ten residues. Rnally. boundary positions are generally diosen from alanine, senne. threonine. 

aspartic add. asparagine. glutamine. glutamic add. arginine. lysine histidine. vattne. Isoleudne. 

leudne. phenylalanine, tyrosine, tryptophan, and methionine. The rotamer set for ead, boundary 

position thus potentially Includes every rotamer for these seventeen residues (assurring cysteine. 

glydne and proflne are not used, although they can be). Additionally, in some preferred embodiments. 

a set of 1 8 naturally occurring amino adds (aM except cysteine and proOne. whid, are known to be 

particularty disrupth/e) are used. 

Thus as will be appredated by those in the art. there Is a computational benefit to classifying the 
residue positions, as It decreases the number of calculations. It should also be noted that there may 
be situations where the sets of core, boundaiy and surface residues are altered from those described 
above: for example, under some drcumstances. one or more amino adds Is either added or 

16 



wo 00/69901 



PCT/USOO/13764 



subtracted from the set of allowed anvno adds. For example, some proteins which dimerize or 
multimerize. or have Ggandlnnding sites, may contain hydrophotnc surface residues, etc. In addition, 
residues that do not allow helix ''capping* or the favorable interaction with an a>helix dipole may be 
subtracted from a set of allowed residues. This modification of amino add groups is done on a 
residue by residue basis. 

In a prefenred en^odiment, proline, cysteine and gtydne are not Induded in the list of possible anruno 
add side chains, and thus the rotamers for these side chains are not used. However, in a preferred 
embodiment, when the variable residue position has a ^ angle (that is. the dihedral angle defined by 1) 
the carbonyl carbon of the preceding amino add; 2) the nitrogen atom of the current residue; 3) the a- 
carbon of the current residue; and 4) the carbonyl carbon of thjB cun-ent residue) greater than 0*, the 
position is set to glydne to minimize backbone strain. 

Once the group of potential rotamers is assigned for each variable residue position, processing 
proceeds as outlined in U.S.S.N. 09/127.926 and PCT US98/07254. This processing step entails 
analyzing interactions of the rotamers with each other and with the protein backbone to generate 
optimized protein sequences. Simplistically, the processing initially comprises the use of a number of 
scoring functions to calculate energies of interactions of the rotamers, either to the backbone itself or 
other rotamers. Prefen-ed PDA scoring functions include, but are not limited to, a Van der Waals 
potential scoring function, a hydrogen IxDnd potential scoring function, an atomic solvation scoring 
function, a secondary structure propensity scoring function and an electrostatic scoring function. As is 
further described below, at least one scoring function is used to score each position, although the 
scoring functions may differ depending on the position dasslfication or other considerations, like 
favorable interaction with an a-hefix dipole. As outlined below, the total energy wWch is used in the 
calculations is the sum of the energy of each scoring function used at a particular pos'ttion, as is 
generally shown in Equation 1 : 

Equation 1 

Eteta, = nE^ + nEa. + nEh4«^ + nE« + nEetec 

In Equation 1, the total energy is the sum of the energy of the van der Waals potential (E^^J, the 
energy of atonic solvation (E J, the energy of hydrogen bonding (Eh.bonrfing). ^le energy of secondary 
structure (E J and the energy of electrostatic interaction (^^, The tenn n is either 0 or 1 , depending 
on whether the term Is to be considered for the particular residue position. 

As outlihed In U.S.S.N.S 60/061 .097, 60/043,464, 60/054,678, 09/127,926, 60/104,612. 60/158,700, 
09/419,351, 60/181630. 60/186,904, U.S patent appKcation. entitled Protein Design Automation For 
Pmtein Libraries (Bled: April 14. 2000; inventon Bassil Dahiyat) and PCT US98A)7254, any 
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combination of these scoring functions, either alone or in comt)ination. nnay be used. Once the 
scoring functions to be used are identified for each variable position, the preferred first step In the 
computational analysis comprises the determination of the interaction of each possible rotamer with all 
or part of the remainder of the protein. Thai is, the energy of Interaction, as measured by one or mote 
of the scoring functions, of each possible rotamer at each variable residue position with either the 
backbone or other rotamers. Is calculated. In a preferred embodiment, the interaction of each rotamer 
with the entire remainder of the protein, I.e. both the entire template and alt other rotamers, is done. 
However, as outlined above, it Is possible to only nrwdel a portion of a protein, for example a domain of 
a larger protein, and thus in some cases, not ail of the protein need be considered. The tenn "portion", 
or similar grammatical equivalents thereof, as used herein, with regard to a protein r^ers to a fragment 
of that protein. This fragment may range in size from 6-10 amino acid residues to the entire amino 
add sequence minus one amino acid. Accordingly, the terni "portion", as used herein, with regard to a 
nucleic refers to a fragment of that nucleic acid. This fragment may range In size from 6-10 
nucleotides to the entire nudeic acid sequence minus one nucleotide. 

In a preferred embodiment, the first step of the computational processing is done by calculating two 
sets of interactions for each rotamer at every position: the interaction of the rotamer side chain with the 
template or backbone (the "singles" energy), and the interaction of the rotannerside chain with aU other 
possible rotamers at every other position (the "doubles" energy), v^rtiether that position is varied or 
floated. It should be undeistood that the backbone In this case includes both the atonre of the protein 
structure backbone, as well as the atoms of any fixed residues, wherein the fixed residues are defined 
as a particular conformation of an annino add. 

Thus, -singles" (rotamer/template) energies are calculated for the Interaction of every possible rotamer 
at every variable residue position with the backbone, using some or all of the scoring functions. Thus, 
for the hydrogen bonding scoring function, every hydrogen bonding atom of the rotamer and every 
hydrogen bonding atom of the backbone is evaluated, and the Ehb is calculated for each possible 
rotamer at every variable position. Sirrilarty. for the van der Waals scoring fundion, every atom of the 
rotamer is compared to every atom of the template (generally exduding the backbone atoms of Its 
own residue), and the E,^ is calculated for each possible rotamer at every variable residue position. 
In addition, generally no van der \Afeials energy is calculated If the atoms are conneded by three bonds 
or less. For tiie atorrilc solvation scoring fundion. the surface of the rotamer is measured aganst the 
surface of the template, and the E„ for each possible rotamer at eveiy variable residue position is 
calculated. The secondary strudure propenaty scoring fundion is also considered as a singles 
energy, and thus the total singles energy may contain an E« terot As will be appredated by those In 
the art. many of these energy terms will be close to zero, depending on the physical distance between 
the rotamer and the template position; that fe. the farther apart the two moieties, the lower the energy. 
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For the calculation of "doubles' energy (rotamer/rotanner), the interaction energy of each possible 
rotamer is compared with eveiy possible rotanner at all other variable residue positions. Thus, 
■doubles" energies are calculated for the interaction of every possible rotamer at every variable 
residue position with every possible rotamer at every other variable residue position, using some or all 
of the scoring functions. TTius. for the hydrogen bonding scoring function, every hydrogen bonding 
atom of the first rotamer and every hydrogen bonding atom of every possible second rotamer is 
evaluated, and the E„b is calculated for each possible rotamer pair for any two variable positions. 
Similarly, for the van der Waals scoring function, every atom of the first rotamer is compared to every 
atom of every possible second rotamer, and the is calculated for each possible rotamer pair at 
every two variable residue positions. For the atonnic solvation scoring function, the surface of the first 
rotamer is measured against the surface of every possible second rotamer, and the E„ for each 
possible rotamer pair at every two variable residue positions is calculated. The secondary structure 
propensity scoring function need not be mn as a "doubles' energy, as it is considered as a component 
of the "singles" energy. As will be appreciated by those in the art, many of these double energy terms 
will be close to zero, depending on the physical distance between the first rotamer and the second 
rotamer; that is, the farther apart the two moieties, the loyver the energy. 

In addition, as will be appreciated by those In the art, a variety of force fields that can be used in the 
PDA calculations can be used, including, but not limited to, Dreiding I and Dreiding II [Mayo et al. J. 
Phys. Chem. 94:8897 (1990)]. AMBER [Weiner et aL. J. Amer. Chem, Soc. 106:765 (1984) and 
Weiner et al. J. Comp. Chem. 106:230 (1986)], MM2 (Allinger. J. Chem. Soc. 99:8127 (1977), 
UljeforB et al.. J. Com. Chem. 8:1051 (1987)]; MMP2 ISprague et al.. J. Comp. Chem. 8:581 (1987)]; 
CHARMM [Brooks et aL, J. Comp. Chem. 106:187 (1983)]; GROMOS; and MM3 [Alfinger et al.. J. 
Amer, Chem. Soc. 1 1 1 :8551 (1989)]. OPLS-AA [Jorgensen et al.. J. Am. Chent Soc. 
118:11225-11236 (1996); Jorgensen. W.L.; BOSS, Version 4.1; Yale University: New Haven. CT 
(1999)]; OPLS [Jorgensen et al.. J. Am. Chem. Soc.110:1657ff (1988); Jorgensen et al., J Am. Chem. 
Soc. 1 12:4768ff (1990)]; UNRES (United Residue Forcefield; Uwo et al,. Protein Science 2:1697-1714 
(1993); Liwo et al., Protein Science 2:1715-1731 (1993); Liwo et al.. J. Comp. Chem. 18:849-873 

(1997) ; liwo et al., J. Comp, Chem, 18:874-884 (1997); Uwo et al., J. Comp. Chem. 19:259-276 

(1998) ; Fbtx»field for Protein Structure Prediction (Liwo et aL, Proc. Natl. Acad. Sd. U.SA 
96:5482-5485 (1999)]; ECEPP/3 [Liwo et aL. J Protein Chem. 13(4):375-80 (1994)]; AMBER 1.1 force 
field (Weiner. et al.. J. Am. Chem. Soc. 106:765-784); AMBER 3.0 force field (U.C. Singh et aL. Proc. 
Natl. Acad. Sd. U.S A. 82:755-759); CHARMM and CHARMM22 (Brooks et aL. J, Comp. Chem. 
4:187-217); cvfB.O [Dauber-Osguthorpe. et aL, Proteins: Structure, Function and Genetics. 4:31-47 
(1988)1; cff91 (Maple, et al,. J. Comp. Chem. 15:162-182); also, the DISCOVER (cvff and cf!91) and 
AMBER forcefields are used in the INSIGHT molecular modeOng package (BiosyrtVMSl. San Diego 
CaDfomia) and HARMM is used in the QUANTA molecular modeling package (Biosynn/MSI. San Diego 
California), all of which are expressly incorporated by reference. 
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Once the singles and doubles energies are calculated and stored, the next step of the computational 
processing may occur. As outlined in U.S.S.N. 09/127,926 and PCT US98/07254. preferred 
embodiments utilize a Dead End Elimination (DEE) step, and preferably a Monte Carlo step. 

PDA. viewed broadly, has three connponents that may be vaned to alter the output (e.g. the primary 
libraiy): the scoring functions used in the process; the filtering technique, and the sampling technique. 

In a preferred embodiment, the scoring functions may be altered. In a preferred embodiment, the 
scoring functions outlined above may be biased or weighted in a variety of ways. For example, a bias 
towards or away from a reference sequence or family of sequences can be done; for example, a bias 
towards wild-type or homolog residues may be used. Similarly, the entire protein or a fragment of It 
may be biased; for example, the active site may be biased towards wild-type residues, or domain 
residues towards a particular desired physical property can be done. Furthermore, a bias towards or 
against increased energy can be generated. Additional scoring funcUon biases include, but are not 
limited to applying electrostatic potential gradients or hydrophobidty gradients, adding a substrate or 
binding partner to the calculation, or biasing towards a desired charge or hydrophobidty. 

In addition, in an attemative embodiment, there are a variety of additional scoring functions that may 
be used. Additional scoring functions include, but are not fimited to torsional potentials, or residue pair 
potentials, or residue entropy potentials. Such additional scoring functions can be used alone, or as 
functions for processing the library after it is scored initially. For example, a variety of functions 
derived from data on binding of peptides to MHC (Major HistocompatibiUty Complex) can be used to 
rescore a fibrary in order to efirrtnate proteins containing sequences which can potentially bind to 
MHC, I.e. potentially immunogenic sequences. 

In a preferred embodiment, a variety of filtering techniques can be done, induding. but not Omited to. 
DEE and its related counterparts. Additional filtering techniques Indude, but are not Bnrtted to branch- 
and-bound techniques for finding optimal sequences (Gordon and Mayo. Structure Fold. Des. 7:1089- 
98, 1999), and exhaustive enunreration of sequences. 

As will be appredated by those in the art. once an optimized sequence or set of sequences is 
generated, a variety of sequence space sampfing metiiods can be done, eittierin addition to the 
preferred Monte Cario methods, or instead of a Monte Cario search. That Is, once a sequence or set 
of sequences is generated, preferred methods utifize sanpling techniques to allow the generation of 
additional, related sequences for testing. 

These sampDng rnetfiods can include ttie use of anino add substitutions, insertions or deletions, or 
recombinations of one or more sequences. As outOned herein, a prefen-ed embodiment utilizes a 



20 



wo 00/69901 



PCTAJSOO/13764 



Monte Cario search, which is a series of laased. systematic, or random jumps. However, there are 
other sampfing techniques that can be used, including Boltzman sampUng, genetic algorithm 
techniques and simulated anneaOng. In addition, for all the sampling techniques, the kinds of jumps 
allowed can be altered (e.g. random jumps to random residues, biased jumps (to or away fixim wild- 
type, for example), jumps to biased residues ^o or away from slnralar residues, for example, etc.). 
Jumps where multiple residue positions are coupled (two residues always change together, or never 
change together), jumps where whole sets of residues change to other sequences (e.g.. 
recombination). Sintlarty, for all the sampling techniques, the acceptance criteria of whether a 
sampfing jump is accepted can be altered. 

In addition, it should be noted that the preferred methods of the invention result in a rank ordered fist of 
sequences; that Is, the sequences are ranked on the basis of some objective criteria. However, as 
outBned herein, it is possible to create a set of non-ordered sequences, for example by generating a 
probabifity table directly (for example using SCMF analysis or sequence alignment techniques) that 
lists sequences without ranking them. The sampling techniques outlined herein can be used in either 
situation. 

In a preferred embodiment. Boltzman sampOng is done. As wiil be appreciated by those In the art, the 
temperature criteria for Boltzman sampfing can be altered to allow broad searches at high temperature 
and narrow searches dose to local optima at low temperatures (see e.g.. Metropolis et al., J. Chem. 
Phys. 21:1087, 1953). 

In a preferred embodiment, the sampling technique utifizes genetic algorithms. e.g., sudi as those 
described by Holland (Adaptation in Natural and Artifical Systems, 1975. Ann Art)or, U. Michigan 
Press). Genetic algorithm analysis generally takes generated sequences and recombines them 
computationally, similar to a nucleic add recombination event, in a manner similar to "gene shuffling". 
Thus the -jumps" of genetic algorithm analysis generally are multiple position jumps. In addition, as 
outlined below, correlated multiple jumps may also be done. Such jumps can occur with different 
crossover positions and more than one recombination at a time, and can involve recombination of two 
or more sequences. Furthermore, deletions or inseitions (random or biased) can be done. In 
addition, as outlined below, genetic algorithm analysis may also be used after the secondary libraiy 
has been generated. 

In a prefen-ed embodiment, the sampling technique utilizes simulated anneafing. e,g., such as 
described by Wrkpatricketal. [Science. 220:671-680 (1983)]. Simulated anneafing alters the cutoff for 
accepting good or bad jumps by altering the temperature. That is. the stringency of the cutoff is 
altered by altering the temperature. This allows broad searches at high temperature to new areas of 
sequence space, altering with narrow searches at low temperature to explore regions in detail. 
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In addition, as outfined below, these sampBng mettiods can be used to furtlier process a first set to 
generate additional sets of lA proteins. 

The computationa! processing results in a set of optinized lA protein sequences. These optimized lA 
protein sequences are generally significantly different from the wild type insulin sequence from which 
the backbone was taken. That is. each optimized lA protein sequence preferably comprises at least 
about 2-50% variant anino adds from the starting or wild type sequence, with at least about 25% 
being preferred, with at least about 15-20% changes being more preferred and at least 4-15% being 
pailiculariy preferred. 

In a preferred embodiment, the lA proteins of the invention have 1. 2. 3. 4. 5, 6. 7, 8. 9. 10. 11. 12. 13. 
14. 15, 16, 17. 18. 19. 20, 21. 22. 23. 24. 25. 26. 27. 28. 29, or 30 different resWues from the human 
insulin sequence. 

In one aspect of this embodiment, the lA protein of the invention has at least one different residue 
from the human insulin sequence. Preferred lA protein sequences comprising a substitution of one 
amino acid resldye are shown In Figures 4B. 4C. 40. 4F, and 4G. 

In another aspect of this embodiment, the lA protein of the invention has at least two different residues 
from the human insufin sequence. Preferred lA protein sequences comprising a substitution of two 
anino acid residues are shown in Figures 3C, 3F. 4A, and 4E. 

In another aspect of this embodiment, the lA protein of the Invention has.at least three different 
reslduesfromthehumanlnsullnsequence. a preferred lA protein sequence comprising a substtuUon 

of three amino add residues is shovm in Figure 3E. 

in another asped of this embodiment, the lA protein of the invention has at least four different 
residues from the human insuDn sequence. Preferred lA protein sequences comprising a substitution 
of four artno add residues are shown In Rgures 3C. 3F. 4A. and 4E. 

in another asped of ttiis embodiment, the lA protein of the invention has at least five different residues 
from «ie human insulin sequence. 

In another asped of this embodiment, the lA protein of the Invention has at least six different residues 
from the human insulin sequence. Aprefened lA protein. sequence comprising a substitution of six 
airino add residues is shown in Figure 5A. A preferred lA protein sequence comprising a subshhiton 
of two amino add residues and a deletion of four amino add residues is shown in Figure 3G. 
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In another aspect of this emtradiment. the lA protein of the invention has at least seven different 
residues from the human insulin sequence. 

In another aspect of this embodiment, the lA protein of the invention has at least eight different 
residues from the human insulin sequence. 

In another aspect of this emlsodiment, the lA protein of the invenUon has at least nine different 
residues from the human insulin sequence. 

In another aspect of this embodiment, the I A protein of the invention has at least ten different residues 
from the human Insulin sequence. 

In another aspect of this embodiment, the I A protein of the invention has at least eleven different 
residues from the human insulin sequence. 

In another aspect,of this embodiment the lA protein of the invention has at least twelve different 
residues from the human insulin sequence. 

In another aspect of this embodiment, the lA protein of the invention has at least thirteen different 
residues from the human insulin sequence. 

In another aspect of this embodiment, the lA protein of the Invention has at least fourteen different 
residues from the human insufin sequence; A preferred lA protein sequence comprising a substitution 
of fourteen amino acid residues is shown in Figure 58. 

in another aspect of this embodiment, the lA protein of the invention has at least fifteen different 
residues from the human insulin sequence. 

In another aspect of this embodiment, the lA protein of the invention has at least sixteen different 
residues from the human insuUn sequence. A preferred lA protein sequence comprising a substitution 
of sixteen amino acid residues is shown in Figure 5C. 

In another aspect of this embodiment the lA protein of the invention has at leait twenty different 
residues from the human insuOn sequence. A prefen-ed lA protein sequence comprising a substitution 
of twenty amino add residues Is shown In Rgure 3A. 

Thus, in the.broadest sense, the present invention is directed to lA proteins that have insuBn activity. 
By Insulin activity" or "lA" herein is iroant that the lA protein exhibits at least one. and preferably 
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more, of the biological functions of an Insulin, as defined below. In one embodiment, the biological 
funcUon of an lA protein is altered, preferably Improved, over the corresponding biological activity of an 
insufin. 

By 'protein' herein is meant at least two covalently attached amino adds. wWch indudes proteins, 
polypeptides. oOgopepBdes and peptides. The protein may be made up of naturally occuning amno 
adds and peptide bonds, or synthetic peptidonimetic stnidures. I.e.. -analogs" sudi as peptdds [see 
Simon et al.. Proc. Natl. Acd. Sd. U.S A 89(20:9367-71 (1992)1. generally depending on the method 
of synthesis. Thus "amino add", or 'peptide residue', as used herein means both naturally occurring 
and synthetic amino adds. For example, homo-phenylalanine. dtrulfine. and noreleudne are 
considered airino adds for the purposes of the invenUon. "Amino add" also indudes inino add 
residues such as profine and hydroxyproBne. In addition, any amino add representing a component of 
the lA proteins can be replaced by the same amino add but of the opposHe diirality. "mus. any amino 
add naturally occurring In the L-configuraBon (whidi may also be refeired to as the R or S. depending 
upon the stnicture of the chemical entity) may be replaced with an amino add of the same chemical 
structural type, but of the opposite diirality. generally referred to as the D- amino acid but which can 
additionally be referred to as the R- or the S-. depending upon its composition and chemical 
configuration. Such derivatives have the property of greatly increased stabffity. and therefore are 
advantageous in the formulation of compounds whi* may have longer In vWo half lives, when 
administered by oial. Intravenous. Intramuscular, intraperitoneal, topical, redal. intraocular, or other 
routes. In the preferred embodiment the amino adds are In the (S) or L-conHguraUon. Ifnon- 
naturelV occurring side diains are used, non-amino add substituents may be used, for example to 
prevent or retell in vivo degradations. Proteins induding non-naturally occurring amino adds may be 
synthesized or in some cases, made recombinantly; see van Hest et al.. FEBS Lett 428:(1-« 68-70 
May 22 1998 and Tang et al.. Abstr. Pap Am. Chem. S218:U138-U138 Part 2 August 22. 1999. both of 
which are expressly incorporated by reference herein. 

Additionally, modified amino adds or diemical derivatives of amino adds of consensus or fragments of 
lA proteins. acconJIng to the present Invention may be provided, which polypeptides contain additional 
diemical moledes or modified amino acids not nomially a part of the protein. Covalent and non- 
covalent modifications of the protein are thus induded within the scope of the present invention. Such 
modifications may be introduced into an lA polypeptide by reading targeted amino add residues of the 
polypeptide with an organic derivatizing agent that is capable of reading with selected side dia^^^ 

temtinal residues. The following examples of diemical derivatives are provided by way of niustration 

arxl not by wray of fimitation. 
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Aromatic amino acids may be replaced with D- or L-naphytalanine, D- or L-Phenylglycine, D- or L-2- 
thieneylalanlne. 0- or L-1-. 2-. 3- or 4-pyreneylalanine, D- or L-3-thleneylalanine. D- or L-(2-pyrid'myl)- 
alanirje, D- or L-(3-pyridinyl)-alanine, D- or L-(2-pyrazinyl)-alanine. D- or L-(4-isopropyl)-phenyl9lycine. 
0-(trifluoromethyl)-phenylglycine. D-(tfifluorQmethy«)-phenylalanine. D-p-fluorophenylalanine. D- or L- 
p-biphenylphenylalanine. D- or L-p-methoxybiphenylphenylalanine. D- or L-2-indole(aIlcyl)alanines, and 
D- or L-alltylainlnes where alkyi may be substituted or ansubstituted methyl, ethyl, propyl, hexyl, butyl, 
pentyl. isopiDpyl, Iso-butyl, sec^sotyl, Iso-pentyl. noh-acidic antino adds, of C1-C20. 

Acidic amino acids can be substituted with non-caitoxylate amino acids while maintaining a negative 
charge, and derivatives or analogs thereof, such as the non-limiting examples of (phosphono)alanine, 
(phosphono)glycine. (phosphono)leucine. (phosphono)isoleudne. (phosphono)threonine, or 
{phosphono)serine; or sulfated (e.g., -SOjH) threonine, serine, tyrosine. 

Other substitutions may include unnatural hyioxylaled amino adds that may be made by combming 
"alkyI- wtth any natural amino acljJ. The term "allcyr as used herein refers to a branched or 
unbranched saturated hydrocarbon group of 1 to 24 carbon atoms, such as methyl, ethyl, n-propyl. 
isoptopyl, n-butyl, Isobutyl, t-butyl. octyl. decyl. tetradecyl. hexadecyt. eicosyl. tetracisyl and the like. 
Preferred alkyI groups herein contain 1 to 12 carbon atoms. Also induded within the definition of an 
aikyi group are cycloalkyi groups sudi as C5 and C6 rings, and heterocydic rings with nitrogen, 
oxygen, sulfur or phosphoms. AlkyI also indudes heteroalkyi, with heteroatoms of suHiir, oxygen, and 
nitrogen being preferred. AlkyI includes substituted alkyI groups. By -subsBhrted alkyi group" herein Is 
meant an alkyI group further comprising one or more substitution moieties. A preferred heteroalkyi 
group is an aBcyl amine. By "alkyl amine" or grammatical equivalents herein is meant an aikyi group as 
defined above, substituted with an amine group at any position. In addition, the alkyI amine may have 
other substitution groups, as outOned above for alkyI group: The amine may be primary (-NH/1), 
secondary (-WHRJ. or tertiary (-NR,). Basic amino adds may be substituted with alkyI groups at any 
position of the naturally occumng amino adds lysine, arginine, omrthine. citrulBne, or (guanldino)- 
acetic add. or other (guankJino)alky»-acetic acids, where Wis define as above. Nitrite derivatives 
(e.g.. containing the CN-mofety In place of COOH) may also be substituted for asparaglne or 
glutamine, and methionine sulfoxide may be substituted for methionine. Methods of preparation of 
such peptWe derivatives are well known to one skilted In the art 

In addltton, any amide linkage In aiiy of the lA polypeptides can be replaced by a ketomethylene 
moiety. Such derivatwes are expeded to have the property of increased stability to degradation by 
enzymes, and therefore possess advantages for the formulatton of compounds which may have 
increased in vivo half lives, as administered by oral. Intravenous, intramuscular, intraperitoneal, topical. 
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rectal, intraocular, orotherroutes. 

Additiona! amino add nfiodtfications of amino adds of [A polypeptides of the present Invention may 
indude the following: Cysteinyl residues may be reacted with alpha-haloacetates (and conBsponding 
amines), such as 2-chloroacetic acid or chloroacetamide, to give carboxymethyl or 
carboxyamidomethyl derivatives. Cysteinyl residues may also be derivatized by reaction with 
compounds such as bromotrifluoroacetone, alpha-brom6-beta-(5Hmidozoyl)proplonlc acid, 
chloroacelyl phosphate, N-alkylmaleimldes, 3-nitrD-2-pyridyt disulfide, methyl 2-pyridyl disulfide, p- 
chloromercuribenzoate, 2-chloromercurHl-nitiophenol, or chlort>-7-nitrobenzo-2-oxa-1.3-diazDle. 

Histidyl residues may be derivatized by reaction with compounds such as diethylprocarbonate e.g.. at 
pH 5.5-7.0 because this agent is relatively specific for the histidyl side chain, and para-bromophenacyl 
bromide may also be used; e.g., where the reaction is preferably performed In 0.1M sodium 
cacodylate at pH 6.0. 

Lysinyl and amino terminal residues may be reacted with compounds such as sucdnic or other 
carboxylic acid anhydrides. Derivatization with these agents is expected to have the effect of reversing 
the chaige of the lysinyl residues. Other suitable reagents for derivatizing alpha-amino-containing 
residues indude compounds such as imidoesters/e.g,. as methyl picolinimidate; pyridoxal phosphate; 
pyridoxal; chloroborohydride; trinitrobenzenesulfonic add; O-methylisourea; 2.4 pentanedione; and 
transaminase-catalyzed reaction with glyoxylate. 

Arginyl residues may be modified by reaction with one or several conventional reagents, among them 
phenylgiyoxal. 2,3-butanedione. 1,2-cyctohexanedk)ne, and ninhydrin according to Icnown method 
steps, Derivatization of arginine residues requires that the reaction be perfomned in alkaline conditions 
because of the high pKa of the guanidine functional group. Furthemiore. these reagents may react 
with the groups of lysine as well as the arginine epsilon-amino group. 

The spedfic modificatton of tyrosyl residues per se Is well-known, such as for Introdudng spectral 
labels Into tyrosyl residues by reactton with aromatic dtazonium compounds ortetranitronn^ane. N- 
acetylln^Jizol and t^nltromethane may be used to fomn 0-acetyl tyrosyl spedes and 3-^Utro 
derivatives, respectively. 

Cart)Oxyl side groups (aspartyl or glutamyl) may be selectively modified by reacUon with cartxxJiintides 
(R-N-C-N-R') such as 1-cydohexyl-3-(2-morpholinyl- (4-ethyl) carbodlimide or 1.ethyl-3-<4-azonia-4,4- 
dimethylpentyl) cart)odiimide. Furthemiore aspartyl and glutamyl residues may be converted to 



26 



wo 00/69901 



PCTAJSOO/13764 



asparaginyl and glutaminyl residues by reaction wHh ammonium tons. 

Glutaminyl and asparaginyl residues may be frequently deamidated to the corresponding glutamyl and 
aspartyl residues. Alternatively, these residues may be deamidated under ntildly acidic conditions. 
Either torn of these residues falls within the scope of the present invention. 

The Insulin may be from any nun*er of organisms, with insuOns from mammals being particularly 
preferred. Suitable mammals Include, but are not finnted to. rodents (rats, mice, hamsters, guinea 
p^s, etc.). primates, fami aramals Cmcluding sheep, goats, pigs. cows, horses, etc) and in the most 
preferred embodiment, from humans ( the sequence of which is depicted in Figure 1). As will be 
appreciated by those in the art. InsuHns based on insuOns from mammals other than humans may find 
use in animal models of human disease. The GenBank accession numbers for a variety of 
mammalian insuHn species is as follows: bovine. IPBO; dog. IPDG; sheep. INSH; cat. INCT; pig. 
IPPG; mouse. INMS1. INMS2: rat, IPRT1. 1PRT2; horse. IPHO; rabbit. INRB; guinea pig. IPGP; 
hamster. INHY; goat, INGT. chimpanzee. A421 79; green monkey, 842179; and human IPHU. 

The lA proteins of the invention exhibit at least one Wotogical function of an insulin. By InsuHn" or 
herein is meant a wild type insulin, an allefic variant thereof or a hybrid fomied from a combination of 
an A-chain of one spedes and a B-chain of another spedes. Thus, insufin refers to all fbm^ of insulins 
that are active In accepted insulin assays. 

The lA proteins of the invention exhibit at least one biological function of an insuOn By -Wotogical 
function- or 'biological property" herein is meant any one of this Foperties or functions of an InsuBn. 
induding. but not Hmited to, the abiOty to bind to a naturally occurring insulin receptor; the abifity to bind 
to a recomWnanl insuOn recepton the ability to bind to a cell comprising an insulin receptor, the ability 
to Induce autophosphorylation of an insuOn receptor [Combettes^ouverain and Issad. Diabetes 
Metab. 24(6):477-89 (1998)]; the ability to regulate a Na+/K=-ATPase (Sweeney and KUp, Mol. Cell. 
Biodiem. 182(1-2):121^3 (1998)J; the ability to regulate the Ras adivationfinactivation cyde [Ceresa 
and Pessin. Mol. Cell. Biochem. 182(1-2):2« (1998)1; the aWBty to treat type I diabetes, the ability to 
treat type 2 diabetes; the ablBty to treat an insulin-dependent dlsorden the aWWy to treat 
hyperglycerta; the abWty to treat familial hyperproinsuBnemia; the ability to regulate carbohydrate 
metaboDsm; the abiOty to regulate skeletal muscle protein turnover (Grizard et al., Reprod. Nutr. Dev. 
39(1):61-74 (1999)1; the abiBty to modulate intracellular processes Induding, but not Bnfted to cellular 
metaboBsm,,cell proBferation, and cell differentiation IRizzo and Romero, J. Basic COn. Physiol. 
Pharmacol 9(2-4):167-95 (1998)1; the abiBty to form muffimers; the abiBty to complex anc; and the 
ability to bind a phenoBc presenralive. 
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All of these I A proteins will exhibit at least 20% of the receptor binding or Uologlcal activity as the wild 
type insulin. More preferred are lA proteins that exhiUt at least 50%, even more prefenred are I A 
proteins that exhibit at least 90%, and nwst preferred are lA proteins that exhibit more than 100% of 
the receptor binding or biotogica! activity as the wild type insulin. Biological assays, receptor binding 
assays, and other assays for measuring insulin activity, stability, conformation, potency, etc. are 
described e.g., in the Examples, in US patents 4,421 .685; 4.992,417; 5,008.241; 5,506.202; 
5,514.646; 5,559,094; 5,618.913; 5.621.073; 5,663,291; 5.700.662; and 6,034.054 and further by. 
Marld et al. [Hoppe Seylers 2. Phsiol. Chem. 360(1 1):1619-32 (1979)]; Hu et al. [Biochemistry . 
32(10):2631-5 (1993)1; Schwartz et al. [Proc. Natl. Acad. Sd. U.SA 84(18):6408-11 (1987)]; Kitagawa 
et ai. [Biocherristry 23(7):1405-13 (1984)]; Kobayashi et al. [Biochem. Biophys. Res. Commun. 
107(1):329-36 (1982)]; Shoelson et al. [Biochemistiy 31(6):1757-67 (1992)1; Hua et al. [Nature 
354:238-241 (1991)]. Bao et al. [Proc. Natl. Acad. Sd. U.SA 94:2975-2980 (1997)1; Kitagawa et al. 
[Biocherrastry 23:1405-1413 (1984)1; Nakagawa and Tager [Biochemistiy 31:3204-3214 (1992)1; 
Brader et al. (Biocheinstiy 30:6636-6645 (1991)]; and references dted therein, all of which are 
expressly incorporated as reference. 

In one embodiment, at least one biological property of the lA protdn is altered when connpared to the 
same property of insulin. As outfihed above, the invention provides lA nucleic adds encoding lA 
polypeptides. Ttie lA polypeptWe pr^erably has at least one property, which is substantially different 
from the same property of the corresporiding naturally occuning insulin. The property of the lA 
polypeptide is the result . the PDA analysis of the present invention. 

The term "altered property" or grammatical equivalents thereof In the context of a polypeptide, as used 
herein, refer to any characteristic or attribute of an lA polypeptide that can be selected or detected and 
compared to the corresponding property of a naturally occurring Insulin protein. These properties 
Indude. but are not lintited to oxidath/e stability, prolonged shelf-life; themial stability, alkaline stability, 
pH activity profile, resistance to proteolytfc degradalfan. kinetic assodatfon (K«J and dissodatfon (K^) 
rate, protein folding, Indudng an immune response, ability to bind to an insulin receptor, ability to be 
secreted, ability to oligomerize, ability to form dimers. ability to form tetramers, ability to fomri 
hexamers. ability to complex zinc, ability to bind phenolic preservative, ability to modulate potency, 
preference for the R state, preference for the T state, ability to signal, abfltty to modulate cell 
proliferation, ability to be modified by phosphorylation or glycosylation. ability to treat hypeiplycenrtf a. 
ability to treat diab^es mellitus. abaity to treat type 1 diabetes. abiTtty to treat type 2 diabetes, ability to 
treat an insulin-dependent dfeorder, ability to treat an insulin-independent disorden ablTrty to achieve 
hormonal homeostasis in a patient; ability to be cleared from the blood of a diab^c patient. 
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Unless otherwise specified, a substantial change in any of the above-listed properties, when 
comparing the property of an lA polypeptide to the property of a naturally occurring insulin protein is 
preferably at least a 20%. more preferably. 50%. more pr^erably at least a 2-fbld IncrBase or 
decrease. 

A change in oxidative stability is evidenced by at least about 20%. more preferably at least 50% 
increase of activity of an lA protein when exposed to various oxidizing conditions as compared to that 
of insulin. Oxidative stabifity is measured by known procedures. 

A change in alkaline stability is evidenced by at least about a 50% or greater increase or decrease 
(preferably Increase) in the half life of the activity of an lA protein when exposed to increasing or 
decreasing pH conditions as compared to that of insulin. Generally, alkaline stability is measured by 
known procedures. 

A change in thermal stability is evidenced by at least about a 50% or greater increase or decrease 
(pr^erably increase) in the half life of the activity of an lA protein when exposed to a relatively high 
temperature and neutral pH as compared to that of insulin. Generally, themnal stability is measured by 
known procedures. In a preferred embodiment, an lA protein of the invention has increase themnal 
stability when compared to the human insulin. Such an lA protein preferably has an amino acki 
sequence which comprises substitution of one or more amino acid residues when compared to the 
amino acid sequence of human insulin. In one aspect of this embodiment, the amino ackJ sequence of 
the lA protein comprises 2, 3, 4. 5. 6. 7. 8. 9, 10, or more substltutk>ns when compared to the amino 
acid sequence of human Insulin. 

Similarly, lA proteins, for example are experimentally tested and validated in in vivo and in in vitro 
assays. Suitable assays include, but are not limited to, e.g.. examining their binding affinity to natural 
occurring or variant receptors and to high affinity agonists and/or antagonists. In addition to cell-free 
biochemical affinity tests, quantitative comparison are made comparing kinetic and equilibrium binding 
constants for the natural receptor to the naturally occurring insulin and to the lA proteins. The kinetic 
association rate (Kon) and dissociation rate (J<^), and tiie equilibrium binding constants (Kd) can be 
determined using surface plasmon resonance on a BIAcore instrument fbltowlng tiie standard 
procedure in tiie Rterature pearce et al., Biochemistiy 38:81-89 (1999)]. Comparing the binding 
constant between a natural receptor and its corresponding naturally occurring insulin with the binding 
constant of a natural occurring receptor and an lA protein are made in order to evaluate the sensitivrty 
and specificity of the lA protein. Preferably, binding affinity of tiie lA protein to natural receptors and 
agonists increases relative to the naturally occurring insulin, white antagonist affinity decreases. lA 
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proteins with higher affinity to antagonists relative to the Insulin may also be generated by the methods 
of the invention. 

In one preferred embodiment, an lA protein of the invention has modulated potency when compared to 
the human insulin. The term "modulated potency" herein means that an lA protein leads to an 
increase (more potent) or a decrease (less potent) in response when compared to the human insulin. 
The response includes any response due to the biological property of an lA protein as defined above. 
In one aspect of this embodiment, the modulated potency is brought about by a faster acting lA 
protein, a slower acting lA protein, a longer acting lA protein, a shorter acting lA protein, or causing a 
similar response as the human insulin at a lower concentration. Such an lA protein preferably has an 
amino acid sequence which comprises substitutions of four or more amino acid residues when 
compared to the arrtino acid sequence of human insulin. In one aspect of this embodiment, the amino 
acid sequence of the lA protein comprises 5, 6, 7. 8. 9, 10, or more substitutions when compared to 
the amino acid sequence of human insulin. 

In one preferred embodiment, an lA protein of the invention shows increased hexamer fonmation when 
compared to the human insulin. The term "increased hexamer formation' herein means that an I A 
protein forms a hexamer at a lower concentration than the human insulin or that the hexamer fornied 
by an lA protein Is more stable than those fornied by human insulin. Such an lA protein preferably has 
an amino acid sequence which comprises substitution of one or more amino acid residues when 
compared to the amino acid sequence of human insulin. In one aspect of this embodiment, the amino 
acid sequence of the I A protein comprises 2, 3. 4. 5. 6. 7. 8, 9. 10. or more substitutions when 
compared to the amino acid sequence of human insulin. 

In one pr^erred embodiment, an lA protein of the invention has an Increased pr^erence for the R 
state when compared to the human insulin. The terni "increased preference for ttie R state" herein 
means a greater portion of an lA protein is In tiie R state than in tfie T state when conr^ared to human 
insulin. Such an lA protein preferably has an amino acid sequence which comprises substitution of 
three or more an^no acid residues when compared to tiie amino acid sequence of human insulin. In 
one aspect of this embodiment, tiie amino acid sequence of the lA protein comprises 4, 5. 6, 7, 8, 9, 
10, or more substitutions when compared to the antino acid sequence of human Insulin. 

As described In tiie art, both zinc and a phenoRc preservative are essential to achieve a complex that 
is stable aria capable of rapid dissociation and onset of action. The hexamer complex consists of two 
zinc ions per hexamer of insulin, and at least tiiree molecules of a phenolic preservative. "Phenolic 
preservative" as used herein refers to clorocresol. nr>-cresol, phenol, or mixtures ttiereof. In a 
preferred embodiment, an lA protein does not bind phenolic preservatives or binds less that insulin. 
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High insulin concentrations lead to ollgomerization. thefomration of dimers, tetramers and hexamers. 
In a preferred embodiment, an lA protein oligomerizes at a lower concentration as the corresponding 
wild type Insulirl. 

Insulin preparations for treatment of diabetes are administered in adequate therapeutic amounts of the 
active compound to achieve hormonal homeostasis. In a preferred embodiment, an lA protein is more 
potent than natural insulin in vitro. In one aspect of this embodiment, it is believed that the therapeutic 
amount of the lA protein to achieve homeostasis in diabetic patient or in animal mod^ is less than the 
therapeutic amount of the wild type insulin. Various animal models are available, for example the BB 
mouse INakbookda et a!., Diabetologia 14(3):199-207 (1978)] and the NOD (non-obese^iabetic) 
mouse in which diabetes develops spontaneously (Prochazka et al.. Science 237(4812):286-9 (1987)1. 

Insulin clearance from the blood is mediated by the insulin receptor on cells. In a preferred 
embodiment, an lA protein binds to an insulin receptor more tightly than wild type Insulin. In one 
aspect of this embodiment, it is believed that I A protein is cleared from the blood of a patient at a faster 
rate than wild type insulin. As a consequence thereof, it is further believed that in treatment of 
diabetes that vascular toxicity associated with the growth-promoting effects of circulating insuUn may 
be lessened by the use of lA proteins. 

As described above, one biological function of an lA protdn Is the abifity of the lA protein to bind to 
cells comprising an insulin receptor. GenBank Accession numbers for InsuOn receptors OF^ are 
available for various species, e.g.. human. INHUR. P06213. NP_000199; mouse. A34157. P15208; 
andrat.A36080.P15127. 

In a preferred embodiment, the assay system used to detenrone lA is an in vitro system using cells 
that either express endogenous insulin receptors or cells stably transfected with the gene encoding the 
hurnan insuBn receptor or an insulin receptor from another species, known to bind human Insulin, e.g., 
mouse. In this system, cell proBferation Is measured as a function of BrdU incorporation, which is 
Incorporated into the nucleic add of proOferating cells. A decrease above background of at least about 
10%, with at least about 20% being preferred, with at least about 30% being more prefered and at 
least about 50%. 75% and 90% being especially prefen-ed is an indication of lA., 

In another preferred embocfiment, autophosphoiylation of tfie InsuUn receptor, which occurs after 
insulin tending Is detemnined using known phosphorylation assays. 

In a preferred embodiment, the antigenic profile in the host animal of tiie lA protein is similar, and 
preferably identical, to the antigenic profile of the host InsuHn; that is. tfte lA protein does not 
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significantly stimulate the host organism (e.g. the patient) to an immune response; that is, any immune 
response Is not clinically relevant and there is no allergic response or neutrafization of the protein by 
an antibody. That is, in a preferred embodiment, the lA protein does not contain additional or different 
epitopes from the insiifin. By 'epitope" or "determinanr herein is meant a portion of a protein which will 
generate and/or bind an antibody. Thus, in most instances, no significant amount of antibodies are 
generated to a lA protein. In general, this is accomplished by not significantly altering surface 
residues, as outlined below nor by addii^ any amino acid residues on the surface which can become 
glycosylated, as novel glycosylatlon can result in an Immune response. 

The lA proteins and nucleic adds of the invention are distinguishable from naturally occurring insulins. 
By 'naturally occurring" or "wild type" or grammatical equivalents, herein is meant an amino acid 
sequence or a nucleotide sequence that is found in nature and includes allelic variations: that is, an 
amino acid sequence or a nucleotide sequence that usually has not been intentionally modified. 
Accordingly, by -non-naturally occurring" or "synth^c" or "recombinant" or grammatical equivalents 
thereof, herein is meant an amino acid sequence or a nucleotide sequence that is not found in nature; 
that is. an amino add sequence or a nucleotide sequence that usuaily has been intentionally modified. 
It is understood that once a recombinant nudeic add is made and reintroduced Into a host cell or 
organism, it will replicate non-recombinantly. I.e., using the In vivo cellular machinery of the host cell 
rather tharf/n vitro manipulations, however, such nucleic adds, once produced recombinantly, 
although subsequently repficated non-recombinantly. are still considered recombinant for the purpose 
of the invention. A representative amino add sequences of a naturally occurring human insufin is 
shown in Rgure 1 . it should be noted that unless othenvise stated, alt positional numbering of lA 
proteins and I A nucleic adds is based on these sequences. That is, as will be appredated by those In 
the art. an alignment of insuOn proteins and lA proteins can be done using standard programs, as is 
ouUined below, with the Identification of "equivalenf positions between tiie two proteins. Thus, the lA 
proteins and nucleic adds of the Invention are non-nalurally occurring; that is, they do not exist in 
nature. 

Thus. In a preferred embodiment, the lA protein has an amino add sequence that differs from a wild- 
type insuOn sequence by at least 2% of ttie residues. That is, the lA proteins of the invention are less 
than about 98% Identical to an insuPn amino add sequence. Accordingly, a protein is an "lA protein" if 
the overall homology of the protein sequence to the amino add sequence shown in Figure 1A or 
Figure 1B is preferably less than about 98%. more preferably less than about 95%. even more 
preferably less than about 90% and most preferably less ttian 85%. In some embodiments the 
homology wifl be as low as about 75 to 80%. In oti>er embodiments the homology will be as low 50- 
70%. Staled differently, based on the human insuOn sequence of 51 residues (see Rgure IB), lA 
proteins have at least about 1 residue that differs from the human insufin sequence (2%). with lA 
proteins having from 2 residues to upwards of 25 residues being different from the human insulin 
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sequence. Preferred I A proteins have 1-20 different residues with from about 2 to about 10 being 
particularly preferred (that is, 4-20% of the protein is not identical to human insulin). 

Homology in this context means sequence similarity or identity, with identity being preferred. As is 
known In the art, a hunr\ber of different programs can be used to identify whether a protein (or nucleic 
acid as discussed be^ow) has sequence Identity or similarity to a known sequence. Sequence Identity 
and/or similarity is d^emtined using standard techniques known in the ait. including, but not linruted to. 
the local sequence identity algorithm of Sn^th & Waterman, Adv. Appl. Math., 2:482 (1981). by the 
sequence identity alignment algorithm of Needleman & Wunsch, J, Mol. Biol.. 48:443 (1970), by the 
search for similarity method of Pearson & Lipman, Proc. Natl Acad. Sd. U.SA.. 85:2444 (1988). by 
computerized implementations of these algorithms (GAP. BESTFIT. FASTA. and TFASTA in the 
Wisconsin Genetics Software Package, Genetics Computer Group, 675 Science Drive. Madison, Wl), 
the Best Fit sequence program described by Devereux et aL. Nucl. Add Res., 12:387-395 (1984), 
preferably using the default settings, or by inspection. Preferably, percent identity is calculated by 
FastDB based upon the following parameters: mismatch penalty of 1 : gap penalty of 1: gap size 
penalty of 0.33; and joining penalty of 30, "Cun^ent Methods in Sequence Comparison and Analysis." 
Macromolecule Sequendng and Synthesis, Selected Methods and Applicatfons, pp 127-149 (1988), 
Alan R. Liss. Inc. 

An example of a useful algorithm is PILEUP. PILEUP creates a multiple sequence alignment from a 
group of related sequences using progressive, paiiwise alignments. It can also plot a tree showing the 
dustering relationships used to create the alignment PILEUP uses a simplification of the progressive 
alignment m^hod of Feng & Doolitlle, J. Mol. Evol. 35:351-360 (1987); the method Is similar to that 
described by Higgins & Sharp CABIOS 5:151-153 (1989). Useful PILEUP parameters induding a 
default gap weight of 3.00. a default gap length weight of 0.10, and weighted end gaps. 

Another example of a useful algorithm is the BLAST algorithm, described in: Altschul ^ aL. J. Mol. 
Biol. 215. 403^10, (1990); Altschul et aL. Nudeic Adds Res. 25:3389-3402 (1997); and Kariin et al„ 
Proc. Natl. Acad. Sd. U.S.A. 90:5873-5787 (1993). A particularty useful BLAST program is the WU- 
BLAST-2 program which was obtained from Altschul el aL. Methods in Enzymology, 266:460-480 
(1996); http://blasl.wustl/edu/blast/ README.htmll. WU-BLAST.2 uses several search param^ers. 
most of which are s^ to the default values. The adjustable parameters are set with the following 
values: overtap span =1 , overtap fradion = 0.125. word threshold (D = H- "Hie HSP S and HSP S2 
parameters are dynamic values and are established by the program itsdf depending upon the 
composition of the particular sequence and composition of the particular database against which the 
sequence of interest is being searched; however, the values may be adjusted to Increase sensitivity. 
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An additional useful algoritlim is gapped BLAST as reported by Altschul et al.. Nud. Acids Res.. 
25:3389-3402. Gapped BLAST uses BLOSUM-62 substitution scores; threshold T parameter set to 9; 
the two-hit m^hod to trigger ungapped extensions: charges gap lengths of k a cost of 10+k; set to 
1 6. and X, set to 40 for database search stage and to 67 for the output stage of the algorithms. 
Gapped alignments are triggered by a score corresponding to -22 bits. 

A % airino add sequence identity value is detenrined by the number of matching identical residues 
divided by the total number of residues of the longer" sequence hi the algned region. The "longer" 
sequence Is the one having the most actual residues in the afigned region (gaps introduced by WU- 
Blast-2 to maximize the aOgnment score are ignored). 

In a similar manner, "percent (%) nudeic add sequence Identity" with resped to the coding sequence 
of the polypeptides identified herdn Is defined as the percentage of nucleotide residues in a candidate 
sequence that are Identical with the nucleotide residues in the coding sequence of the cell cycle 
protein. A preferred method utilizes the BLASTN module of WU-BLAST-2 set to the default 
paramrters, with overtap span and overlap fraction set to 1 and 0.125. respectively. 

The aHgnment may indude the Introduction of gaps in the sequences to be aligned. In addition, for 
sequences whidi contain either more or fewer amino adds than the protdn encoded by the sequence 
of Rgure 1 . It is underetood that In one embodiment, the percentage of sequence identity will be 
determined based on the number of identical amino adds in relation to the total number of amino 
adds. Thus, for example, sequence identity of sequences shorter than that shown in Figure 1 . as 
discussed below, will be deternined using the number of amino adds in the shorter sequence, in one 
entodiment. In percent Identity calculaBons relative wdght Is not assigned to various manifestations 
of sequence variation, such as. Insertions, deletions, substitutions, etc. 

In one embodiment only identities are scored positively (+1) and all forms of sequence variation 
induding gaps are assigned a value of "0". which obviates tiie need for a weighted scale or 
parameters as described below for sequence similarity calculaBons. Percent sequence Identity can be 
calculated, for example, by dividing the number of matdiing Identical residues by the total number of 
residues of the "shortei* sequence In the aligned region and multiplying by 100. The "longer" 
sequence Is the one having ttie most actual residues In the afigned region. 

Thus. I A proteins of the present invention may be shorter or longer than the anino add sequence 
shown In Rgure 1. Thus. In a preferred embodiment, induded witWn the definition of lA proteins are 
portions or fragments of the sequences depided herein. Fragments of lA proteins are considered lA 
proteins if a) they share at least one antigerac eiMtope; b) have at least the Indicated homdogy; c) and 
preferably have lA biological activity as defined her^n. 
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In a preferred embocfiment. as is more fully outlined below, the lA proteins include further anwno acid 
variations, as compared to a wild type insulin, than those outUned herein. In addition, as outOned 
herein, any of the variations depcted herein may be combined in any vay to torn additional novel lA 
proteins. 

In addition. lA proteins can be made that are longer than those depicted in the figures, for example, by 
the addition of epitope or purification tags, as outfined herein, the addition of other fusion sequences, 
etc. For example, the lA proteins of the invention may be fused to other therapeutic proteins such as 
lL-1 1 or to other proteins such as Fc or serum albuirin for pharmacokinetic purposes. See for 
example U.S. Patent No. 5.766.883 and 5.876.969, both of which are expressly incorporated by 
reference. 

In a preferred embodiment, the lA proteins are globally redesigned for improved stabiOty. In one 
aspect of this embodiment, the lA protein comprises variable residues in core residues. 

In one embodiment, the variable core positions are altered to any of the other 19 anino adds. In a 
preferred eirtjodiment, the variable core residues are chosen from Ala. Val. Phe. He. Leu. Tyr , Trp 
and Met. In another preferred embodiment, the variable core residues are chosen from Ala, Val, Leu, 
Me. Phe. Tyr. and Trp. In another preferred embodiment, the variable core residues are chosen from 
Ala. Val. leu, He, and Gly. 

Human Insulin core residues are as follows: positions A2, A3, A16. B11, B15, and B24, whereby "A" 
refers to a residue in the A-chain of insuOn and the number identifies the position within the A-chain. 
Accoidingly. 'B' refere to a residue in the B-chain of insulin and the number identifies the position 
within the B-chain. In the context of the mature insufin, the A-chain comprises residues 1-21 and the 
B-chain residues 22-51 of the amino add sequence shown In Figure IB. In sonw embodiments, when 
referred expDdtly to the B^^ain. residues 22-51 are also referred to B1 - B30, respectively. 
Accordingly, in a preferred eirtwdiment, lA proteins have variable positions selected from positions 
A2. A3, A16, B1 1. B15, and B24. 

In a preferred embodiment, the lA protein of the invention has a sequence that differs from a vflld4ype 
human insufin pftitein in at least one amino add position selected from positions A2, A3, A16. B11. 
B15.orB24. 

Preferred amino add substitutions for each position are as follows: Position A3: He; position A16: Tyn 
position B11 : He; and position B24: Tyr. TTie preferred amino adds in position A2 is lie (wild type) and 
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B15 is Leu (wild type). Any combination of the aijove-listed amino add sulKtitutions is possible and 
should result in an lA protein that is more stable than insuUn. 

In a preferred embotfiment lA proteins have variable positions selected solely from core residues of 
human Insufin. Alternatively, at ieast a majority (51%) of the variable positions are selected from core 
residues, with at least about 75% of the variable positions being preferably selected from core residue 
positions, and at least about 90% of the variable positions being particularly preferred. A specifically 
prefeired embodiment has only core variable positions altered as compared to human insuHn. 

In another aspect of this prefen-ed embodiment, globally redesigned lA proteins comprise variable 
residues in boundary residues. 

Human insulin boundary residues are as follows: positions A5. A15. A17, A19. A21 . B2. B3. 84. 88, 
812, 814 818. 822. 826. and 828. Accordingly, in a preferred embodiment, lA proteins have variable 
positions selected from these positions. 

In a preferred embodiment, the I A protein of the invention has a sequence thai differs from a wild4ype 
human insufin protein in at least one amino add position selerted from positions A5. A15. A17. A19. 
A21. 82. 83. 84, 88. 812, B14. B18. 822. 826. or B28. 

Preferred amino add substitutions for each position are as follows: Position A5: Glu. and Arg: position 
A15- Glu Leu. and Aig: position A17: Lys, Trp. Gin, and Tyr; position A19: Phe; position A21: Asp. 
Gin. arui Arg; position 82: Lys. position 84: Phe and Tyr. position 88: Lys. Leu. Glu; position 812: Aig 
and Lys; position 814: Gki and Trp; position 818: Lys; position 822: Gin; position 826: Phe; and 
posWon KS: Phe and Asn. The preferred amino adds in position 83 is Asn (wild type). Any 
combination of the above-Bsted anino add substitutions is possible and should result In an lA protein 
that is more stable than insufin. 

In another aspect of this Referred embodiment, globally redesigned lA proteins comprise variable 
residues in surfeice residues. 

Human insuOn surface residues are as follows: positions A1, A4. A8. AS. A10. A12. A13. A14. A18. 81, 
85. 86, 89, 810, 813. 816. 817, 821. 825 827. 829. and 830. Accordingly, in a preferred 
entbodiment. lA prolans have variable positions selected from these positions. 

In a preferred embodiment, the lA protein of the invention has a sequence that differs from a wild4ype 
human insulin protein in at least one amino add position selected from positions A1 . A4. A8, A9, AID. 
A12. A13. A14, A18, B1. 85, 86, 89, 810. 813, 816, 817. 821. 825, 827. 829. and 830. 
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po^ible and should resultin an lAproleln thatis more stable than insufin. 

,„ a preferred eml«diment. the lA protein «>n,p,ises a^no add substitutions from core residues, 
and/or boundaiy residues, and/or surface residues. 

,n a orefened aspect of this embodiment, the lA protein comprises the amino acid sequence shown In 
insu«n sequence) and comprises A1-N. AIO^. Aie-Y. Bt^. B25^. and BTT-D. 

,„ another aspect of this embodiment the lApro.eincenplsesthean.no add sequent 
FigureSB. THs sequence sh««s1 4 mufation (27-28% divergence f.mtt,e^^^^^^^^ 
s^uen^) and comprises A14a.A1CH5. A16-Y.A17.Y.A19^. B1-D. B2-K.B4-F. Bll-l. B12-R. B14- 
W. B2S-N. B26-F. and B27-D. 

,„ one preferred aspect of this embodiment, the lA pmtein comprises the artno add sequence shown 
in Rgure 5C. This sequence shows 16 mutation (31-34% divergence 

sequence) and comprisesA1-N.A10-Q.A16-Y.A17-Y.A19^.B1-D.B2-K.B4-F.B8-L.Bl1M.B 
B14-W. B25-N. B26^=. B27-D. and B28-N. 

The following cysteine residues are linked via disulfide bonds with one another in human insuOn: AS- 
A11 A7-B7 andA20-B19. In a preferred embodiment. lA proteins are provided which comprise a 

disulfide bond replacement In one aspect of this embodiment, the following residues were considered 
in the PDA calculation: A2. A3. AS. A6.A7.A11.A15.A16.A19.A20.B2.B7.B11.B15.B18.B19 

B22. and B24. Accordingly. In a prefened embodiment. lA pmteins have variable positions selected 

from these positions. 

in a preferred embodiment the 1 A protein of the Invention has a sequence that differs fhmi a wlld^ 
human InsuBn protein In at least one amino add position selected from positions A2. A3. AS. A6. A7. 
All. A15. A16. A19. A20. 82. B7. B11. BIS. B18. B19. B22. and B24. 
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Preferred anvno add substitutions for each position are as follows: Position A2:Leu; position A3: lie; 
position A5: Arg, Val, and Glu; position A6: Ala; position A7: Ala and Sen position A11 : Ala; position 
A15: Leu. Glu, and Arg; position A16: lie; position A19: Phe; position A20: Asp. Ala. and Set; position 
B2: Asp. Asn. Glu. Gin, Lys. and Arg; position B7: Tyr and Hsp; position B15: Leu; position B19: Ala; 
and position B22: Gin. Any connbination of the atx)ve-Osted anruno add substitutions Is possible and 
should result in an lA protein that is more stable than insulin. Preferred anwno adds at position B1 1 , 
B18, B24 are Leu, Val , and Phe, respectively (all wild type). 

In one aspect of this embodiment, the I A protein comprises the amino add sequence shown in Figure 
3A. This sequence shows 20 mutation (39-40% dh^ergence from the wild type insulin sequence) and 
comprises A1-N. A2-I. A6«A, A7«S. A10-Q. A11-A. A16-I. A17-Y. A19-F. A20-D. B1-D. B4-F, B7-Y. 
B11-I, B12-R. B14-W, B19-A. B25-N. B26-F. and B27-D. 

In another aspect of this embodiment, the following residues were considered for a disulfide bond 
replacement A3. A7. B2. B4. B7, and B11. Accordingly, in a preferred embodiment, lA proteins have 
variable positions selected from these positions. 

In a preferred embodiment the I A protein of the Invention has a sequence ttiat differs from a wild-type 
human insufin protein in at least one attino add position selected from positions A3, A7, B2, B4, B7, 
andBII. 

Preferred anwno add substitutions for each position are as follows: Position A3: lie and Ala; position 
A7: Ala and Ser; position 82: Thr, Asp, Asn, Glu, Gin. Lys. and Arg; position B4: Phe, Tyr, Glu, Lys, 
and Arg; position B7: Ala, Ser, Asp, Asn, Tyr and Hsp; and position B11: He. 

In one aspect of tills embodiment, the lA protein comprises the amino add sequence shown in Rgure 
3B, This sequence shows 4 mutation (8% divergence from the wild type insulin sequence) and 
comprises A7-S. B2-E. B4-Y. and B7-Y. 

In another aspect of this embodiment, the lA protein comprises the amino add sequence shown in 
Figure 3C. TWs sequence shows 2 nuilation (8% divergence from the wild type insuOn sequence) and 
comprises A7-S. and B7-D. 

In a preferred aspect of this embodiment, the lA protein comprises the amino add sequence shown In 
Figure 3D. TWs sequence shows 4 mutation (8% divergence from the wild type insufin sequence) and 
comprises A7-S. B2-T, B4-Y, and B7-Y. 
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In another aspect of this embodiment, the lA protein comprises the amino acid sequence shown in 
Figure 3E. This sequence shows 3 mutation (6% divergence from the wild type insulin sequence) and 
comprises A7-S. B4-Y, and B7-Y. 

In one aspect of this embodiment, the I A protein comprises the amino acid sequence shown in Figure 
3F. This sequence shows 2 mutation (4% divergence from the wild type insuFm sequence) and 
comprises A7-S, B7-E. 

In another preferred aspect of this embodiment, the I A protein comprises the amino acid sequence 
shown in Figure 3G. This sequence shows 2 mutation and 4deletions at positions 81 to 84 (12% 
divergence from the wild type insulin sequence) and comprises A7-E and B7-E. 

Insulin ofigomerizes to dimers, tetramers and hexamers. The hexamers are biologicaWy inactive, 
however, they stabilize the monomeric compounds, paiticulariy in phamiaceutical compositions, 
wherein insuUn fonns a hexamer which is complexed with zinc ions and bound phenolic presen^atives. 
Upon administration into the bloodsfream it is betieved that the (active) monomer dissociates from the 
hexamer complex. In some applications it may be desired to occlude phenoKc preservatives from the 
hexamer complex, but still retain the capat»lity to form hexamers. These lA protein complexes find 
use as slow-acting insulins In contrast to fomiutations which favor the monomeric state. 

Thus, in one prefen-ed embodiment, PDA design is used to generate lA proteins that promote 
hexamer formation, but occlude phenol txndlng. In one aspect of this embodiment, lA proteins are 
generated that are stable and fo^rm hexamers in the absence of phenolic preseivatives. Some of 
these lA proteins may form hexamers that are more stable than the human insulin bound to a phenofic 
compound. In Ms embodiment, the PDB entry Iwav was chosen. For the PDA calculation, the 
entire insulin hexamer complex, consisting of 6 A-chains (chains 1 . 3, 5, 7. 9, and 1 1 in hexamer) and 
6 B-chains (chains 2, 4, 6, 8. 10. and 12 in the hexamer) was used. 

In one embodiment, the following residues were considered in the PDA calculation: Leu-BI 7 (chain 6), 
Va^B2 (chain 8). His^5 (chain 5). Leu-B6 (chain 8). Leu^16 (chain 11). Leu-BII (chain 12). and Ala- 
814 (chain 12). Acooidingly. in a preferred en*odiment. lA proteins have variable positions selected 
from these positions. 

In a preferred embodiment, the lA protein of the invention has a sequence that differs from a wildrtype 
human insulin protein In at least one amino acid position selected from positions B17, B2, 85. B6, 
A16, 811, and 814. 
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Preferred atrtno add substitutions for each position are as follows: Position Leu-B17 (chain 6): Val 
and Glu; . position Val-B2 (chain 8): Phe and Glu; position His-B5 (chain 5): Phe, Trp, and Leu; 
position Leu-A16 (chain 1 1): Phe and Tyn position Leu-B11 (chain 12): lie; and position Ala-B14 (chain 
12):Trp. Phe, and Leu. Ttie preferred amino add at position Leu-B6 (chain 8) is Leu (wild type). Thus, 
prefen-ed substitutions occur at positions B14 and B5. 

In a further embodiment, the B14 and 85 positions were analyzed with respect to possible 
substitutions that promote hexamer formation. 

In a preened embodiment, the lA protein conprises a substitution at airtno acid residue B14. 

In one aspect of this embodiment, the lA protein comprises the amino add sequence shown in Figure 
4C. This sequence shows only 1 mutation (2% divergence from the wild type insulin sequence) and 
comprises B14^. This 1A protein does not bind effidently a phenol preseivative. however, it still forms 
a hocamer. 

In another aspect of this embodiment, the lA protein comprises the amino add sequence shown In 
Figure 4D. This sequence shows onV 1 mutation (2% divergence from the wild type insuTin sequence) 
and comprises B14AW. This lA protein does not bind effidently a phenol presenrative. however, it still 
forms a hexamer. 

In another aspect of this embodiment, the lA protein comprises the amino add sequence shown in 
Rgure 4F. This sequence shows only 1 mutation (2% divergence from the wild type insufin sequence) 
and comprises B14-Y. This lA protein does not bind effidently a phenol presewative, however. It still 
forms a hexamer. 

In another aspect of this embodiment, the lA protein comprises the amino add sequence shown in 
Figure 4G. This sequence shows only 1 mutation (2% divergence from the wild type insulin sequence) 
and comprises B14^. This lA protein does not bind effidently a phenol preseivative. however, it still 
forms a hexamer. 

In a another prefened embodiment, the lA protein comprises a substitution at amino acid residue B5. 

In one aspect of tWs embodiment the lA protein comprises the amino add sequence shown in Rgure 
48. This sequence shows only 1 mutation (2% divergence from the wild type InsuKn sequence) and 
comprises 85-F. Wis lA protein does not bind effidently a phenol preseivative. however. It still forms 
a hexamer. 
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in another preferred embodiment, the lA protein comprises substitutions at amino acid residuesBS and 
B14. 

In one aspect of this embodiment, the lA protein comprises the amino add sequence shown in Rgure 
4A. This sequence shows 2 mutations (42% divergence from the wild type insulin sequence) and 
comprises B5-F and B14-F. This I A protein does not bind efficiently a phenol preseivative, however, it 
Still forms a hexamer. 

In another aspect of this embodinrent, the lA protein comprises the amino add sequence shown In 
Rgure 4E. This sequence shows 2 mutations (42% divergence from the wild type insulin sequence) 
and comprises B5-F and B14-W. This lA protein does not bind effidently a phenol presewative. 
however, it still forms a hexamer. 

In a preferred embodiment, the I A proteins of the invention are human insuBn conformers. By 
"conformed herein is meant a protein that has a protein backbone 3D structure that is virtually the 
same but has significant differences in the amino add side chains. That is, the lA proteins of the 
invention define a confonner set. wherein all of the proteins of the set share a backbone structure and 
yet have sequences that differ by at least 3-5%. The three dimensional backbone structure of an lA 
protein thus substantially con-esponds to the three dimensional badcbone structure of human insulin 
•Backbone" In this context means the non^lde chain atoms: the nitrogen, carbonyl carijon and 
oxygen, and the a-cirbon. and the hydrogens attached to the nitrogen and a^arbon. To be 
considered a confonner. a protein must have backbone atoms that are no more than 2 A from the 
human insuFin structure, with no more than 1 .5 A being preferred, and no more than 1 A being 
particulariy preferred. In general, these distances may be determined in two ways. In one 
embodiment, eadi potential confonner Is ccystalUzed and its three dimensional structure determined. 
Alternatively, as the former is quite tedious, the sequence of each potential conformer is run in the 
PDA program to deterrwne whether It is a confonner by examining the energy of the sequence. 

lA proteins may also be identified as being encoded by lA nudeic adds. In the case of the nudeic 
add. the overall homology of the nucldc add sequence Is commensurate with amino add homology 
but takes into account the degeneracy in the genetic code and codon Was of different organisms. 
Accoidingly. the nudeic add sequence homology may be either lower or Wgher than that of the protein 
sequence, with lower homology being prefe^ed. 

In a preferred embodiment, an lA nudeic add encodes an lA protein. As will be appredated by those 
in the art, due to the degeneracy of the genetic code, an extremely large number of nudeic adds may 
be made, an of which encode the lA proteins of the present invention. Thus, having identified a 
particular an^no add sequence, those skilled In tfie art couW make any number of different nucleic 
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adds, by simply modilying the sequence of one or more codons in a way which does not change the 
airino add sequence of the lA protein. 

In one embodiment, the nudeic add homology is deteitrined through hybridization studies. Thus, for 
example, a nucleic add. which hybridize under high stringency to a nudeic add sequence encoding 
insulin (afl. see GenBank accession numbers J00265 M10039) and encode an lA protein is 
considered an lA gene. 

High stringency conditions are known in the act; see for example Sambrook et al.. Molecular Cloning: 
A Laboratoiy Manual. 2d Edition. 1989. and Short Protocols in Molecular Biology, ed. Ausubel. et al.. 
both of whidi are hereby incorporated by reference. Stringent condlUons are sequence^ependent 
and will be different in different drcumslances. Longer sequences hybridize spedfically at higher 
temperatures. An extensive guide to the hybridization of nudeic adds is found in Hjssen. Techniques 
in Biodiemistcy and Molecular Biology-Hybridization with Nudeic Add Probes. "Ovemew of prindples 
of hybridization and the strategy of nudeic add assays" (1993). Generally, stringent condMons are 
seleded to be about 5-1 0'C.lower than the thermal melting point (J J for the spedfic sequence at a 
defined ionic strength and pH. The T„ is the temperatuiJe (under defined Ionic strength. pH and 
nucleic add concentration) at whidi 50% of the probes complementary to the target hybridize to the 
target sequence at equiObrium (as the target sequences are present in excess, at T„. 50% of the 
probes are occupied at equilibrium). Stringent conditions will be those in which the salt concentration 
is less than about 1 .0 M sodium ion. typically about 0.01 to 1 .0 M sodium Ion concentration (or other 
salts) at pH 7.0 to 8.3 and the temperature Is at least about 30*C for short probes (e,g. 10 to SO 
nucleotides) and at least about 60*C for long probes (e.g. greater than 50 nudeotides). Stringent 
conditions may also be achieved with the adcfition of destabilizing agents such as forman«de. 

In another embodiment, less stringent hybridization conditions are used; for example, moderate or low 
stringency conditions may be used, as are known In the art; see Sambrook et al., supra; Ausubel et 
al., supra, and Hjssen, supra. 

The lA proteins and nudeic adds of the present inventiori are recombinant As used herein, "nudeic 
add" may refer to either DNA or RNA, or molecules which contain both deoxy- and ribonudeotides. 
The nucleic adds include genonic DNA. cDNA and oDgonucieotides including sense and anti-sense 
nucleic adds. Such nucleic adds may also contain modifications in the ribose-phosphate backbone to 
bicrease stabifity and half life of such molecules in physiotogical environments. 

The nudeic add nray be double stranded, single stranded, or contain portions of botii double stranded 
or single stranded sequence. As will be appredated by those In the art, the def«ction of a single 
strand CWatson") also defines the sequence of the other strand fCrick"). By the temi "recombinant 
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nucleic add" herein is meant nucleic acid, originally fonned in vitro, in general, by the manipulation of 
nucleic add by endonucleases, in a fonn not nonnally found in nature. Thus an isolated lA nudeic 
add, in a linear form, or an expression vector fomned in vitro by ligating DMA molecules that are not 
nonnally joined, are both considered recomttnant for the purposes of this invention. It is understood 
that once a recombinant nudeic add is made and reintroduced into a host cell or organism, it will 
replicate non-recombinantly. i.e. using the in vivo cellular machinery of the host cell rather than in vitro 
manipulations; however, such nudeic adds, once produced recombinantly, although subsequently 
repBcated non-recombinantly. are still considered recombinant for the purposes of the invention. 

Similariy. a "recombinant protein" is a protein made using recombinant techniques, i.e. through the 
expression of a recombinant nudeic add as depicted above. A recomlwnant protein is distinguished 
from naturally occumng protein by at least one or more characteristics. For example, the protein may 
be isolated or purified away from some or all of the proteins and compounds with which it is nonnally 
assodated in its wild type host, and thus may be substantially pure. For example, an isolated protein 
is unaccompanied by at least some of the material with which it is normally assodated in its natural 
state, preferably constituting at least about 0.5%. more preferably at least about 5% by weight of the 
total protdn in a given sample. A substantially pure protein comprises at least about 75% by weight of 
the total protein, with at least about 80% being prefen-ed, and at least about 90% being particularly 
preferred. The definition Includes the production of an lA protein from one organism in a different 
organism or host cell. Alternatively, the protein may be made at a significantly higher concentration 
than is nonnally seen, through the use of an indudble promoter or high expression promoter, such that 
the protein is made at increased concentration levels. Furthermore, all of the I A proteins outlined 
herein are in a form not nomialV found in nature, as they contain amino add substitutions, insertions 
and deletions, with substitutions being prefen-ed, as discussed below. 

Also induded \Mthin the definition of lA proteins of the present invention are anino add sequence 
variants of the I A protein sequences outlined herein and shown in the Figures. That is, the lA proteins 
may contain additional variable positions as compared to human insulin. These variants fall into one 
or more of three dasses: substitutional, insertional or deletional variants. These variants ordinarily are 
prepared by site spedfic mutagenesis of nudeotides In the DMA encoding an lA protein, using 
cassette or PGR mutagenesis or other techniques well known in the art, to produce DMA encoding the 
variant, and thereafter expressing the DNA in recomtMnant cell culture as outlined above. However, 
variant lA protein fragments having up to about 100-150 residues may be prepared by in vitro 
synthesis using established techniques. Spedfically. lA proteins comprising only an^no add 
sequences for the A^aln may be prepared by in vitro synttiesls. Simllarty. lA proteins comprising 
only amino add sequences for the B-chain may be prepared by in vitro synthesis. The individually 
prepared lA protein A-chain and lA protein B-chain can be combined after synthesis to form the 
mature lA protein. Techniques for the synttiesis of insufin A-chains and B-chains and subsequent 



43 



wo 00/69901 



PCT/USOO/13764 



formation of mature insulin, i.e.. the insulin comprising an A-chain and a B-chain appropriately linked 
via the disulfide bonds descrit)ed herein, are known in the art The same techniques apply to the 
formation of mature lA proteins starting out with individual lA A-chains and lA B-chains. Amino acid 
sequence variants are characterized by the predetennned nature of the variation, a feature that sets 
them apart from naturally occurring allelic or Interspecies variation of the lA protein amino add 
sequence. The variants typically exhibit the same quafitative biological activity as the naturally 
occurring analogue, although variants can also be selected which have modified characteristics as wiU 
be more fully outOned betow. 

While the site or region for Introducing an anrano add sequence variation is predetermined, the 
mutation per se need not be predetennned. For example, in order to oplinwze the performance of a 
mutation at a given site, random mutagenesis may be conducted at the target codon or region and the 
expressed I A variants saeened for the optimal combination of desired activity. Techniques for making 
substitution mutations at predeternined sites in DNA having a known sequence are well known, for 
example. Ml 3 primer mutagenesis and PGR mutagenesis. Screening of the mutants is done using 
assays of lA protein activities. 

Amino add substitutions are typically of single residues; insertions usually will be on the order of from 
about 1 to 20 amino adds, although considerably larger insertions may be tolerated. Deletions range 
from about 1 to about 20 residues, although in some cases deletions may be much larger. 

Substitutions, deletions, insertions or any combination thereof may be used to arrive at a final 
derivative. GeneralV these changes are done on a few amino adds to minimize the alteration of the 
molecule. However, larger dianges may be tolerated in certain drcumstances. When small 
alterations in the characteristics of the lA protein are desired, substitutions are generally made in 
accordance witii the following chart: 

Chart I 

Qriolnai Residue Fxemplary Substitutions 



Ala 

is G^His 



Ser 

Lys 



Asp 



Glu 



Cys Ser, Ala 



Gin 



Asn 



Glu Asp 
Gly 

His Asn. Gin 

lie Leu.Val 

Leu "«'Val 

Lys Arg. Gin. Glu 

Met Irf"'"® ^ 

Phe Met,Leu.Tyr 



Ser 



Thr 
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Thr 
Trp 
Tyr 
Val 



Ser 
Tyr 

Trp, Phe 
lie, Leu. 



Substantial changes in function or immunological identity are made by selecting substitutions that are 
less consewative than those shown In Chart I. For example; substitutions may be nriade which wore 
significantly affect the structure of the polypeptide backbone In the area of the alteration, for example 
the alpha-heBcal or beta-sheet stmcture; the charge or hydrophobldty of the molecule at the taiget 
site; or the bulk of the side chain. The substitutions which in general are expected to produce the 
greatest changes in the polypeptide's properties are those in which (a) a hydrophific residue, e.g. seiyl 
or threonyl, is substituted for (or by) a hydrophobic residue, e.g. leucyl. isoleucyl. phenylalanyl, valyl or 
alanyl; (b) a cysteine or proBne is substituted for (or l>y) any other residue; (c) a residue having an 
electropositive side chain, e.g. lysyl, arginyl. or histidyl, is substituted for (or by) an electronegative 
residue, e.g. glutamyl or aspartyl; or (d) a residue having a bulky side chain, e.g. phenylalanine, is 
substituted for (or by) one not having a side chain, e.g. glycine. 

The variants typically exhibit the same qualitative biological activity and will elicit the same Immune 
response as the original lA protein, although variarits also are selected to modify the characteristics of 
the lA proteins as needed. Alternatively, ttie variant may be designed such that the bological activity 
of the lA protein Is altered. For example, glycosylation sites may be altered or removed. Sintilarty. the 
biological function may be altered; for exannple. in some instances it may be desirable to have more or 
less potent insulin activity. 

The lA proteins and nucleic adds of the invention can be made In a number of ways. Individual 
nucleic adds and proteins can be made as known in the art and outfined below. Alternatively, libraries 
of lA proteins can be made for testing. 

In a preferred embodiment, sets or fibraries of lA proteins are generated from a probability distribution 
table. As outiined herein, there are a variety of metiiods of generating a probability distribution table, 
induding using PDA. sequence aUgnments, fbrcefieW calculations such as SCMF calculations, etc. In 
addition, the probability distribution can be used to generate information entropy scores for each 
position, as a measure of the mutational frequency observed in the library. 

In ttiis embodiment, the frequency of eadi amino add residue at each variable position In the fist Is 
kJentified. Frequendes can be thresholded. wheran any variant frequency lower than a cutoff is set to 
zero. This cutoff Is preferably 1%. 2%. 5%. 10% or 20%, with 10% being particularty preferred. These 
frequendes are then bialt into the lA protein Bbrary. That is, as above, these variable positions are 
colleded and all possible combinations are generated, but the amino add residues that "fiir the library 
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are utiDzed on a frequency basis. Thus, in a non-frequency based library, a variable position that has 
5 possible residues will have 20% of the proteins comprising that variable position with the first 
possible residue, 20% with the second, etc. However, in a frequency based library, a variable position 
that has 5 possible residues with frequencies of 10%, 15%. 25%. 30% and 20%, respectively, will have 
1 0% of the proteins comprising that variable position with the first possible residue, 1 5% of the 
proteins with the second residue, 25% with the tNrd. etc. As will be appreciated tyy those in the art, the 
actual frequency may depend on the method used to actually generate the proteins; for example, 
exact frequencies may be possible when the proteins are synthesized. However, when the frequency- 
based primer system outlined below is used, the actual frequencies at each position will vaiy, as 
outlined below. 

As will be appreciated by those in the art and outlined herein, probability distribution tables can be 
generated In a variety of ways. In addition to the methods outlined herein, self^nsistent mean field 
(SCMF) methods can be used In the direct generation of probability tables. SCMF is a detenrinistic 
computational method that uses a mean field description of rotamer interactions to calculate energies. 
A probability table generated in this way can be used to create libraries as described herein. SCMF 
can be used in three ways: the fi-equencies of amino adds and rotamers for each amino acid are listed 
at each position; the probabilities are determined directly from SCMF (see Delarue et la. Pac. Symp. 
Biocomput. 109-21 (1997). expressly incorporated by reference). In addition, highly variable positions 
and non-variable positions can be identified. Alternatively, another method is used to determine what 
sequence is jumped to during a search of sequence space; SCMF is used to obtain an accurate 
energy for that sequence; tWs energy is then used to rank it and create a rank-ordered list of 
sequences (sitwiar to a Monte Cario sequence list). A probabifity table showing the frequencies of 
arnno adds at each position can then be calculated from this Rst (Koehl et al. J. Mol. Biol. 239249 
(1994); Koehl et aL. Nat. Struc. Biol. 2:163 (1995); Koehl et al.. Curr. Opin. Struct. Biol. 6:222 (1996); 
Koehl et al.. J. Mol. Bio. 293:1 183 (1999); Koehl et al., J. Mol. Biol. 293:1161 (1999); Lee J. MoL BioL 
236:918 (1994); and Vasquez Biopolymers 36:53-70 (1995); all of which are expressly incorporated by 
reference. Similar methods Include, but are not fimited to, OPLS-AA (Jorgensen, et al. J. Am. Chem. 
Soc. (1996). V 1 1 8. pp 1 1225-1 1236; Jorgensen. W.L; BOSS, Version 4.1 ; Yale University: New 
Haven. CT (1999)); OPLS (Jorgensen, et al., J. Am. Chem. Soc. (1988). v 110, pp 1657ff; Jorgensen, 
et al.. J Am. Chem. Soc. (1990). v 112. pp 4768ff); UNRES (United Reskiue Forcefield; Liwo, et al.. 
Protein Science (1993). v 2, pp1697-1714; Uwo, et al., Protein Sdence (1993), v 2. pp1715-1731; 
Uwo, et al.. J, Comp. Chem. (1997), v 18, pp849.873; Liwo, et al. J. Comp. Chem. (1997). v 18, 
pp874^84; Uwo, et aL, J. Comp. Chem. (1998), v 19, pp259-276; Forcefield for Protein Structure 
Prediction (Uwo. et al. Proc. Natl. Acad. Sd. USA (1999). v 96. pp5482-5485); ECEPP/3 (Uwo et aL. 
J Protein Chem 1994 May;13(4):375-80); AMBER 1 .1 force field (Weiner, et aL. J. Am. Chem. Soc. 
v106. pp765-784); AMBER 3.0 force field (U.C. Singh et aL. Proc. Natl. Acad. Sd. USA. 82:755-759); 
CHARMM and CHARMM22 (Brooks, et al., J. Comp. Chem v4, pp 187-217); cvff3.0 
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(Dauber-Osguthorpe, et al..(1988) Proteins: Structure. Function and Genetics. v4.pp31-47): cfOI 
(Maple, et al.. J. Comp. Chem. v15. 162-182); also, the DISCOVER (cvff and cffM) and AMBER 
forcefieWs are used in the INSIGHT molecular modeling package posymrt^flSI. San Diego Cafifomia) 
and HARMM is used In the QUANTA molecular modeling package (BlosymrtVISI, San Diego 
CaBfomia). 

In addition, as outUned herein, a preferred method of gerierating a probaWnty distribution table Is 
through the use of sequence alignment programs. In addition, the probability table can be obtained by 
a combination of sequence alignments and computational approaches. For example, one can add 
amino adds found in the alignment of homologous sequences to the result of the computation. 
Preferable one can add the wiW type amino acid identity to the probabifity table if it is not found in the 
computation. 

As win be appreciated, an lA protein ibrary created by recombining variable positions and/or residues 
atthe variable position may not be in a rank-ordered fist. In some embodiments, the entire fist may 
just be made and tested. Altematively. in a preferred embodiment, the lA protein Bbrary is also in the 
form of a rank onlered fist. This may be done for several reasons, including the size of the Bbrary is 
still too big to generate experimentally, or for predicHve purposes. This may be done in several ways. 
In one embodiment, the Bbraiy is ranked using the scoring ftindions of PDA to rank the Bbrary 
members. Alternatively. statlsUcal methods couW be used. For example, the Bbrary may be ranked by 
frequency score; that is. proteins containing the most of high frequency residues could be ranked 
higher, etc. This may be done by adding or multiplying the frequency at each variable position to 
generate a numerical score. Simlarty. the Bbrary different positions couW be weighted and then the 
proteins scored; for exanqile. those containing certan residues couU be artntrariiy ranked. 

In a prefened embodiment, the different protein members of the lA protein Obraiy may be chemically 
synthesized. This is particularly useful when the designed protans are short, preferably less than 150 
amino acids in length, with less than 100 amino adds being preferred, and less than 50 an^no adds 
being particularly preferred, although as is known in the art, longer proteins can be made chemical|y w 
enzymatically. See for example Wilken et al. Cun^. Opin. Biotechnol. 9:412-26 (1998). hereby 
expressly incorporated by reference. 

In a prefen-ed embodiment, particularly for longer proteins or proteins for which large samples are 
desired, the Ubiary sequences are used to create nudeic adds such as DNA wtuch encode the 
menlber sequences and which can then be doned into host c^. expressed and assayed, if desired. 
Thus, nudeic adds, and parlicularty DNA, can be made which encodes each member protein 
sequence. This is done using well known procedures. The dioice of codons. suitable expression 
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vectors and suitable host cells will vary dependjng on a number of factors, and can be easily optinvzed 
as needed. 

In a preferred embodiment, multiple PGR reactions with pooled oligonucleotides is done, as is 
generally depicted in Figure 6. In this emtxxiiment overlapping oligonucleotides are synthesized 
which correspond to the full length gene. Again, these oligonucleotides may represent all of the 
different amino adds at each variant position or subsets. 

In a preferred embodiment, these oligonucleotides are pooled in equal proportions and multiple PGR 
reactions are performed to create full length sequences containing the combinations of mutations 
defined by the library. In addition, this may be done using error-prone PGR methods. 

In a preferred embodiment the different oligonucleotides are added In relative amounts con-esponding 
to the probability distribution table. The multiple PGR reactions thus result In full length sequences 
with the desired combinations of mutations in the desired proportions. 

The total nun^er of oligonucleotides needed is a function of the number of positions being mutated 
and the number of mutations being considered at these positions: 

(number of oUgos for constant positions) + M1 + M2 + M3 +... Mn = Cotal number of ofigos required), 
where Mn Is the number of mutations considered at position n In the sequence. 

In a preferred embodiment, each overlapping oligonucleotide comprises only one position to be varied; 
in altemate embodiments, the variant positions are too close together to allow tWs and multiple 
variants per oligonucleotide are used to allow complete recomtMnation of all the possibilities. That is, 
each oligo can contain the codon for a single position being nurtated, or for more than one position 
being mutated. The multiple positions being mutated must be close in sequence to prevent the oligo 
lengtti from being impractical. For multiple mutating positions oh an oligonucleotide, particular 
combinations of mutations can be Included or excluded In the Gbrary by including or excluding the 
o6gonucleotide encoding that combination. For example, as discussed herein, there may be 
correlations between variable regions; ttiat is, when position X is a certain residue, position Y must (or 
must not) be a particular residue. These sets of variable positions are sometimes refen-ed to herein as 
a "cluster". When the clusters are comprised of residues close togetiier. and thus can reside on one 
oligonucleotide primer, tiie clusters can be set to Uie "good" correlations, and efininate the bad 
combinations tfiat may decrease tiie effectiveness of ttie library. However, If the residues of the cluster 
are far apart In sequence, and thus win reside on cfifferent oligonucleotides for synthesis, it may be 
desirable to either set the residues to the "good" correlation, or eliminate them as variable residues 
entirely. In an altemative embodiment, the library may be generated in several steps, so ttiat the 
cluster mutations only appear together. This procedure, i.e. the procedure of Identifying mutation 
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dusters and either placing them on the same otigonucleotides or eliminating them from the library or 
Gbraiy generation In several steps presen/ing clusters, can considerably enrich the experimental libraiy 
with properly folded protein. Identification of clusters can be carried out by a number of ways, e.g. by 
using known pattern recognition methods, comparisons of frequencies of occurence of mutations or by 
using energy analysis of the sequences to be experimentally generated (for example, if the energy of 
interaction Is high, the positions are conrelated). These con'elations may be positional con'elations (e.g. 
variable positions 1 and 2 always change together or never change together) or sequence correlations 
(e.g. if there is residue A at position 1. there is always residue B at po^on 2). See: Pattern discovery 
in Biomolecular Data: Tools. Techniques, and Applications; edited by Jason T.L. Wang, Bruce A. 
Shapiro, Dennis Shasha. New York: Oxford University, 1999; Andrews. Harry 0. Introduction to 
mathematical techniques in pattern recognition; New Yoric. Wiley-lnteiscience 11972]; Applications of 
Pattern Recognition; Editor, K.S. Fu. Boca Raton, Fla. ORG Press, 1982; Genetic Algorithnre for 
Pattern Recognition; edited by Sankar K. Pal. Paul P. Wang. Boca Raton: ORG Press. c1996; Pandya. 
Abhijit S., Pattern recognition with neural networics In G++ / Abhifit S. Pandya, Robert B. Macy. Boca 
Raton, Fla.: CRC Press, 1996; Handbook of pattern recognition & computer vision / edited by C.H. 
Chen, LF. Pau. P.S.P. Wang. 2nd ed. Singapore; River Edge, N.J.: Worid Scientific, c1999; 
Friedman. Introduction to Pattem Recognition: Statistical. Stwctural, Neural, and Fuzy Ugic 
Approaches; River Edge, N.J.: World Scientific, c1999, Series title: Series in machine perception and 
artificial inteingence; vol. 32; all of which are expressly Incorporated by reference. In addition, 
programs used to search for consensus motife can be used as well. 

In addition, correlations and shuffling can be fixed or optimized by altering the design of the 
ongonucleotides; that is. by deciding where the oligonucleotides (primers) start and stop (e.g. where 
the sequences are "cut"). The start and stop sites of ofigos can be set to maximize the number of 
dusteis that appear in single oligonucleotides, thereby enriching ttie Obrary wtth higher scoring 
sequences. Different oligonucleotide start and stop site options can be computationally nrwdeled and 
ranked accoiding to number of clusters tiiat are represented on single oligos, or ttie percentage of tite 
resulting sequences consistent with the predicted library of sequences. 

The total number of oDgonucleotides required increases when multiple mutable positions are encoded 
by a single oligonucleotide. The annealed regions are the ones that remain constant, Le. have the 
sequence of the reference sequence. 

OBgonudeotides with insertions or deletions of codons can be used to create a library expressing 
different lengtti proteins. In particular computational sequence screening for insertions or deletions can 
result in secondaiy libraries defining different lehgtti proteins, which can be expressed by a libraiy of 
pooled oligonudeotide of different lengths. 
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In a preferred embodiment the lA protein Gbrary is done by shufffing the family (e.g. a set of variants); 
that is, some set of the top sequences Of a rank-ordered list is used) can be shuffled, either with or 
without error-prone PGR. "Shuffling" in this context means a recombination of related sequences, 
generally in a random way. It can indude "shuffling" as detined and exempDfied In U.S. Patent Nos, 
5.830,721; 5.811.238; 5.605.793; 5,837.458 and PCT US/19256, all of which are expressly 
incorporated by reference In thar entirety. This set of sequences can also be an artificial set; for 
example, from a probabifity table (for example generated using SCMF) or a Monte Cario set Sin^larly. 
the "family" can be the top 10 and the bottom 10 sequences, the top 100 sequence, etc. This may also 
be done using error-prone PGR. 

Thus, in a preferred embodiment, in silico shuffling is done using the computationaf methods 
described herein. That is, starting with either two libraries or two sequences, random recombinations 
of the sequences can be generated and evaluated. 

In a preferred embodiment en^or-prone PGR is done to generate the lA protein fibrary. See U.S. 
Patent Nos. 5.605.793. 5.811.238. and 5.830.721. all of which are hereby incorporated by reference. 
This can be done on the optimal sequence or on top members of the Bbrary, or some other artificial set 
or family. In this embodiment, the gene for the optimal sequence found In the computational screen of 
the primaiy fibrary can be synthesized. Error prone PGR is then performed on the optimal sequence 
gene in the presence of oligonucleotides that code for the mutations at the variant positions of the 
Bbrary (l^as oligonucleotides). The addition of the oligonucleotides will create a Uas favoring the 
incorporation of the mutations in the fibrary. Alternatively, only oBgonucleotides for certain mutations 
may be used to bias the fibrary. 

In a preferred embodiment gene shuffling with error prone PGR can be perfomied on the gene for the 
optimal sequence, in the presence of bias oligonucleotides, to aeate a DMA sequence fibrary that 
reflects the proportion of the mutations found in the lA protein fibrary. The choice of the bias 
oBgonucleotides can be done in a variety of ways; they can be chosen on the basis of their frequency, 
i.e. ofigonudeotides encoding high mutational frequency positions can be used; alternatively, 
oligonucleotides containing the most variable positions can be used, such that the diversity is 
increased; if the secondary Bbrary is ranked, some nunUser of top scoring positions can be used to 
generate Was ofigonudeotides; random positions may be chosen; a few top scoring and a few low 
scoring ones may be chosen; etc. What is important is to generate new sequences based on preferred 
variable positions and sequences. 

In a preferred embodiment PGR using a wfid type gene or other gene can be used, as is 
schematically depided in Rgure 7. in this embodiment a starting gene is used; generally, although 
this is not required, the gene is usually the wild type gene. In some cases it may be the gene encoding 
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the global optintized sequence, or any other sequence of the fist, or a consensus sequence obtained 
e.g. from aligning homologous sequences from different organisms. In this embodiment, 
oligonucleotides are used that correspond to the variant positions and contain the different amino 
adds of the fibrary. PGR is done using PGR prlnrrers at the tenroni. as is known In the ait. Wis 
provides two benefits; the first is that this generally requires fewer ofigonudeotides and can result In 
fewer errors. In addition, it has experimental advantages in that if the vwld type gene is used, it need 
not be synthesized. 

In addition, there are several other techniques that can be used, as exempfified in the figures, e.g. 
Figures 8-10. In a preferred embodiment, Ogation of PGR products is done. 

In a prefen-ed embodiment, a variety of additional steps may be done to the lA protein library; for 
example, further computational processing can occur, different I A protein fibraries can be recombined. 
or cutoffs from different Ubraries can be combined. In a preferred embodiment, an lA protein library 
may be computationally remanipulated to form an addtional lA protein library (sometimes referred to 
herein as "tertiary fibraries"). For example, any of the lA protein fibrary sequences may be chosen for a 
second round of PDA, by freezing or fixing some or all of the changed positions in the first fibrary. 
Alternatively, only changes seen in the last probability distribution table are allowed. Alternatively, the 
stringency of tiie probability table may be altered, either by increasing or decreasing ttie cutoff for 
Inclusion. Sin^larly, the I A protein fibrary may be recombined experimentally after the first round; for 
example, the best gene/genes from ttie first screen may be taken and gene assembly redone (using 
techniques outfined below, multiple PGR. error prone PGR. shuffling, etc.). Alternatively, ttie fi^gments 
from one or more good gene(s) to change probabifities at some positions. This Uases the search to an 
area of sequence space found In the first round of computational and experimental screening. 

In a preferred embodiment, a tertiary fibrary can be generated firom combining different lA protein 
fibraries. For example, a probabifity distribution table from a first lA protein fibrary can be generated 
and recombined. either computationally or experimentally, as outilned herein. A PDA lA protein fibrary 
may be combined with a sequence afignntent lA protein fibrary. and either recondwned (again, 
computationally or experimentally) or just the cutoffs from each joined to make a new tertiary fibrary. 
The top sequences from several fibraries can be recombined. Sequences from the top of a fibrary can 
be combined witfi sequences from the bottom of ttie fibrary to more broadiy sample sequence space, 
or only sequences distant from the top of the fibrary can be comtnned. lA protein fibraries that analyzed 
different parts of a protein can be corr^ned to a teitiaiy fibraiy ttiat treats the combined parts of the 
protein. 

In a preferred emt)odiment, a tertiary fibrary can be generated using correlations in an lA protein 
fibrary. That is. a residue at a first variable position may be correlated to a residue at second variable 
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position (or correlated to residues at additional positions as well). For exaniple. two variablB positions 
may sterically or electrostatically interact, such that if the first residue is X. the second residue must be 
Y. This may be either a positive or negative con-elation. 

Using the nucleic adds of the present Invention v4iich encode an lA protein, a variety of expression 
vectors are made. The expression vectors may be either self-repficating extrachromosomal vectors or 
vectorcvrfiich integrate Into a host genome. Generally, these expression vectois include 
transcriptional and translational regulatory nucleic add operably Onked to the nuddc add encoding the 
I A protein. The term "control sequences" refers to DNA sequences necessary for the expression of an 
operably linked coding sequence in a particular host organism. The control sequences that are 
suitable for prokaryotes. for example, indude a promoter, optionally an operator sequence, and a 
ribosome binding site. Eukaryotic cells are known to utifize promoters, polyadenylatioh signals, and 
enhancers. 

Nudeic add is "operably linked" when it is placed into a fundional relationship witti anottier nucleic 
add sequence. For example. DNA for a presequence or secretory leader is operably finked to DMA 
for a polypeptide if it is expressed as a preprotein that partidpates In the secretion of the polypeptide; 
a promoter or enhancer is operably linked to a coding sequence if it affects the transcription of the 
sequence; or a ribosome binding site Is operably finked to a coding sequence if it is positioned so as to 
fadOtate translation. 

In a preferred embodiment, when the endogenous secretory sequence leads to a low level of secretion 
of ttie naturally occurring protein or of ttie lA protein, a replacement of tiie naturally occuning secretory 
leader sequence Is desired. In this embodiment, an unrelated secretory leader sequence Is operably 
linked to an lA protein encoding nucleic add leading to increased protein secretion. Thus, any 
secretory leader sequence resulting in enhanced secretion of the lA protein, when compared to the 
secretion of insufin and its secretory sequence, is desired. Suitable secretory leader sequences that 
lead to the secr^on of a protein are know In the ait. 

In anottier prefeaed embodiment, a secretory leader sequence of a naturally occurring protein or a 
protein Is removed by technkjues known in the art and subsequent expresston results In IntraceUular 
accumulation of the recombinant protein. 

Generally, "operably linked" means that the DNA sequences being linked are contiguous, and. In the 
case of a secretory leader, contiguous and in reading phase. However, enhancers do not have to be 
contiguous. Unking is accompOshed by figation at convenient restriction sites. Ifsuch sites do not 
exist, the synttietic ofigonudeotide adafrtors or linkers are used In accordance with conventional 
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practice. The transcriptional and translational regulatory nucleic add will generally be appropriate to 
the host cell used to express the fusion protein; for example, transcriptional and translational 
regulatory nucleic add sequences from Bacillus are preferably used to express the fusion protein in 
Bacillus. Numerous types of appropriate expression vectors, and suitable regulatory sequences are 
known in the art for a variety of host cells. 

In general, the transcriptional and translational regulatory sequences may include, but are not Gnvted 
to, promoter sequences, ribosomal binding sites, transcriptional start and stop sequences, 
translational start and stop sequences, and enhancer or activator sequences. In a prefen-ed 
embodiment, the regulatory sequences include a promoter and transcriptional start and stop 
sequences. 

Promoter sequences encode either constitutive or indudble promoters. The promoters may be either 
naturally occurring promoters or hybrid promoters. Hybrid promoters, which combine elements of 
more than one promoter, are also known in the art, and are useful in the present invention. In a 
prefered embodiment, the promoters are strong promoters, allowing high expression in cells, 
particulariy mammalian cells, such as the CMV promoter, particularty In combination with a Tet 
regulatory element 

In addition, the expression vector may comprise additional elements. For exairple, the expression 
vector may have two replication systems, thus allowing it to be maintained in two organisms, for 
example in mammalian or insect cells for expression and in a prokaryotic host for cloning and 
ampBllcation. Furthermore, for Integrating expression vectors, the expression vector contains at least 
one sequence homologous to the host cell genome, and preferably two homologous sequences which 
flank the expression construct The integrating vector may be directed to a spedfic locus in the host 
cell by selecting the appropriate homologous sequence for Indusion in the vector. Constructs for 
integrating vectors are well known in the art 

In addition, in a preferred embodiment, the expression vector contiains a selectable nnaricer gene to 
allow the selection of transtbrmed host cells. Selection genes are well known in the art and will vary 
with the host cell used. 

A preferred expression vector system is a retroviral vector system such as is generally described In 
PCTAJS97/01019 and PCTAJS97/01048, both of which are hereby expressly incorporated by 
reference. 
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In a preferred embodiment, the expression vector comprises the components described above and a 
gene encoding an lA protein. In this aspect, only one species of an lA protein will be expressed In the 
cell comprising the expression vector. 

In one preferred embodiment, it is desired to express two different lA proteins (variant A and variant B) 
and thus, two expression vectors, one comprising a gene coding lA protein variant A, the other one 
comprising a gene coding for lA protein variant B are introduced into the same host cell. This allows 
formation of a preferred lA protein dimer. 

\n one aspect of this embodiment, an expression vector is constmcted that comprises two lA genes 
encoding two different lA proteins (variant A and variant B). In one aspect of this embodiment, a 
polydstronic gene can be constructed as is known in the art for co-expression in a host cell. 

As will be appreciated by those in the art, ail combinations are possible and accordingly, as used 
herein, the combination of components, comprised by one or more vectors, which may be retroviral or 
not. Is refen-ed to herein as a "vector composition". 

The lA nucleic adds are Introduced Into the cells either alone or in combination with an expression 
vector. By Introduced into " or grammatical equivalents herein is meant that the nucleic adds enter 
the cells in a manner suitable for subsequent expression of the nuddc add. The method of 
introduction Is largely dictated by the targeted cell type, discussed below. Exemplary methods include 
(Ca3P04)2 precipitation, liposome fusion, Upofeclin®, electroporation, viral Infection, etc. The lA 
nucleic adds may stably integrate into the genome of the host cell (for example, with retroviral 
introduction, outlined below), or may exist either transiently or stably in the cytoplasm (i.e. through the 
use of traditional plasmds, utiOzing standard regulatory sequences, selection maricers. etc.). 

The !A proteins of the present Invention are produced by culhiring a host cell transformed witti an 
expression vector containing nudeic add encoding an lA protein, under the appropriate conditions to 
induce or cause expression of the lA protein. The conditions appropriate for lA protein expression will 
vary with the choice of the expression vedor and the host cell, and will be easily ascertained by one 
skilled in the art through routine experimentation. For example, the use of constitutive promoters In 
the expression vector wilt require optinrazing the growth and profiferatton of the host cell, while the use 
of an Indudble promoter requires the appropriate growth conditions for indudlon. In addition, In sonne 
embodiments, the timing of the han^est is important For example, the baculovlral systems used In 
insect cell expression are lytic viruses, and thus han/est time seledlon can be crudal for product yieW. 
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Appropriate host cells Include yeast, bacteria, archebacteria, fungi, and insect and animal cells, 
including mammalian cells. Of particular interest are Dmsophila melangaster ceWs, Saccharomyces 
cerevisiae and other yeasts. E coli, Bacillus stAtilis, SF9 cells, CI 29 cells. 293 cells. Neurospora, 
BHK, CHO, COS. Pichia Pastoris, etc. 

In a preferred embodiment, the I A proteins are expressed In mammalian cells. MammaBan expression . 
systems are also known in the art. and include retroviral systems. A mammaOan promoter is any DMA 
sequence capable of b'nding ntammaGan RNA polymerase and initiating the downstream (3*) 
transcription of a coding sequence for the fusion protein into mRNA. A promoter will have a 
transcription initiating region, which is usually placed proximal to the 5' end of the coding sequence, 
and a TATA box. using a located 25-30 base pairs upstream of the transcription initiation site. The 
TATA box is thought to direct RNA polymerase II to begin RNA synthesis at the con'ect site. A 
mammalian promoter will also contain an upstream promoter element (enhancer element), typically 
located within 100 to 200 base pairs upstream of the TATA box. An upstream promoter element 
deterrnnes the rate at which transcription is initiated and can act in either orientation. Of particular 
use as mammafian promoters are the promoters from mammalian viral genes, since the viral genes 
are often highly expressed and have a broad host range. Examples include the SV40 eariy promoter, 
mouse mammary tumor virus LTR promoter, adenovirus major late promoter, herpes simplex virus 
promoter, and the CMV promoter. 

Typically, transcription termination and polyadenylation sequences recognized by mammalian cells are 
regulatory regions located 3' to the translation stop codon and thus, together with the promoter 
elements, flank the coding sequence. The 3' terminus of the mature mRNA is fomned by site-specific 
post-translational cleavage and polyadenylation. Examples of transcription terminator and 
polyadenlytion signals include those derived fomi SV40. 

The methods of introducing exogenous nucleic add Into mammalian hosts, as well as other hosts. Is 
well known in the art, and wiU vary with the! host cell used. Techniques include dextran-mediated 
transfedion, calcium phosj^ate precipitation, polybrene mediated transfedion, protoplast fusion, 
electroporation, viral infection, encapsulation of the polynucleotide(s) in Oposomes, and direct 
microinjection of the DNA into nuclei. As outlined herein, a particularty prefered method utifizes 
retroviral infection, as outlined in PCT US97/01019, incorporated by reference. 

As will be appreciated by those in the art, the type of manvnalian cells used in the present invention 
can vary widely. Basically, any mammalian cells may be used, with mouse, rat, prinnate and human 
cells being particularty preferred, although as will be appreciated by those in the art, modifications of 
the system by pseudotyping allows all eukaryotic cells to be used, preferably higher eukaryotes. As is 
more fully described below, a screen will be set up such that the cells exhibit a selectable phenotype in 
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the presence of a btoactive peptide. As is more fully described t>elow, cell types impOcated in a wide 
variety of disease conditions are particulariy useful, so long as a suitable screen may be designed to 
allow, the selection of cells that exhibit an altered phenotype as a consequence of the presence of a 
peptide within the cell. 

Accordingly, suitable celt types include, but are not linrvted to, tumor cells of all types (particularly 
melanoma, myeloid leukemia, carcinomas of the lung, breast, ovaries, colon, kidney, prostate, 
pancreas and testes), cardiomyocytes, endothelial cells, epithelial cells, lymphocytes (T-cell and B 
cell), mast cells, eosinophils, vascular intimal cells, hepatocytes, leukocytes including mononuclear 
leukocytes, stem cells such as haemopoetic. neural, skin, lung, kidney, liver and myocyte stem cells 
(for use in screening for differentiation and de-differentiation factors), osteoclasts, chondrocytes and 
other connective tissue cells, keratinocytes, melanocytes, liver cells, kidney cells, and adipocytes. 
Particulariy preferred cells are pancreatic beta cells or variants thereof. Suitable cells also include 
known research cells, including, but not limited to. Jurist T cells. NIH3T3 cells. CHO, Cos, etc. See 
the ATCC ceil fine catalog, hereby expressly incorporated by reference. 

In one embodiment, the cells may be additionally genetically engineered, that is. contain exogeneous 
nucleic add ottier than tiie lA nucleic add. 

In a preferred embodiment, the lA proteins are expressed in bacterial systems. Bacterial expression 
systems are well known in the art. 

A suitable bacterial promoter is any nudeic add sequence capable of binding bacterial RNA 
polymerase and initiating the downstream p*) transcription of the coding sequence of the I A protein 
into mRNA. A bacterial promoter has a transcription initiation region which is usually placed proximal 
to tiie 5* end of tiie coding sequence. This transcription initiation region typically indudes an RIMA 
polymerase binding site and a transcription initiation site. Sequences encoding metabolic pathway 
enzymes provide particularly useful promoter sequences. Examples include promoter sequences 
derived from sugar metabolizing enzymes, such as galactose, lactose and maltose, and sequences 
derived from biosynttietic enzyn^es such as tryptophan. Promoters from bacteriophage nwy also be 
used and are known in the art. In addition, synthetic promoters and hybrid promoters are also useful; 
for example, the fac promoter is a hybrid of the trp and lac promoter sequences. Furthermore, a 
bacterial promoter can indude naturally occurring promoters of non-bacterial origin that have the ability 
to bind bacterial RNA polymerase and initiate transcription. 

In addition to a functioning promoter sequence, an efiident ribosome binding site is desirable. In E 
CO//, the ribosome binding site is called ttie Shine-Delgamo (SD) sequence and indudes an initiation 



56 



wo 00/69901 



PCT/USOO/13764 



codon and a sequence 3-9 nucleotides in length located 3-11 nucleotides upstream of the initiation 
codon. 

The expression vector may also include a signal peptide sequence that provides for secretion of the lA 
protein in bacteria. The signal sequence ty|xca!ly encodes a signal peptide comprised of hydrophobic 
amino adds which direct the secretion of the protein from the cell, as Is well known in the art The 
protein Is either secreted Into the growtii media ^ram-positive bacteria) or into the perlplasntic space, 
located between the inner and outer membrane of the cell ^ranvnegative bacteria). For expression in 
bacteria, usually bacterial secretory leader sequences, operably linked to an lA protein encoding 
nucleic add, are preferred. 

The bacterial expression vector may also include a selectable marker gene to allow for the selection of 
bacterial strains that have been transformed. Suitable selection genes indude genes which render the 
bacteria resistant to dmgs such as ampidIRn, chloramphenicol, erythromydn. kanamydn, neomydn 
and tetracycline. Selectable markers also include biosynthetic genes, such as those in the histidine, 
tryptophan and leudne biosynthetic pathways. 

These components are assembled into expression vectors. Expression vectors for k)acteria are well 
known in the art, and indude vectors for Bacillus subtilis, E. coll, Streptococcus cremoris, and 
Streptococcus IMdans, among others. 

The bacterial expression vectors are transfomied into bacterial host cells using techniques well known 
In the art, such as caldum chloride treatment, electroporation, and others. 

In one embodiment, I A proteins are produced in insect cells. Expression vectors for the 
transfomnation of insect cells, and in particular, baculovirus-based expression vectors, are well known 
in the art. 

In a preferred embodiment, lA proteins are produced in yeast cells. Yeast expression systems are 
well known in the art and indude expression vectors for Saccharomyces cerevisiae, Candida albicans 
and C. mattosa, Hansenula polymorpha, Kluyveromyces fragilis and K. lacGs, Pichia gwllerimondTi and 
P. pastoris, Sdiizosaccharomyces pombe, and Yarwwia lipolytica. Preferred promoter sequences for 
expression in yeast indude the indudble GAL1.10 promoter, the promoters from alcohol 
dehydrogenase, enolase, glucoWnase, glucose-6-phosphate isomerase, glyceraldehyde-3-phosphate- 
dehydrogenase, hexokinase, phosphofructokinase. 3-phosphoglycerate mutase. pyruvate kinase, and 
the add phosphatase gene. Yeast seledable markers Indude ADE2. HiS4, l£U2. TRP1 , and ALG7. 
which confers resistance to tunicamydn; the neomydn phosphotransferase gene, which confers 
resistance to jG418; and the CUP1 gene, which allows yeast to grow in the presence of copper ions. 
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In addition, the lA polypeptides of the invention nnay be further fused to other proteins. If desired, for 
example to increase expression or stabilize the protein. 

In one enfibodiment, the lA nucleic adds, proteins and antibodies of the invention are labeled with a 
label other than the scaffold. By "labeled" herein is nieant that a compound has at least one element, 
isotope or chemicat conpound attached to enable the detection of the compound. In general, labels 
fall Into three classes: a) isotopic labels, which may be radioactive or heavy isotopes; b) immune 
labels, which may be antibodies or antigens; and c) colored or fluorescent dyes. The labels may be 
incorporated into the compound at any position. 

Once made, the lA proteins may be covalently modified. One type of covalent modification includes 
reacting targeted amino add residues of an I A polypeptide vwth an oiganic derivatizing agent that is 
capable of reacting with selected side chains or the N-or G-temninai residues of an lA polypeptide. 
Derivatization with bifunctional agents is useful, for instance, for crosslinking an lA protein to a water- 
insoluble support matrix or surface for use in the method for purifying anti-IA protein antibodies or 
screening assays, as is more fully described below. Commonly used crossOnking agents include, e.g., 
1,1-bis(diazpacetyl>-2-phenylethane, glutaraldehyde, N-hydroxysucdnlmide esters, for example, esters 
with 4-azldosaRcylic add. homobifunctional invdoesters, including disucdnimidyl esters such as 3,3 - 
dithiobis(sucdnirridylpropionate), bifunctional maleimides such as bis^^maleimido-1.8-octane and 
agents such as methyl-34(p-a2idophenyOdithio]propioiniidate. 

Other modifications include deamidation of glutaminyl and asparaginyl residues to the con^esponding 
glutamyl and aspartyl residues, respectively, hydroxylation of proline and lysine, phosphorylation of 
hydroxyl groups of seryl or threonyl residues, methylation of the "-aniino groups of lysine, arginine. and 
histidine side chains |T.E. Creighton, Proteins: Stnicture and Molecular Properties, W.H. Freeman & 
Co.. San Frandsco, pp. 79-86 (1983)1, acetylation of the N-terminal amine, and amidation of any C- 
teminal cart^oxyl group. 

Another type of covalent modification of the I A polypeptide Included within the scope of this invention 
comprises altering the native glycosylation pattern of the polypeptide. "Altering the native glycosylation 
pattem" Is Intended for purposes herein to rr^ean deleting one or more cartjohydrate mdeties found in 
native sequence lA polypeptide, and/or adding one or nwre glycosylation sites that are not present In 
the native sequence lA polypeptide. 

Addition of glycosylation sites to lA polypeptides may be accompGshed by altering the amino add 
sequence thereof. The alteration may be made, for example, by the addition of. or substitution by, one 
or more serine or threonine residues to the native sequence I A polypeptide (for O-linked glycosylation 
sites). The lA protein amino add sequence may optionally be altered through changes at the Dt^lA 
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level, particularly by mutating the DMA encoding the lA polypeptide at preselected bases such that 
codons are generated that will translate Into the desired amino acids. 

Another means of increasing the nunriber of carbohydrate nrK)ieties on the !A polypeptide Is by 
chenical or enzymatic coupling of glycosides to the polypeptide. Such methods are described in the 
art, e.g., in WO 87/05330 pubGshed 11 September 1987. and in Aplin and WHston, CRC Crit. Rev. 
Biochem.. pp. 259-306 (1981). 

Removal of carbohydrate moieties present on the I A polypeptide may be accomplished chemically or 
enzymaticaily or by mutational substitution of codons encoding for amino add residues that serve as 
targets for glycosylation. Chemical deglycosylation techniques are known in the art and described, for 
instance, by HaWmuddin. et aL, Arch. Biochenrv Biophys.. 259:52 (1987) and by Edge et al.. Anal. 
Biochem.. 118:131 (1981). Enzymatic cleavage of carbohydrate moieties on polypeptides can be 
achieved by the use of a variety of endo-and exo-glycosidases as described by Thotakura et al.. Meth. 
EnzymoL. 138:350 (1987). 

Such derivatized moieties way improve the solubility, absorption, permeability across the blood brain 
barrier, bfologlcal half life, and the like. Such nnolelies or modifications of lA polypeptides nray 
alternatively eliminate or attenuate any passible undesirable side effect of the protein and the like. 
Moieties capable of mediating such effects are disctosed. for example, in Remington's Phamnaceutical 
Sciences, 16th ed.. rwiack Publishing Co.. Easton. Pa. (1980). 

Another type of covalent modification of lA proteins comprises linking the lA polypeptide to one of a 
variety of nonproteinaceous polymers, e.g.. polyetiiylene glycol, polypropylene glycol, or 
polyoxyalkylenes. in the manner set forth In U.S. Patent Nos. 4.640,835; 4.496.689; 4.301.144; 
4.670.417; 4.791,192 or 4.179,337. 

lA polypeptides of the present invention nray also be modified in a way to form chinneric nnolecules 
comprising an lA polypeptide fused to another, heterologous polypeptide or annino add sequence. In 
one embodiment, such a chimeric nrwlecule comprises a fusion of an lA polypeptide with a tag 
polypeptide which provides an epitope to which an anti-tag antibody can selectively bind. The epitope 
tag is generally placed at the amino-or carboxyWemtinus of the I A polypeptide. The presence of such 
epitope-tagged fomis of an I A polypeptide can be detected using an antibody against tine tag 
polypeptide. Also, provision of \he epitope tag enables the lA polypeptide to be readily purified by 
affinity purification using an anti-tag antibody or anotiier type of affinity matrix ftat binds to the epitope 
tag. In an altemative embodinrrent. the chimeric molecule may conoprise a fusion of an lA polypeptide 
with an immunoglobulin or a particular region of an immunoglobulin. For a tMvalent form of the 
chimeric molecule, such a fusion could be to the Fc region of an IgG molecule. 



59 



wo 00/69901 



PCTAJSOO/13764 



Various tag polypeptides and their respective antibodies are well known In the art. Examples include 
poly-histidine (poly-his) or poly-histidine-glycine (poly-his-gly) tags; the flu HA tag polypeptide and its 
antilwdy 12CA5 [Re!d et al.. Mol. Ceil. Biol. 82159-2165 (1988)]; the c-myc tag and the 8F9. 3C7, 
6E10, G4, B7 and 9E10 antibodies thereto [Evan et al., Molecular and Cellular Biology. 5:3610-3616 
(1985)]; and the Herpes Simplex vims glycoprotein D (aO) tag and its antibody [Paborsky et al.. 
Protein Engineering. 3{6):547-553 (1 990)]. Other tag polypeptides include the Flag-peptide [Hopp et 
al., BioTechnology 6:1204-1210 (1988)]; the KT3 epitope peptide (Martin et al.. Science 255:192-194 
(1992)1; tubuUn epitope peptide [Skinner et al.. J. Biol. Chem. 266:15163-15166 (1991)]; and the T7 
gene 10 protein peptide tag [Lutz-Freyernuith et al. Proc. Natl. Acad. Sd. U.S.A. 87:6393-6397 
(1990)]. 

In a preferred embodiment, the lA protein Is purified or isolated after expression. lA proteins may be 
isolated or purified in a variety of v^ays known to those skilled in the art depending on what other 
components are present in the sample. Standard purification methods include electrophoretic. 
molecular, immunological and chromatographic techniques, including ion exchange, hydrophobic, 
affinity, and reverse-phase HPLC chromatography, and chromatofocusing. For example, the lA 
protein may be purified using a standard anti-library antibody colunrm. Ultrafiltration and diafillration 
techniques, in conjunction with protein concentration, are also useful. For general guidance in suitable 
purification techniques, see Scopes. R.. Protein Purification, Springer-Vertag. NY (1982). The degree 
of purification necessary will vary depending on the use of the lA protein. In some instances no 
purification Will be necessary. 

Once made, the lA proteins and nucleic adds of the invention find use in a number of applications. In 
a preferred embodiment, the lA proteins are administered to a patient to treat aninsufin^ssoctated 
disorder. 

By Insurm associated disorder", -insufin-dependent disorder-. Insulin responsive disorder, 
-condition" or similar grammatical equivalents thereof, herein is meant a disorder that can be 
ameliorated by the administration of a phamiaceutical composition comprising an insufin or an lA 
protein, including, but not limited to type I diabetes, type 2 diabetes; hyperglycemia; diabetes mellitus; 
famifial hyperprolnsuBnemia; disorder of cartiohydrate metabofism; disorder of skeletal muscle protein 
tumoven and disoiders associated with any of the above-fisted conditions. 

In a preferred embodiinent a therapeutically effective dose of an lA protein Is administered to a patient 
In need of treatment. By "therapeuticalV effective dose" herein Is meant a dose that produces the 
effects for which It is administered. The exact dose win depend on the purpose of the treatment, and 
virill l>e ascertainable by one skilled in the art using known techniques. As is known In the art. 
adjustments for lA protein degradation, systemic versus kxafized deOveiy, and rate of new protease 
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synthesis, as well as the age. body weight, general health, sex, diet, time of adnrvnistration, drug 
interaction and the seventy of the condition nnay be necessary, and will be ascertainable with routine 
experimentation by those skilled in the art. 

A "patient" for the purposes of the present invention includes both hun^ns and other animals, 
particularty mammals, and organisms. Thus the methods are appGcable to both human therapy and 
veterinary appfications. In the preferred embodiment the patient is a mammal, and In the most 
prefen-ed embodiment the patient is human. 

The term "treatment" in the instant Invention is meant to include therapeutic treatment, as well as 
prophylactic, or suppressive measures for the disease or disorder. Thus, for example, In the case of 
type I diabetes, successful administration of an lA protein prior to onset of the disease results in 
"treatment" of the disease. As another example, successful administration of an I A protein after 
clinical manifestation of the disease to combat the symptoms of the disease comprises "treatmenf of 
the disease. Treatment" also encompasses administration of an lA protein after the appearance of 
the disease in order to eradicate the disease. Successful administration of an agent after onset and 
after clinical symptoms have developed, witti possible abatement of clinical symptoms and perhaps 
amelioration of the disease, comprises "treatment" of the disease. 

Those "in need of treatment" include mammals, in particular humans, already having the disease or 
disorder, as well as those prone to having the disease or disorder, including those in which the disease 
or disorder Is to be prevented. Specifically, an "individual at risl^ for type I diabetes or diabetes 
melBtus, herein means 0) an Individual having a blood relative with type 1 diabetes or diabetes 
melDtus; 00 auto-antibody-positive individuals without overt type 1 diabetes (these auto-antibodies 
inckude cytoplasmic islet cell auto-antibodies, insulin antibodies and glutanmc add decarboxylase 
auto-antibodies; (m) individuals with histocompatibility (HLA) type DR3 or DR4DQW8; Ov) individuals 
with glucose abnomiaOties such as a loss of first phase insufin secretion on glucose tolerance tests. 

In another embodmeht, a therapeutically effective dose of an lA protein, an lA gene, or an lA protein 
antibody is adntinistered to a patient having a disease involving inappropriate expression of insulin. A 
"disease involving inappropriate expression of insulin" within the scope of the present invention is 
meant to include diseases or disorders characterized by an overabundance of insufln. This 
overabundance may be due to any cause, including, but not fimlted to, overexpresslon at the 
molecular level, prolonged or accumulated appearance at the site of action, or increased activity of 
insuGn relative to nonmal. Included within this definition are diseases or disorders characterized by a 
reduction of insufin. This reduction may be due to any cause, including, but not Omited to, reduced 
expression at tiie molecular level, shortened or reduced appearance at the site of action, or decreased 
activity of insufin relative to normal. Such an overabundance or reduction of InsuOh can be measured 
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relative to normal expression, appearance, or activity of insuGn according to, but not limited to, the 
assays described and referenced herein. 

The administration of the lA proteins of the present invention, preferably In the form of a sterile 
aqueous solution, can be done in a variety of ways, including, but not Gmited to, orally, subcutanjsousty. 
Intravenously, intranasally, transdennatly, intraperitoneaily, intramuscularly, intrapulmonary, vaginally, 
rectally, or intraocularly. In some instances, for example, in the treatment of wounds, inflammation, or 
multiple sclerosis, the lA protein may be directly appGed as a solution or spray. Depending upon the 
manner of introduction, the phamnaceutical composition may be formulated in a variety of ways. The 
concentration of the therapeutically active lA protein in the formulation may vary from about 0.1 to 100 
weight %. In another prefen-ed err^odiment, the concentration of the lA protein Is in the range of 1- 
1 GO U/nH, with dosages from 0.05, 0.1 . 0.2, and 0.5 U/l^ of body weight being prefen^ed. 

The pharmaceutical compositions of the present invention comprise an I A protein in a form suitable for 
administration to a patient. In the preferred embodiment, the pharmaceutical compositions are in a 
water soluble fomi, such as being present as pharmaceutically acceptable salts, which is meant to 
include both add and base addition salts. "Pharmaceutically acceptable add addition salT refers to 
those salts that retain the biological effectiveness of the free bases and that are not biologically or 
otherwise undesirable, formed with inorganic adds such as hydrochloric add, hydrobromic add, 
sulfuric acid, nitric add. phosphoric add and the Gke, and organic adds such as acetic add, propionic 
add, glyceric add, pyruvic add, oxalic add, maleic add, malonic add, succinic add, fumaric add, 
tartaric add, dtric add, benzoic add, dnnamic add, mandeGc add. methanesulfonic add, 
ethanesulfonic add, p-toluenesulfonic add, saGcyQc add and the fike. Tharmaceutically acceptable 
base addition salts" indude those derived from inorganic bases such as sodium, potassium, Othlum, 
ammonium, caldum, magnesium, iron, zinc, copper, manganese, aluminum salts and the Gke. 
Particularly prefen^ed are the ammonium, potassium, sodium, caldum, and nragnesium salts. Salts 
derived from phannaceutically acceptable organic non-toxic bases indude salts of primary, secondary, 
and tertiary amines, substituted amines Inducing naturally occuning substituted amines, cyclic anrunes 
and basic ion exchange resins, such as isopropylamine, trimethylamine, diethylanrune, triethylanrtine, 
tripropylan^ne, and ethanolamine. 

The phamnaceutical compositions may also indude one or wore of the following: earner proteins such 
as serum albunin; buffers such as NaOAc; fillers such as mscrocrystalGne cellulose, ladtose, com and 
other starches; binding agents; sweeteners and other flavoring agents; coloring agents; and 
polyethylene glycol. Addithfes are well known in the art, and are used in a variety of formulations. 

In addition, in one embodiment, the lA proteins of the present invention are formulated using a process 
for phamnaceutical compositions of reconrA)inant insuGn as described in U.S. Patent Nos. 5,843,886, 
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6,034.054. 4,992.417, 5,506,202. 5,559,094. 5,700,662. 5.618.913. 5.514,646. and 5.514.646. which, 
hereby, are expressly incorporated In its entirety. 

In a further enr\bodiment, the lA proteins are added in a nticellular formulation; see U.S. Patent 
No.5,833,948, heretiy expressly incorporated by reference in its entirety. 

Combinations of pharmaceutical compositions may be administered. Moreover, the compositions may 
be adnrunistered in combination with other therapeutics. 

In one embodiment provided herein, antibodies, including but not Omited to nranodonal and polyclonal 
antibodies, are raised against lA proteins using methods known in the art In a prefenred ennbodiment, 
these anti-IA protein antibodies are used for invnunotherapy. Thus, methods of immunotherapy are 
provided. By "immunotherapy" is meant treatment of an IFN-P related disorders with an antibody 
raised against an I A protein. As used herein, immunotherapy can be passive or active. Passive 
immunotherapy, as defined herein, is the passive transfer of antibody to a recipient (patient). Active 
immunization is the Induction of antibody and/or T-cell responses in a recipient (patient). Induction of 
an immune response can be the consequence of providing the recipient with an lA protein antigen to 
which antibodies are raised. As appreciated by one of ordinary skill in the art, the I A protein antigen 
may be provided by Injecting an 1 A polypeptide against which antibodies are desired to be raised into a 
redpent, or contacting the redpient with an lA protein encoding nudeic add. capable of expressing 
the lA protein antigen, under conditions for expression of the lA protein antigen. 

In another preferred embodiment, a therapeutic compound Is conjugated to an antibody, preferably an 
antMA protein antibody. The therapeutic compound may be a cytotoxic agent. In this method, 
targeting the cytotoxic agent to e.g.. pancreatic tumor tissue or cells, results in a reduction in the 
number of affllded cells, thereby redudng symptoms assodated with cancer, and lA protein related 
disorders. Cytotoxic agents are numerous and varied and indude. but are not limited to, cytotoxic 
dmgs or toxins or active fragments of such toxins. Suitable toxins and their conrespoiidlng fragments 
Indude diptheiia A chain, exotoxin A chain, ridn A chain, abrin A chain, curdn, crotin, phenomydn, 
enomydn and the Oke. Cytotoxic agents also Indude radiochemicals made by conjugating 
radidsotopes to antibodies raised against cell cyde proteins, or binding of a radionucfide to a chelating 
agent that has been covalently attached to the antibody. 

In a preferred embodiment. lA proteins are administered as therapeutic agents, and can be formulated 
as outlined above. Sinularly. !A genes (including both the fulWength sequence, partial sequences, or 
regulatory sequences of the lA protein coding regions) can be administered In gene therapy 
appOcations. as Is known In the art These lA genes can indude antisense appOcations, either as gene 
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therapy (ue. for incorporation into the genome) or as antisense compositions^ as will be appreciated by 
those in the art 

In a preferred embodiment, the nucleic add encoding the lA proteins may also be used in gene 
therapy. In gene therapy applications, genes are introduced into cells in order to achieve in vivo 
synthesis of a therapeutically elfecdve genetic product, for example for replacement of a defective 
gene. "Gene therapy" Includes both conventional gene therapy where a lasting effect is achieved by a 
single treatment, and the adn^nistration of gene therapeutic agents, which involves the one time or 
repeated adntinistration of a therapeutically effective DMA or mRNA. Antisense RNAs and DNAs can 
be used as therapeutic agents for blocking the expression of certain genes in vivo. It has already 
been shown that short antisense oFigonudeotides can be imported into cells where they act as 
inhibitors, despite their low intracellular concentrations caused by tiieir restricted uptake by the cell 
membrane. [Zamecnik et al.. Proc. Natl. Acad. Sd. U.SA 83:4143-4146 (1986)]. The 
oligonudeotides can be modified to enhance their uptake, e.g. by substituting their negatively charged 
phosphodiester groups by uncharged groups. 

There are a variety of techniques available for introdudng nudeic adds into viable cells. The 
techniques vary depending upon whether the nudeic add is transfen-ed into cultured cells in vffro, or in 
vfvD in the cells of the intended host. Techniques suitable for the transfer of nudeic add into 
mammaDan cells in vitro include the use of liposomes, electroporation, nticroinjection. cell ajsion. 
DEAE-dextran, the caldum phosphate predpitation method, etc. The cunently preferred in vivo gene 
transfer techniques indude transfection witti viral (typically retroviral) vectors and viral coat proteln- 
Uposome mediated transfection jDzau et al.. Trends in Biotechnology 11205-210 (1993)1. In some 
situations it is desirable to provide the nudeic add source with an agent that targets the target cells, 
such as an antibody spedfic for a cell surface membrane protein or the target cell, a ligand for a 
receptor on the target cell. etc. Where fiposomes are employed, proteins which bind to a cell surface 
membrane protein assodated with endocytosis may be used for targeting and/or to fedlitate uptake, 
e.g. capsid proteins or fragments thereof tropic for a particular cell type, antibodies for proteins wHch 
undergo intemalization in cycling, proteins that target intracellular k)calizalion and enhance intracellular 
halWife. The technique of receptor-mediated endocytosis is described, for example, by Wu et al., J. 
Biol. Chem. 262:44294432 (1987); and Wagner et al.. Proc. Natl. Acad. Sd. U.SA 87:3410-3414 
(1 990). For review of gene marking and gene therapy protocols see Anderson et al.. Sdence 
256:808-813 (1992). 

In a preferred embodiment, lA genes are administered as DNA vacdnes, either single genes or 
combinations of I A genes. Naked DNA vacdnes are generally known in the art Brower. Nature 
Biotechnology. 16:1304-1305 (1998). Methods for the use of genes as DNA vacdnes are well known 
to one of ordinary skill in the art. and indude pladng an lA gene or portion of an lA gene under the 



64 



wo 00/69901 



PCTAJSOO/13764 



control of a promoter for expression in a patient in need of treatment. The lA gene used for DMA 
vaccines can encode fulWength I A proteins, but more preferably encodes portions of the lA proteins 
including peptides derived from the lA protein. In a preferred embodiment a patient is immunized with 
a DMA vaccine comprising a plurality of nucleotide sequences derived from an lA gene. Similarly, it is 
possible to immunize a patient vwth a plurality of lA genes or portions thereof as defined herein. 
Without being bound by theory, expression of the polypeptide encoded by the DMA vacdne, cytotoxic 
T-cells, helper T-ceils and antibodies are induced which recognize and destroy or eliminate cells 
expressing IFN-P proteins. 

In a preferred embodiment, the DMA vaccines include a gene encoding an adjuvant molecule with the 
DMA vaccine. Such adjuvant molecules include cytokines that increase the Immunogenic response to 
the I A polypeptide encoded by the DMA vacdne. Additional or alternative adjuvants are known to those 
of ordinary skill in the art and find use in the invention. 

The following examples serve to more fully describe the manner of using the above-described 
invention, as well as to set forth the best modes contemplated for carrying out various aspects of tfie 
invention. It is understood that these examples in no way seive to Bnnit the tme scope of this invention, 
but rather are presented for illustrative purposes. All references cited her«n are incorporated.by 
reference in their entirety. 

EXAMPLE 1 

DESIGN AND CHARACTERIZATION OF NOV EL lA PROTEINS BY PDA 

Summary: Sequences for novel insulin activity proteins (lA proteins) were designed by simultaneously 
optin^'zing residues In the buried core of the protein, residues on the surface of the protein and 
residues at tiie boundary of the protein using Protein Design Auton^tion (PDA) as described in 
WO98/47089. U.S.S.Nos. 09/058.459. 09/127,926. 60/104.612. 60/158.700, 09/419.351, 60/181.630. 
60/186.904. and U.S patent appOcation. entitled Protein Design Automation For Protein Ubraries 
(RIed: April 14, 2000; Inventor. Bassil Dahiyat). aU of which are expressly incorporated by reference In 
their entirety. 

Several core designs (6 positions), boundary designs (15 positions) and surfeoe designs (?2 positions) 
were completed, with 43 residues considered in total, confesponding to 20^ sequence possibiOties. 
Residues unatposed to solvent were designed in order to mli^rraze changes to the molecular surface 
and to Bnit the potential for antigenicity of designed novel protein analogues. 
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Calculations required from 12-19 hours on 16 Silicon Graphics R10000 CPU's. The global optimum 
sequence from each design was selected for characterization. From 1-20 residues were changed from 
human inisulin in the designed proteins, out of 51 residues total. 

COMPUTATIONAL PROTOCOLS 

Temolate sfructure preparation: 

For this study the ctystal structure of human Insulin as deposited In the PDB data bank was used [PDB 
record 1TRZ; CIszak and Strtth, Biochemistry 33(6):1512-7 (1994)]. aszak and Smith solved the 
structure by x-ray crystallography to a resolution of 1 .6 Angstrom. The structure of the T3R3 human 
Insulin hexamer is complex ed with two zinc ions. The asymmetric unit consists of a TR dimer and the 
insuOn hexamer Is generated by the action of the crystallographic 3-foW axis. The conformation of one 
InsuGn trimer is neariy identical to that of the TB hexamer, while the other trimer approximates that of 
the R6 hexamer, except for the three N-tenninal B-chain residues that adopt an extended rather than 
an alpha-hefical conformation. Each of the two zinc ions, which lie on the crystallographic 3-fold axis 
and exhibit two different, disordered coordination geometries, is coordinated by the inrudazole groups 
of three symmetry-related BIO-Hs residues. The coordination sphere of the zinc in the T3 trimer is 
either tetrahedral, with the fourth site filled by a chloride Ion. or octahedral, completed by three water 
molecules. The coordination of the zinc in the 12 A narrow channel in the R3 trimer Is tetrahedral, with 
either a second chloride ion or a water molecule completing the coordination sphere [Qszak and 
Snuth, supra). 

PDA calculations were performed for the A-chain and B-chain simultaneously. The zinc Ions and all 
water molecules as well as all hydrogen atoms that are present in the PDB file 1TRZ or 1 wav were 
removed from the stmcture prior to the PDA calculation. 

Desiqri strategies: 

Core residues were selected for design since optimization of these positions can improve stability, 
although stablfization.has been obtained from modifications at other sites as well. Core designs also 
n^nlmize changes to the molecular surface and thus Bmlt the designed protdn's potenfial for 
antigenicity. Other PDA calculations Involved core, boundary and surface residues. 

PDA calculations 

All PDA calculations were perfomied with solvation model 2. Solvation model 2 Is the solvation model 
d^cribed by Street and Mayo [FoW. Design 3253-258 (1998)]. If possible. Dead End Efimination 
(DEE) was run to completion to find the PDA ground slate. The hexamer calculations were performed 
with botii solvation 1 and 2. ghnng similar results. 
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TNs was done for the PDA calpulations for the A-chain and B-chain of insulin In the designs for CO 
disulfide rep{acement{ti) mutations promoting hexamer fomiation, and (iii) global redesigns for 
improved stability. For the trz_08' calculation (see below), DEE was aborted after the rotamer 
sequence space was reduced to less than 10^^ sequences. The DEE calculation was for all the given 
Core calculation followed by l\4onte Carlo (MC) minimization and a list of the 1,000 or up to 10.000 
lowest energy sequences was generated. 

Before the PDA calculations were started an initial preparation of the structure was performed. For the 
A-chain and B-diain, all side chains were nninimlzed with Biograf for 50 steps using conjugate gradient 
procedure without a Coulomb potential, this is followed by an additional 50 steps of conjugate gradient 
nunimization without a Coulomb potential for the complete stnicture of the A-chain using Biograf. This 
minimization procedure was chosen to remove initial bad contacts in the structure. 

The PDA calculations for all the designs were run using the a2h1p0 rotamer library. This library is 
based on the backtjone-dependent rotamer library of Dunbrack and Karplus (Dunbrack and Karplus, J. 
Mol. Biol. 230(2):543-74 (1993); hereby expressly incorporated by reference) but includes more 
rotamers for the aromatic and hydrophobic amino adds; Xi and )(2 angle values of rotamers for all the 
aromatic aoino acids and Xi angle values for all the other hydrophobic anrvno adds were expanded ± 
1 standard deviation about the mean value reported in the Dunbrack and Karplus library. Typical PDA 
parameters were used: the van der Waals scale factor was set to 0.9, the H-bond potential welWepth 
was set to 8.0 kcal/mol. the solvation potential was calculated using type 2 solvation with a nonpolar 
burial energy of 0.048 kcal/mol and a nonpolar exposure multipOcation factor of 1.6, and the secondaiy 
structure scale factor was set to 0.0 (secondary stnjcture propensities were not considered). 
Calculations required from 12-24 hours on 16 SiGcon Graphics R10000 CPU's. 

Monte Carlo analysis 

Monte Carlo analysis of the sequences produced by PDA shows the ground state (optimal) amino add 
and airino adds allowed for each variable position and their firequendes of occun-ence (see Tables 1- 
4). 

EXAMPLE 2 
Disulfide Replacement (Design cvsl) 

THs PDA calculation considers potential cysteine knockouts and all positions within dose proximity 
(SA) to these cysteines to allow for flexibility in the structure and to accommodate additional mutations. 

By visual inspection, the following residues were identified as belonging to the boundary of the protein: 
A5-Gln. A15-Gln. A19-Tyr, B2-Val, B18-Val, and B22-Arg. The following residues were identified as 
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belonging to the core of the protein: A2-lle. A3-Vai. A16-Leu, B11-Leu, B15-Leu, and B24-Phe. The 
following cysteine residues were included in this design: A6; A7, A1 1 . A20. 87, and 819. 

Thus, the following positions were included in the PDA design: 

A2 A3 AS A6 A7 All A15 A16 A19 A20 B2 B7 Bll B15 B18 B19 B22 B24 
lie Val Gin Cys Cys Cys Gin Leu Tyr Cys Val Cys Leu Leu Val Cys Arg Phe 

The boundary and cysteine residues were allowed to change to any amino add, except Cys, Pro, and 
Gly. Core residues were allowed to change to any PHOBIC anino add (Ala, Val, Leu, lie, Phe. Tyr, 
Trp, Met) and the PDA core solvation potential was used including surface area calculation. 

Using Monte Carlo technique a Ust of low energy sequences was generated. The analysis of the 
lowest 1000 protein sequences generated by Monte Carlo leads to the mutation pattern are shown In 
the Table 1: If no probability is shown, a value of 1 was assigned (xlOOO =1 000). 



t 


consensus 


orobabilities (xioOO) 


2 (A2) : 


LEU 


905 












3 (A3): 


ILE 


811 


VAL 


188 


ILE 


811 




5 (A5): 


ARG 


306 


VAL 


84 


GLU 


218 


GLN 


6 (A6) : 


ALA 


999 












7 (A7) : 


SER 


765 


ALA 


234 


SER 


765 




11 (All) ; 


ALA 


999 












15(A15) : 


LEU 


388 


LEU 


388 


GLU 


362 


GLN 


16 (A16) : 


ILB 


566 


LEU 


383 


ILE 


566 




19(A19) : 


PHE 


999 












20(A20) : 


ASP 


538 


ALA 


119 


SER 


332 


ASP 


23 (B2) : 


GLU 


386 


ASP 


70 


TISN 


71 


GLU 


28 (B7): 


TYR 


808 


HSP 


171 


TYR 


808 




32 (Bll) : 


LEU 


968 












36 (BIS) : 


LEU 


999 












39(B18) : 


VAL 


999 












40(B19) : 


ALA 


999 












43(B22) : 


GLN 


985 












45(B24) : 


PHE 


963 













ARG 306 



ARG 52 



GLN 235. LYS 117 ARG 85 



Thus, any protein sequence showing nmitations at the positions according to Table 1 will potentially generate 
a more stable and Mth^e lA protein. In particular those protein sequences found among the list of the lowest 
101 MC generated sequences (data not shown) have a high potential to result in a more stable and active 
lA protein. A prefen-ed lA sequence derived from the PDA design *cysr is shown in Figure 3A. This 
sequence shows 20 mutations when compared to the wild type insufin: G-A1-N. I-A2-L, C-A6-A. C-A7-S, 
l-A10^. C-A1 1-A, L-A16-I. E-A17-Y. Y-A19F. C-A2043. F-B1-D. CHB4-F, C-B7-Y, L-B1 1-1. V.B12-R. A-BM- 
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W. C-B19-A. F-B25-N, Y-B26-F, and T-B27-D, Cysteines at positions A6, All, and B19 all become Ala, 
indicating a lack of space at these positions. 

Any combination of amino acid substitutions shown in Table 1 is possible and has a high probability to result 
in an lA protein that has fewer disulfide bonds than insulin. 

Dinntfid'^ p^pia>^amont fnAgif|nR rysTTa, cysTTb, cys 77d. cvs77d+. helix24 and cvs-4> 

PDA design 'cysTT' , only A7-Cy s, B7-Cys and all residues within 6 A of their sulftjr atoms (A3;7, 82,4,7,1 1) 
were considered in the calculation. 

By visual inspection, the following residues were Identified as belonging to the boundary of the protein: A7- 
Cys, B2-Val, B4-Gln, and B7-Cys. The following residues were identified as belonging to the core of the 
protein: A3-Val and B1 1-Leu. 

Thus, the following positions were included in the PDA design: 

A3 A7 B2 B4 B7 Bll 
Val Cys Val Gin Cys Lou 

Using Dead-End EUn^nation and Monte Cario technique a list of low energy sequences was generated. 
The analysis of the lowest 1000 protein sequences generated by Monte Carto leads to the mutation pattem 
are shown in the Table 2: 



£ consensus 

3 <A3) VAL 751 

7 (A7) SER 573 

23 (82) GLN 204 

25 (B4) TYR 304 

28 (B7) TYR 467 

32 (Bll) LEU 756 



probabilities (xlOOO) 

ALA 161 VAL 751 ILE 80 

ALA 408 SER 573 

THR 90 VAL 54 

GLN 204 LYS 78 

PHE 207 TYR 304 

ALA 139 SER 107 

LEU 756 ILE 239 



ASP 76 
ARG 103 
GLU 67 GLN 93 
ASP 157 ASN 52 



ASN 51 GLU 185 



LYS 62 
HSP 77 



ARG 154 
TYR 467 



Tluis/any protein sequence Showing mutations at the positions according to Table2 Will po 

a more stable and active lA protein. In particularthose protein sequences found among the list of the lowest 

101 MC generated sequences (data not shown) have a high potential to result in a more stable and active 
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lA protein. A preferred lA sequence derived from the PDA design 'cysyTa* is shown in Figure 3B. This 
sequence shows 4 mutations when cofttpar&l to the wikJ type insulin: C- A7-S , V-B2-E, Q-B4-Y, and C^7-Y. 

PDA design •cys77b' is sintilarto PDA design •cys77a\ however, in this calculation only positions A7 and 
87 were allowed to change to other residues. The other positions had their amino acid identities fixed, 
but were allowed to change conformation. X)ys77b \vas a design of only A7-B7, the reduced conformational 
freedom bloclced the Tyr mutation from B7. A preferred sequence from this design is shown In Figure 3C. 
This sequence shows 2 mutations, C-A7-S and C-B7-D. 

PDA designs 'cys77tf and •cys77d+ * are the minimal sets of mutations that allow the A7 Ser, B7 Tyr pair 
to occur. Preferred sequences from these two designs are shown in Figures 3D and 3E. 

The sequence shown in Figure 3D shows 4 mutations. C-A7-S , V-B2-T. Q-B4^Y. and C-B7-Y. 

The sequence shown In Figure 3E shows 3 mutations, C-A7-S , Q-B4-Y. and C-B7-Y. 

PDA design 'helix 24* is same as PDA design •cys77a* except that the backbone of B-chain residues 1-5 
were set to helical ui of -57" and 50°, respectively. This is an idealized R state w'lth more room. 

A prefen-edsequenceobtainedfromthisdeslgnisshowninFigureSF. This sequence shows 2 mutations. 
OA7-S and OB7-E, 

PDA design 'cys4* included deletion of residues B1-B4 and calculations forppsftipns A3, A7, B3, and B7. 

ApreferredsequenceobtainedframthisdesignisshowninFigure3G. This sequence shows 2 substitutions. 
OA7-E and OB7-E and deletions of residues at positions B1 to B4. This lA protein generates nroreroom 
forA7-B7. 

Other preferred mutations fbrthe A-chain are as follows: lAla. 2l^u, 6Ala, 7Ser. 7Ala, 7Glu, 7Asn, 10Gln. 
1 lAla. leile. 17Tyr, 19Phe, 20Asp, 20Ser, and 20Ala. 

Other preferred mutations for the B-chain are as follows: 1 Asp. 2Glu. 2Thr. 4Tyr. 4Phe. 7Tyr. 7His. 7Asp. 
7Aia. 7Ser. 7G!u. 1 llle, 12Arg. 14Trp, lOAIa. 25Asn. 26Phe, and 27Asp. 

Any combination of arrtno ackl substitutions shown In Table 2 is possible and should result in a more stable 
lA protein. 
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Generally, the steric constraints in the insulin dimer make the disulfide positions difficult to model for side 
chain suljstitution without perturbing the backbone significantly. A6^11 was very sterically constrained 
as was A20<B1 9. The disulfide A7-B7 was an exception, and was replaced with a Ser-Tyr. or Ser-Glu/Asp 
as shown above for various calculations. 



EXAMPLE S 

Insulin Variants that Promote Hexamer Formation: B14 Substitutions 

In preliminary calculations, position 14 on the B chain optimized to become a large hydrophobic residue, 
despite its location near the surface of the protein. In the R6 hexamer this position is nrrastiy buried, and 
is in contact with the phenol molecule that stabiKzes the R-form of the hexamer. In an effort to occlude 
the phenol binding site of the hexameric R6 confonnation of human Insulin (PDB Iwav; Ding et al. Sd. 
China C. Life Sd. 39(2):144-53 (1996); hereby expressly incorporated as reference), orto mimic the allosteric 
effect of phenol binding, redesigns on the residues (from different subunits) sun-ounding one of the phenol 
binding sites (in the absence of phenol) were made. Predicted stmctures of these mutants fill the space 
occupied by the phenol in the hexamer, and mimic the structural placement of the phenol and its cherrtcal 
identity. By ocduding the phenol binding site, phenol binding will be prevented and the link between phenol 
binding and R state formation will be broken. 

All of the PDA Vrav14' design sequences and hexamer sequences were derived from calculations that 
were done on the entire insuOn hexamer complex (PDB file Iwav; Ding et al., supra) which contains 6 sets 
ofinsufinheterodimersCchainsI through 12). Aslngle phenol binding site centered on the chain 1 2. Ala-BI 4 
(theethinsuRnheterodimerinthecomplex.B-chaln) wasdesigned. Thefollowing positions were involved: 

where the first number Is the PDA renumbering of residues and the second number is the chain number 
in the hexamer; the core and boundary residues are spedfied last 



140 


#6 


LeuB17 


boundary 


176 


#8 


ValB2 


boundary 


179 


#8 


HsB5 


boundary 


180 


#8 


LeuBG 


boundary 


271 


#11 


LeuAie 


core 


287 


#12 


UuBII 


core 


290 


#12 


AlaB14 


boundary 



These positions were chosen because they are within 4.6 A from the B14 side chain. 
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Using DEE and Monte Carlo technique a.list of low energy sequences was generated. The analysis of 
the lowest 1000 protein sequences generated t3y Monte Carlo leads to the nujtation pattern are shown In 



the Table 3: 










res# 




consensus 


orobabinties (x1000) 


140 


#6 


LeuB17 


Leu 350 


Val184,Leu 597,GIu218 


176 


#8 


ValB2 


Val587 


Val684. Phe183. Glu133 


179 


#8 


HisB5 


Phe490 


Leu 205, Phe 597,Trp199 


180 


#8 


LeuB6 


Leu 767 


Leu 1000 


271 


#11 


LeuA16 


Phe667 


Leu 108. Phe 713. Tyr 179 


287 


#12 


LeuBII 


lie 514 


Leu 449, lie 551 


290 


#12 


AlaB14 


Trp469 


Leu 157. Phe 228. Trp 615 



Preferred lA sequences obtained fronra this PDA calculation and comprising only B14 substitutions are shown 
in Figures 4C. 4D. 4F, and 4G. 

Further, substitutions B14Trp or B14Phe will prevent formation of the T state formed by the monomer in 
solution, further driving the equiObrium to the R state. This will promote self-association of insulin, since 
the R state is not stable as a monomer. 

Other substitutions at B14 with similar effects are B14Tyr (Figure 4F) and B14 lie (Figure 4G). These 
mutations may require further mutations as described further below for B5 substitutions. 

InsuGn Variants that Promote Hexamer Formation: B5 Substitutions 

In this PDA design the B5 position of the B-chain was analyzed with respect to possible substitutions that 
promote insuGn hexamer formation. 

Substitutions of HisBS to Trp or Phe is expected to have the same occlusion of the phenol Wnding site in 
h«camerjc insulin as a mutation of B14 to Trp or Phe. These B5 substitutions should destabilize the 
monomeric form of insufin (relative to the unfolded monomer) because the side chain of B5 is solvent 
exposed in the monomer and should stal^Oze the R6 fomt of the hexamer (relative to the wild type fonn 
of the R6 hexamer in the absence of phenoO because it has fomved a close interaction with other nonpolar 
side chains from other subunits wWch surround the phenol binding site, similar to the Interactions normally 
madebyphenoL 

Unlike the B14 mutations, the presence of Trp and Phe at B5 is sterically compatible with a T state. 
A preferred lA protein sequence comprising a B5 substitution only is shown In Rgure 4B. 



72 



wo 00/69901 



PCT/ySOO/13764 



Insufin Variants that Promote Hexamer Formation: B14/B5 Double Substitutions 

In this PDA design the B5 and B14 positions of the B-chain were analyzed with respect to possible 
substitutions that promote insulin hexamer formation. The following substitutions were found: 
B-chain. position 14: Phe. Trp. Tyr, or lie; and 
B-chain. position 5: Phe or Trp. 

Any combination of the above substitutions is possible, such as B14-Phe/B5-Phe, B14-Phe/B5-Trp. B14- 
Trp/B5-Phe, B14-Trp/B5-Trp. B14-Tyr/B5-Phe, B14-Tyr/B5-Trp, B14-lle/B5-Phe. and B14-lle/B5-Trp. 
Preferred lA sequences comprising B5 and B14 substitutions are shown in Figure 4A and in Figure 4E. 

Insulin Variants that Promote Hexamer Formation: Summary 

The results suggest that His-B5-Phe and/or Ala-BI 4-Phe should occlude the phenol binding site and prevent 
binding of phenol-fike molecules at these sites in the hexamer. B14-Phe should be more effective at 
preventing phenol binding simply for steric reasons. Bther mutation can be made singly, or they can be 
made as a double mutant, which should be even more effective at occluding the phenol binding site. 

B5-Phe is completely compatible with the T-form of the monomer, but B14-Phe (and the double mutant 
containing BUPhe) should be sterically unable to form the standard T-state. in either the monomer orthe 
hexamer, due to steric clashes with atoms from the first several residues of the B-chain, particularly B6-Leu. 

Thus, in addition to preventing phenol binding, these mutations should promote the formation of the R6 
form of the protein ^hereby mimicking phenol bending) for two distinct reasons. Rrst, the BSPhe mutation 
should destatMlize the monomeric form of the protein (relative to the unfolded monomer) because the side 
chain of B5 is solvent exposed in the monomer and should stabilize the R6 form of the hexamer (relative 
to the wild type form of the R6 hexamer in the absence of phenol) because it has formed a dose interaction 
with other nonpolar side chans from otfier subunits which surround the phenol binding site, sinwiar to the 
interactions normally made by phenol. Second, the B14-Phe mutation is sterically incapable of forming 
the nonrol T-state monomer structure, thereby providing a large driving force for formation of the R6 
hexamer. 

Thus, all the mutants described above under B14 substitutions, B5 substitutions and B14/B5 double 
substitutions should effect the stablttly of the monomer and the nronomenhexamer equflibrium as wen as 
the T6 ^ T3R3 ^ R6 equilibria within the hexamer (although the T state should not be the canonical T 
state as seen In structural studies of the wild type nunomei). 
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EXAMPLE 4 
Global Redesigns for Improved Stabilttv 

In these PDA designs (trz*_06\ tr2_07b* and trz_08'), the insulin monomer was improved for stability. 
No consid eration was given to the possible functional role of the mutations. For example. Gly-AI is always 
designed to be A1*Asn, though this residue should be more conserved for function. 

In general, PDA design ni2j)6' is nrore conservative than PDA designs 1(z_07b* and trzjDS*. In all three 
PDA calculations, residues at the following positions were considered 'core* residues: A2, A3. A1 6, B1 1 , 
B15, and B24; residues at the following positions were considered 'boundary' residues: A5. A15. A17, A19. 
A21. B2, B3. B4. B8. B12, B14, 818. B22, B26. and B28; residues atthefollowing positions were considered 
'surface' residues: A1. A4. A8. A9. AID. A12, A13, A14. A18. B1. B5. B6. B9, B10. B13. B16. B17. B21, 
B25. B27. B29. and B30. 

A prefeired lA protein sequence from the PDA calculation tiz_06' is shown in Figure 6A, This sequence 
shows 6 mutations when compared to the wild type insulin sequence, G-A1-N, I-A10-Q, L-A16-Y, F-B1-D, 
F-B25-N, and T-B27-D. 

A preferred lA protein sequence from the PDA calculation trz_07b' is shown in Figure 5B. This sequence 
shows UmutationswhenconrparBdtothewildtypeinsulinsequence, G-A1-N, I-A10-Q. L-A16-Y,E-A17-Y. 
Y-A19-F. F-B1-D. V.B2-K. Q-B4-F. L-B11-I. V-B12-R. A-B14-W, F-B25-N. Y.B26-F and T-B27-D. 

A preferred lA protein sequence from the PDA calculation trz_08' Is shown in Figure 5C. This sequence 
shows lemutations when compared to thewUd type insulin sequence, G-A1-N. I-A10-Q. L-A16-Y, E-A17-Y. 
Y-A19-F. F-B1-D, V-B2-K. OrBA-F, G-B8-L, L-B1 1-1, V-B12-R, A-B14-W, F-B25-N. Y-B26-F. T-B27-D and 
P-B28-N. 

As note In the PDA design to promote hexamer binding (see Example 3), Tip at position B14 replaces Ala 
in PDA designs tfz_07b* and trz_08'. Ttiere appears to be a small cleft between the B-chain helix and 
the A-chain strand that accommodates Trp binding. This Is the same ddt used In the hexamer design 
in Example 3. 

Other preferred mutations for the A-chain are as follows: 1 Asn. lOGIn, 16Tyr, 17Tyr, 17Lys, 17Trp, and 
19Phe. 
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Other preferred mutations fortheB^:hain are as foltows: 1Asp,2Lys,4Tyr,4Phe. 8Leu.8Lys. 12Arg,14Trp, 
14GIU. 18Lys. 25Asn, 26Phe. 27Asp, 28Asn and 28Phe. 



Using Monte 
1000 protein 



Carlo technique a list of low energy sequences was generated. The analysts of the lowest 
sequences generated tsy Monte Carlo leads to the nrtutation pattern are shown in the Table 



4: 
i 




consensus 


orobabilitv table ( xlOOO) 


1 


(Al) 


ASH 993 


ASN 1000 






2 


(A2) 


ILE 992 


ILE 1000 






3 


(A3) 


ILE 819 


VAL 180 


ILE 820 




4 


(A4) 


GLN 874 


GLN 1000 






5 


(AS) 


GLU 376 


GLU 383 


GLN329 


ARG 287 


8 


(A8) 


GLN 638 


ASP 154 


GLN 846 




9 


(A9) 


ASM 858 


ASN 879 


GLN 121 




10 


(AlO) 


GUI 921 


GLN 1000 






12 


(A12) 


LYS 415 


THR 226 


GLN 332 


LYS 441 


13 


(A13) 


GLU 855 


GLU 1000 






14 


(A14) 


LYS 778 


LYS 846 


ARG 154 




15 


(A15) 


GLU 324 


LEU 216 


GLU 346 


GLN 132 


16 


(A16) 


TYR 979 


TYR 1000 






17 


(A17) 


LYS 399 


TRP 410 


GLN 128 


. LYS 463 


18 


(A18) 


SER 441 


SER 534 


GLU 213 


LYS 130 


19 


(A19) 


PHE 537 


PHE 543 


TYR 457 




21 


(A21) 


ASN 412 


ASP 116 


ASN 470 


GLN 123 


22 


(Bl) 


ASP 929 


ASP 1.0 






23 


(B2) 


LYS 93 0 


LYS 1000 






24 


(B3) 


ASN 998 


ASN 1000 






25 


(B4) 


PHE 672 


PHE 673 


TYR 327 




26 


(B5) 


GLU 610 


GLU 646 


LYS 164 


ARG 190 


27 


(B6) 


GLU 309 


GLU 370 


GLN 266 


LYS 207 


29 


{B8) 


LYS 603 


LEU 245 


GLU 128 


LYS 627 


30 


(B9) 


ARG 593 


GLU 188 


GLN 169 


ARG 643 


31 


(BIO) 


ARG 421 


GLU 256 


GLN 185 


ARG 559 


32 


(BID 


ILE 674 


LEU 322 


ILE 678 




33 


(B12) 


ARG 743 


LYS 219 


ARG 781 




34 


(B13) 


GLU 607 


ASN 368 


GLU 632 




35 


(B14) 


GLU 898 


TRP 102 


GLU 898 




36 


(B15) 


LEU 1000 


LEU 1000 






37 


(B16) 


GLU 775 


GLU 880 


ARG 120 




38 


(B17) 


LYS 728 


LYS 1000 






39 


(B18) 


VAL 652 


VAL 683 


LYS 317 




42 


(B21) 


ARG 623 


GLN 123 


LYS 160 


ARG 717 


43 


(B22) 


GLN 877 


GLN 1000 






45 


(B24) 


TYR 975 


TYR 1000 






46 


(B25) 


ASN 989 


ASN 1000 






47 


(B26) 


PHE 1000 


PHE 1000 






48 


(B27) 


ASP 915 


ASP iooo. 






49 


(B28) 


PHE 873 


ASN 123 


PHE 877 




50 


(B29) 


THR 310 


SER 183 


THR 351 


ASN 219 


51 


(B30) 


SER 606 


ALA 251 


SER 749 





ARG 306 

ARG 124 
ARG 291 



ARG 157 



GLU 247 
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Positions A6. A7, A1 1 , A20, B7, and B19 are cysteines in the wild type sequence and are not included In 
the *trz' series of PDA designs. Positbns B20 and 823 are glycines are were also not included in these 
calculations. Positions that are critical for function can be constrained to wild type and combined with 
mutations from this table to give more stable lA proteins. 

EXAMPLE S 

Measurable Properties that will reflect the Effect of Mutations In lA Proteins: 
The effect of the 814 and/or 85 substitutions can be measured by one or more of the following methods. 

1. Metal binding. 

Zn and Co binding geometry in the hexamer is measured by UV spetnophotometry of the Co species. The 
stabilization of the R state causes a well-characterized spectral shift that depends on whether the 
conformation of insulin within the hexamer is in the R or T state [Nakagawa and Tager. Biochemistry 
31(12):3204-3214 (1992); Braderetal., Biochemistry 30(27):6636-45 (1991)]. 

2. Association. 

The increased tendency to form a hexamer at a given insulin concentration is measured by methods known 
in the ait, such as analytical ultracentrifugation. size exclusion chromatography, peptide circulardichroism. 
and NMR spectrometry. 

3. Confomnatlon. 

The increased Rstateconfonmation is detected using known methods in the art, such as circulardichroism, 
far UV measurement of increased helix content in the B-chain. and Trp fluorimetry. Stmctural changes 
to the R state can be detected by X-ray crystallography or NMER solution of the stnicture. 

4. StabUity. 

Stability is measured using methods known In the art, such as isothermal titratton calorimetry (e.g., Zn 
binding), themnal melts, resistenceto chentical Inactivatbn. and guanidiniumdenaturatk}n. 

The increased stabiRty of the hexamer form of these mutant insulins should provide tonger shelf life and 
resistance to inactivation from environmental stress such as heat, light, or chemical denaturatton. These 
effects will be measurable by accelerated shelf life testing assessed by e.g. , chromatography and by chemical 
denaturatton assessed by e.g.« circular dichroism or fluorescence. 
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A stable kisulin hexamerwithout phenol would also have useful pharmacological properties, such as delayed 
onset of action following adnrunistration subcutaneousty or to the circulation. 

the above described mutants should slow the dissociation ofthehexamer in the bloodstream. Additionally, 
receptor affinity may be reduced. 

5 Phanmacoklnetlcs can be measured by injection of an lA protein into animals or humans followed by a 
glucose challenge at varying tinnes afterward. Delayed onset of glucose blood level control would be 
observed. Receptoraffinitycanbenrieasuredwlthanadipocytetipogenesisassayor/ni^ 

Standard Potencv Assavs: 

(1) Insulin radioreceptor assays are performed, in which the relative potency of an insulin is defined as 
1 0 the ratio of insulin to I A protein required to displace 50% od ^^l-lnsulin specifically bound to insulin receptors 

present on cell mennbranes, e.g.. rat liver plasma membrane fraction. 

(2) Upogenesis assays are performed, e.g., with rat adipocytes, in which relative insulin potency is defined 
as the ratio of insulin to lA protein required to achieve 50% of the maximum conversion of pH]glucose into 
organio-extractable material (i.e. lipids);. 

15 . (3) Glucose oxidation assays are performed in isolated fat cells in which the relative potency of the lA proteins 
is defined as the ratio of insulin to lA protein to achieve 50% of the maximum conversion of gtuoose-1-[^^C] 
into ^^COa. 

(4) Insulin radioimmunoassays are perfonmed, which can determine the immunogenedty of lA proteins 
by measuring the effectiveness by which insulin or an lA protein competes with 125t-lnsulin in binding to 

2 0 specific anti-insulin antibodies. 

(5) Otherassays are performed that measure the binding of insulin oran lA protein to cells known to possess 
specific insulin receptors. 
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CLAIMS 

We claim: 

1 . A norvnaturally occurring lA protein comprising an amino add sequence which oonprises substitution 
of at least one amino acid residue when compared to an amino acid sequence of a naturally occurring human 
insulin and wherein said lA protein has an altered property when compared to the same property of human 
insulin and binds to a cell comprising an insulin receptor. 

2. A non-naturaliy occuning lA protein according to claim 1, wherein said lA protein comprises a 
substitution at a position selected from the group consisting of positions A3. AS, A6, A7, A1 1. A15. A16, 
A19. A20. B2. B7. B15. 819. and 822. 

3. A non-naturally occurring lA protein according to claim 2. wherein substitution is selected from 
the group of A7-S. A7-E. B2-E, 82-T. B4-Y. 87-Y. BA-F, 87-Y, B7-E, and B7-D. 

4. A norvnaturally occurring lA protein according to claim 1, wherein said lA protein comprises 
substitutions of at least four amino acid residues. 

5. A non-naturally occuning lAproteinconformerhavingathreedimensionalbackbonestructurethat 
substantially conresponds to the three dimensional backbone structure of human insulin, wherein the amino 
add sequence of said confonner and said amino acid sequence of said human Insulin are less than about 
98% klentical. 

6. A non-naturally occurring lA protein comprising at least one amino acid substitution as compared 
to human insulin, wherein at least one of said substitutions is selected from the amino add residues at 
positions selected from positions B5 and B14, and wherein said lA protein forms a hexamer in the absence 
of a phenolic preservative. ■ 

7. The non-naturally occurring lA protein according to daim 6. wherein said lA protein comprises 
a substitution selected from the group of B5-F. B5-W, B14-F. B14-W, B14-Y. and B14-I. 

8. The non-naturally occurring lA protein according to claim 1 , wherein said lA protein comprises At 
least 5 substitution at positions selected from the group consisting of positions A1, A10, A16, A17. A19, 
B1 . B2. B4, B8. B1 1 . 812. B14. B25. 826. 827 and 828. 
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9. The non-naturaily occurring lA protein according to claim 8, wherein said substitutions are selected 
from the group of substitutions consisting of A1-N, A10-Q, A16-Y, A17-Y. A19-F. Bl-p. B2-K, B4-F, B8-L, 
Bll-t, B12-R: B14-W, B25-N. B26-F. B27-D, and B28-N. 

1 0. The non-natura!ly occurring lA protein according to claim 1 wherein said I A protein comprises an 
amino acid sequence selected from the group of anino acid sequences shown in Figure 3A, Figure 38, 
Figure SC. Figure 3D, Figure 3E, Rgure 3F, Figure 3G. Figure 4A, Figure 48, Figure4C. Figure 4D, Figure 
4E, Figure 4F, Rgure 4G, Figure 5A. Figure SB. and Figure SC. 

11. A recombinant nucleic add encoding the non-naturally occurring lA protein of claim 1 or 10. 

12. An expression vector comprising the recombinant nucleic acki of claim 1 1 . 

1 3. A host ceil comprising the recombinant nucleic acid of claim 1 1 . 

14. A host cell comprising the expression vector of claim 12. 

15. A method of producing a non-naturally occurring lA protein comprising culturing the host cell of 
claim 13 under conditions suitable for expression of said nucleic add. 

16. The method according to claim 15 further comprising recovering said lA protein. 

17 A pharmaceutical composition comprising an lA protein according to daim 1 or 10 and a 
phanmaceutical carrier. 

18. A method for treating an insulin responsive condition comprising administering an lA protein 
according to claim 1 or 10 to a patient in need thereof. 

1 9. The method according to daim 18, wherein said condition is a disorder of carbohydrate metabolism. 

20. The method according to daim 18. wherein said condition is type 1 diabetes. 

21. The method according to daim 22, wherein said condition is type 2 diabetes . 
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A) . 

1 M2\LVMRLLPL lALLALWGPD PAA7VFVN6HL CGSHLVEALY LVCGERGFFY 
51 TPKTRREAED LQVGQVELGG GPGJIGSLQPL ALEGSLGKRG IVEQCCTSIC 
101 SLYQLENYCN 

B) . 

A-chain B-chain 

1 

1 GIVEQCCTSI CSLYQLENYC NFVNQHLCGS HLVEALYLVC GERGFFYTPK T 

C) . 

ITRZzA 

1 GIVEQCCTSI CSLYQLENYC N 
HHHHTTS HHHHGGG 

1TRZ:B 

i FVNQHLCGSH LVEALYLVCG ERGFFYTPKT 
S HHH HHHHHHHHHG GG 

1TRZ:C 

1 GIVEQCCTSI CSLYQLENYC N 
HHHHSSS HHHHGGG 

17KZ:D 

1 FVNQHLCGSH LVEALYLVCG ERGFFYTPKT 
HHHHHHH HHHHHHHHHG GG S 
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FIGURE 2B 
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A) . A-chain B-chain 

1 . . 

1 SLVEQftSTSfi ASLYQXXNTO M&VI^LIGS HroEgLYLVA GERGFHFDPK T 

B) . A-chain B-chain 

^ 1 

1 GIVEQCSTSI CSLYQLENYC NraNYHLYGS HLVEALYLVC GERGFFYTPK T 

C) . A-chain B-chain 

^ 1 . 

1 GIVEQCSTSI CSLYQLENYC NFVNQHI^GS HLVEALYLVC GERGFFYTPK T 

D)» A-chain B-chain 

„ . . 

1 GIVEQCSTSI CSLYQLENYC NFTNYHLYGS HLVEALYLVC GERGFFYTPK T 

£)• A-chain B-chain 

„„ 1 

1 GIVEQCSTSI CSLYQLENYC NFVNTHLYGS HLVEALYLVC GERGFFYTPK T 

F) . A-chain B-chain 

_ „| . . 

1 GIVEQCSTSI CSLYQLENYC NFVNQHI£GS HLVEALYLVC GERGFFYTPK T 

G) . A-chain B-chain 

1 . 

1 GIVEQCETSI CSLYQLENYC N HlfiGS HLVEALYLVC GERGFFYTPK T 
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A) . A-chain B-chain 

, ^ 

1 GIVEQCCTSI CSLYQLENYC NFVNQPLCGS HLVEPLYLVC GERGFFYTPK T 

B) . A-chain B-chain 

, 

1 GIVEQCCTSI CSLYQLENYC NFVNQFLCGS HLVEALYLVC GERGFFYTPK T 

C) . A-chain B-chain 

1 

1 GIVEQCCTSI CSLYQLENYC NFVNQHLCGS HLVEPLYLVC GERGFFYTPK T 



D). A-chain B-chain 

1 

1 GIVEQCCTSI CSLYQLENYC NFVNQHLCGS HLVEWLYLVC GERGFFYTPK T 



E). A-chain B-chain 

. 1 . „ 

1 GIVEQCCTSI CSLYQLENYC NFVNQPLCGS HLVEWLYLVC GERGFFYTPK T 



F)* A-chain B-chain 

1 . 

1 GIVEQCCTSI CSLYQLENYC NFVNQHLCGS HLVEYLYLVC GERGFFYTPK T 



G). A-chain B-chain 

1 

1 GIVEQCCTSI CSLYQLENYC NFVNQHLCGS HLVEILYLVC GERGFFYTPK T 
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A). A-chain B-chain 

_ 1 . . 

1 MIVEQCCTSfi CSLYQIENYC NDVNQHLCGS HLVEALYLVC GERGFNYDPK T 



B). A-chain ^ B-chain 

1 NIVEQCCTS2 CSLYQYXNFC NDKNFHLCGS HIREWLYLVC GERG FNFD PK T 



C). A-chain B-chain 

.| ^ 

1 NIVEQCCTSQ CSLYQYYNFC NDKNFHLCLS HIREWLYLVC GERG FNFPNK T 
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FIGURE 6 
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